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Executive Summary 


The objective of this study is to investigate the yam geometry of multiaxial preforms. The 
importance of multiaxial preforms for structural composites is well recognized by the 
industry but to exploit their full potential, engineering design rules must be established. 

This study is a step in that direction. In this work the preform geometry for knitted and 
braided preforms were studied by making a range of well designed samples and studying 
them by photo microscopy . The structural geometry of the preforms is related to the 
processing parameters. Based on solid modeling and B-spline methodology a software 
package is developed . This computer code enables real time structural representations of 
complex fiber architecture based on the rule of preform manufacturing. The code has the 
capability of zooming and section plotting. These capabilities provide a powerful means to 
study the effect of processing variables on the preform geometry. The code also can be 
extended to an auto mesh generator for downstream structural analysis using finite element 
method. This report is organized into six section. In the first section the scope and 
background of this work is elaborated. In Section two the unit cell geometries of braided 
and multi-axial warp knitted preforms is discussed. The theoretical frame work of yam path 
modeling and solid modeling is presented in Section three. The thin section microscopy 
carried out to observe the structural geometry of the preforms is the subject in section four. 
The structural geometry is related to the processing parameters in section five. Section six 
documents the implementation of the modeling techniques into the computer code MP- 
CAD. A user manual for the software is also presented here. The source codes and 
published papers are listed in the Appendices. 



Chapter - 1. Introduction 


The important role of preforming in the chain of composite manufacturing processes has 
been well recognized by the composite industry in the recent years. It has been 
demonstrated in many cases that the mechanical properties of a composite can be improved 
by the pre-orientation, pre-shaping and pre-placement of matrices (as in the case of 
commingled thermoplastic composites). The reduction in processing steps due to 
preforming also contributes to the reduction of manufacturing cost of composites. The 
recent rise in popularity of composite production by resin transfer molding (RTM) further 
expand the need for advanced preforming technology. Of the large family of textile 
preforms, the class of multiaxial preforms presents perhaps the most promising solution to 
the economic manufacturing of high damage tolerant structural composites. 

Considering the importance of multiaxial preforms for structural composites, it was 
recognized by NASA and Drexel that engineering design rules must be established in order 
to exploit preforms to their full potential for composites. In order for them to be useful, 
these design rules must be capable of linking preform processing parameters to the 
structural engineering design environment. 

Central to this linkage is the quantification of the structural geometry or the fiber 
architecture of the preforms in terms of processing parameters. This geometric model 
provides a means for the incorporation of materials properties into a preprocessor for down 
stream finite element structural analysis. The fiber geometry for knitted and braided 
preforms were studied by making a range of well designed samples and studying them by 
photo microscopy. 

Based on solid modeling and B -spline methodology a software package was developed at 
the Fibrous Materials Research Center. This computer code enables real time structural 
representations of complex fiber architecture based on the rule of preform manufacturing. 
The code has the capability of zooming and section plotting. These capabilities provide a 
powerful means to study the effect ot processing variables on the preform geometry. The 
code also can be extended to an auto mesh generator for structural analysis. 
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Chapter 2. - Processing Model 


2.1 Background 

The current trend in the composite materials industry is to expand the use of composites 
from secondary non-load bearing to primary load bearing structural application. This 
requires a significant improvement of the damage tolerance and reliability of composites. 

In addition, it is also desirable to reduce the cost and broaden the usage of composites from 
aerospace to automotive and building construction applications. This calls for the 
development of a capability for quantity production and the net or near-net shape 
reinforcement of structural composites. 

In order to improve the damage tolerance and delamination property of composites, a high 
level of through-the-thickness strength is required. The reliability of a composites depends 
on the uniform distribution of the materials and consistency of interfacial properties. The 
structural integrity, handleability and formability of the reinforcing material for the 
composite are critical for large scale automated production of structural shapes. 

As introduced by Ko [2.1], there is a large family of textile preforms available for advanced 
composites ranging from 3-D integrated net-shape structures to thin and medium thickness 
multilayer fabrics. In addition to properties translation efficiency, structural integrity and 
formability, a key requirement for textile preforms for building construction is their 
availability at a reasonable fabrication cost. There are two families of preforms which have 
the potential to meet the demand for structural composites 3-D braid and multiaxial war 
knit (MWK) fabrics. 

In order to fully understand the processes associated with three dimensional fiber network, 
it is necessary to develop a mathematical model of the fiber network. The fiber networks 
under consideration have no value to the realm of engineering unless they can be 
manufactured. In this section the methods of forming 3-D braids and multiaxial warp knits 
(MWK) using existing mechanisms are investigated. The models developed herein relate 
the nature of the fiber networks to the geometric properties of the resulting fabric through 
an understanding of the fabrication machine. 

The interest of this section is to develop processing models of 3-D braid and MWK which 
predict the geometric parameters of the braided and knitted fabric and relate them to the 
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processing variables. In general, the specific geometry of a fabric to be formed is affected 
by the size of the yams and their placement within the machine, and the formation of 
geometrical shapes can be achieved by the proper positioning of longitudinal, transverse 
and through thickness laid-in yam systems. In addition to the geometrical effects of these 
systems, the packing of yams within the structure will also be affected. Furthermore the 
mechanical properties of the resulting structures are affected. 

2.2 3-D braiding 

2.2.1 3-D Braiding Process 

3-D braiding technology is an extension of well established two-dimensional (2-D) braiding 
technology in which fabric is constructed by the intertwining or orthogonal interlacing of 
two sets of yams (braiders and axials) in order to form an integral structure. 3-D braids, 
widely used to reinforce structural composites, provide enhanced properties and the 
possibility for near-net-shape reinforcement, as opposed to conventional weaves which are 
layered to form composites. 

A generalized schematic of a 3-D braiding process is shown in Figure 2.1. Axial yams, if 
present in a particular braid, are fed direcdy into the structure from packages located below 
the track plate. Braiding yams are fed from bobbins mounted on carriers that move on the 
mack plate. The pattern produced by the motion of the braiders relative to each other and 
the axial yams establish the type of braid being formed, as well as the microstructure. 

3-D braids have been produced on traditional hom-gear machines for ropes and packings in 
solid, circular, or square cross-sections. A number of new machines without the traditional 
horn-gears have been developed to create 3-D braids with complex shapes. Track-and- 
column and 2-step braiding processes are two examples of the new developments. The 
mechanism of these braiding methods differs from the tradiuonal hom gear method only in 
the way the carriers are displaced to create the final braid geometry. Instead of moving in a 
continuous Maypole fashion, as does in the solid braider, these 3-D braiding methods 
invariably move the carriers in a sequential, discrete manner. 
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Figure 2.1 Schematic drawing of a generalized 3-D braider. 

Figure 2.2 shows a basic loom setups and the carrier motions during four step braiding in 
a rectangular configuration. The carriers are arranged in tracks and columns to form the 
required shape and additional carriers are added to the outside of the array in alternating 
locations. Four steps of motion are imposed to the tracks and columns during a complete 
braiding machine cycle, resulting in the alternate X and Y displacement of yam carriers. 
Since the track and column both move one carrier displacement in each step, the braiding 
pattern is referred to as lxl. 0° axial reinforcements can also be added to the track-and- 
column braid as desired. The formation of shapes, such as T-beam and I-beam, is 
accomplished by the proper positioning of the carriers and the joining of various 
rectangular groups through selected carrier movements. 
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Step 3 Step 4 


Figure 2.2 Formation of a rectangular 3-D track and column braid 

2.2.2 Quantification of the Braiding Operation 

The understanding and quantification of the braiding process is necessary for later 
idendficauon of the unit ceil geometry. One way of looking at braiding is to consider the 
process as a permutation operator, in which the operator maps from a two-dimensional 
lattice to another two-dimensional lattice and the braided fabric is considered as the trace of 
the map. This approach is taken by Ko and Pastore [2.2], Although it provides a 
mathematical description for braiding motions, considerable simplification has to be made 
to apply the general theory to practice. 

In the present work a direct approach is taken where the motion of yam carriers are traced 
throughout a braiding cycle. This trace along with other relevant processing information is 
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used to generate the path of the yam and the shape of the braided fabric. This approach can 
easily be translated to computer codes to simulate braiding [Chapter 6]. 

2.2.3 Notations and Conventions used in the Model 
Geometry 

The loom is consists of m tracks and n colums. The outermost tracks and columns do not 
move. A loom matrix By is defined such that By=l if a carrier position contains a bobbin. 
Otherwise By=0. The lower right comer bobbin is taken as the reference bobbin and is 
defined to be at position 1,1. 

Braiding cycle 

The four step braiding cycle, as shown in Figure 2.2 can be broken down into the 
following operations: 

Step 1: The odd tracks move left and even tracks move right, (left -1; right +1) 

Step 2: The odd columns move down and even columns move up. (up +1; down -1) 

Step 3: The odd tracks move right and even tracks move left. 

Step 4: The odd columns move up and even columns move down. 

The model is developed for lxl braiding pattern, i.e, all tracks and columns only move 
zero- or one-carrier distance in any steps, but it can easily be extended to 1x2. 1x3, 2x3 or 
even complicated mixed patterns. 

Relational Operators 

The relational operators used in the formulation may not be apparent to all the readers and 
hence are explained here 

Relation 

This is a conditional operator and is used to express an “if-then” relationship. p->q implies 
that if p is true q takes place. 

i.e. (b>10)-»a=2, implies that if b is less than 10 then a is equal to 2. 
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Conjunction “ A " 

The conjunction operator is equivalent to logical AND used in most high level 
programming languages. 

i.e. (b>10) A (c>9)— »a=2, implies if b is less than 10 AND c is less than 9 then a is equal to 
2 

The Position Tensor. P£ 

The positon tensor traces the bobbin as the braiding progresses. 

The subscripts i and j denote the bobbin that is being traced. (i,j) is the position of the 
bobbin at the start of braiding. 


The third subscript k can have two values. T or C. T is for track and C is for column. 


The superscripts p and s stand for period and steps respectively. In four step braiding each 
period consists of four steps. Writing out the third subscript: 


pp.» — 
r :jk - 


l P "c. 


( 2 - 1 ) 


Here. is the track position of the bobbin i.j at the end of p period and s steps. 
Similarly, P^is the column position. If the position tensor is known as the braiding 
progresses the trace of the yam can be easily determined. 


2.2.4 Representation of the Braiding Cycle 

Startup (0 cycle, 0 step) 

At the startup phase all the bobbins are at their inital positons. This is represented as 
follows: 



The inital cycle 

Step 1: In step 1, the odd tracks move left (-1) and the even tracks move right (+1). This 
is formally represented as: 
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. |( B , = l) a (1 < < m) -» P“ + (- l) p " 


The condition By=l in the conjunction states that we trace only the looms which have 
bobbins. Since this is always true, it will be omitted but implied. With this simplification: 


p ut 1 J(l<PK<m)->PE +(-!)* 


Step 2 : In step 2, the odd columns move down (—1) and even columns move up (+1). 


(1 < Put < n )— > P«i +(“ 1 )^ 


Step 3 : In step 3, the odd tracks move right and the even tracks move left. 


P S1 J(l<P$?<m)-*P« -(-1)* 


Step 4 : In step 4, the odd columns move up and even columns move down. 


> 0.4 1 

« T L 


P?c) 1(1 < PijT < n)-» Pyc -(- 1 ) ?:,T 


Generalization to subsequent steps 

Once the motion of the loom for the first braiding cycle (which consists of 4 braiding steps) 
is simulated, the subsequent steps can be progressively generated by equeations (2-8) to (2- 


_ (l<PP- c w <m)->Pr T u +(-l) J 


PP-M 

*ijC 


(1 < PSt < n) — > PS’r + (-1) F 
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(2-10) 



i (l<P£<m)-»pM-(-l) 

pP.2 

r ijC 





DP.3 
r ijT 

P ' 3 ^ « \ ^ DP ’ 3 _ /_n p ?r 


(l<PK<n)->Pg-(-l)' 


( 2 - 11 ) 


An example 

We use the formulation to trace the orbit of bobbin initially at (3,3) on a (5x14) loom. Here 
i=j=3, m=5 and n=14. 


The machine setup and the orbit of the bobbin is shown in Figure 2.3. 
Initialization 



inital cycle 
Step 1 



( 2 - 12 ) 


(2-13) 


(2-14) 


(2-15) 


Step 4 
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( 2 - 16 ) 


(2-17) 


(2-18) 


(2-19) 


( 2 - 20 ) 


We can carry on the process to produce the following trace for the bobbin. The process can 
be carried on to show that the orbit of the bobbin initially at (3,3) is given by: 

(3,3), (1,4), (2,4), (4,2), (5,2), (7,4), (8.4), (10,2), (11,2), (13.4), (13.5), (12.3), 
(10,1), (9,3), (7,5). (6.3), (4,1) (the sequence repeats itself) 

In each ordered pair the first element is the track position at the end of the cycle and the 
second element is the column position. 

The motion of the bobbin is periodic. While to trace the yam it is not necessary to 
determine the period, the information is usetul in tracing the bobbin on the loom and to 
study the braid cross-section. The knowledge also helps in conserving computer memory 
when the formulation is used to simulate the machine motion on a computer. The period is 
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measured by comparing the position of the bobbins at the end of each cycle to the initial 
position. A full period is braided when the following equation becomes true. 

V(i.j.k).P£-P£ (2. 


Equation (2-21) implies that all the bobbins move to their initial positions. 

Initial Position of Bobbin (3,3) 



= Yarn Carrier 


= Yarn Carrier in Orbit of (3, 


■ 

9 

m 

■ 

as 

■ 

w 

SI 

m 

u 


■ 

m 

m 

m 


m 

m 

m 

m 

m 

h 

si 

si 

!H 

m 

m 

9 

■ 

m 

m 

n 

as 

m 

Si 

m 

a 

m 

m 

m 

m 

m\ 

m 

m 

m 

S3 

m 


m 

m 

m 

m 

m 

m 

m 

u 

■ 

m 

■ 

m 

ii 

m 

u 

% 

u 

m 

a 


u 

■1 


Fig 2.3 Machine Setup of 5x14 loom for lxl braiding and the orbit of bobbin initially at 

(3.3) 

We note that the orbit of the element (3,3) do not pass through all possible bobbin 
positions. Out of 51 possible bobbin posiuons the orbit passes through 17 of them. If a 
single yam passes through all possible bobbin posiuons, the braid pattern is called fully 
integrated. In a fully integrated braid, the orbit of a single yam determines the entire 
pattern. The braid shown in Figure 2.3 is not fully integrated. To obtain the complete 
braiding pattern we need to examine the orbit of other elements. Let us examine the orbit of 
outer elements (3,2) and (3,4). The orbit of (3.2) is given by: 

(3,2), (5,4), (6,4), (8,2), (9,2), (11,4), (12.4), (14,3), (12,1), . (11,3), (9,5), (8,3), 
(6,1), (5,3), (3,5), (2,3), (2,2), (3,2) 
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The orbit of (3,4) is given by: 


(3,4), (4,4), (6,2), (7,2), (9,4), (10,4), (12,2), (13,2), (13,3), (11,5), (10,3), (8,1), 
(7,3), (5,5), (4,3), (2,1), (1,2), (3,4) 


We see that the three orbits are mutually exclusive and they together pass through all 5 1 
possible bobbin position. Therefore for our example the orbit of these three elements 
together define the complete brading pattern. 

As an example of an integrated braid we examine a 5x13 braid. Here ra=5, n=13. Again 
using equations (2-2) to (2-11), we determine the orbit of element (3,3) as follows: 


(3,3), (1,4), (2,4), (4,2), (5,2), (7,4), (8,4), (10,2), (11,2), (13,3), (11,5), (10,3), 
(8,1), (7,3), (5,5), (4,3), (2,1), (1,2), (3,4), (4,4), (6,2), (7,2), (9,4), (10,4), (12,2), 

(13.2) , (11,3), (9,5), (8,3), (6,1), (5,3), (3,5), (2,3), (1,2), (3,2), (5,4), (6,4), (8,2), 

(9.2) , (11,4), (12,4), (12,3), (10,1), (9,3), (7,5), (6,3), (4,1) 


As it can be seen from Figure 2.4, the orbit passes through all possible 47 bobbin 
positions. This also shows that slight change in processing parameters may cause a drastic 
change in braid pattern. 

2.2.5 The Yarn Trace 


Once the elements of the position tensor P is known, the yam trace can be calculated. In 
calculating the yam trace, it is assumed that a motion L of the loom will cause the yam at 
the braiding plane to move a distance Lr. r is defined as the braiding ratio. The takeup is 
assumed to take place at every second step. 

With the above assumptions for yam starting from bobbin ij, we have: 

Xr = P^.r (2-22) 

Here. Xy* is the x-coordinate of the yam. The origin of the coordinate system is at the 
braiding plane corresponding to yam 1.1. 
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Initial Position of Bobbin (3,3) 


□□□□□□□□□□□□I 

■□□□□□□□□□□□I 

IGllHEIOli 


= Yam Carrier ^ = Yarn Carrier in Orbit of (3,3) 



Fig 2.4 Machine Setup of 5x13 loom for lxl braiding and the orbit of 

bobbin initially at (3,3) 


Similary, 


Y C1 = P CJ r 
I j r ijc- r 


and. 


Z CJ = 


_ r 1 S | , 

2p + Int< — > 


pick 


(2-23) 


(2-24) 


2.3 Multiaxial Warp Knit. 

Knitted fabrics are interiooped structures wherein the kintting loops are preodued by the 
introduction of the knitting yam either in the cross machine direction (weft knit) or along 
the machine direction ( warp knit). The unique feature of the weft knit structures is their 
conformability. But. the most undesirable feature, from the structural reinforcement point 
of view, of the weft knit structure is their bulkiness which leads to the lowest packing 
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density, or lowest level of maximum fiber volume fraction compared to the other fabric 
preforms. Thus, the multiaxiai warp knit (MWK) 3-D structures are more promising and 
they have undergone a great deal more development in recent years. 

The MWK fabric systems consist of warp (0°), weft (90°), and bias (±0) yams held 
together by a chain or tricot stitch through the thickness of the fabric, as illustrated in 
Figure 2.5. The major distinctions of these fabrics are the linearity of the bias yams; the 
number of axes; and the precision of the stitching process. The way of introducing the bias 
yams is to lay in a system of linear yams at an angle. Depending on the number of 
gui debars available and the yam insertion mechanism, the warp knit fabric can consist of 
predominately uniaxial, biaxial, triaxial or quadraxial yams. The latest development, as 
shown in Figure 2.6, in the impaled MWK is the LIB A system wherein six layers of linear 
yams can be assembled in various stacking sequences and stitched together by knitting 
needles piercing through the yam layers. 

Theoretically, the MWK can be made to as many layers of multiaxiai yams as needed, but 
the current commercially available machines only allow four layers (the Mayer system) of 
0°, 90°, +0. and -9 insertion yams, or six layers (the LIBA system) of 2(90°), 0°, 2(+e), 
and -9 insertion yams to be stitched together. .All layers of insertion yams are placed in 
perfect order each on top of the other in the knitting process. Each layer shows the 
uniformity of the uncrimped parallel yams. To ensure the structural integrity, it is clear that 
the 0° yams cannot be placed in either top or bottom layer. The insertion yams usually 
possess a much higher linear density than the stitch yams, and are therefore the major load 
bearing components of the fabric. 
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Figure 2.6 MWK LIBA system. 
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The MWK fabric preforms having four directional reinforcements similar to quasi-isotropic 
lay-up can be produced in a single step. Besides good handleability and production 
economics, the MWK fabric preforms also provide the conformability to complex shapes, 
the flexibility in the principal yam directions, and the improved through-the-thickness 
strength. The MWK fabrics have now been used for playing field surfaces, rubber coated 
fabrics for lifeboats, rescue floats, inflatable sport boats, heavy-duty tarpaulins, geotextiles 
and filter fabrics, and automotive air bags. 

The key geometric parameters of the MWK fabric preforms, which affect the reinforcement 
capability and the composite processability, include the number of yam axis, the orientation 
of bias yams, total fiber volume fraction, pore size and pore distribution, and percentage 
of stitch fibers to total fiber volume. The process variables adjustable to control the MWK 
micro-structure include the type of knit stitch, the ratio of stitch-to-insertion yam linear 
density, the orientation angle of bias yams, and the thread count The concept of a unit-cell 
is used to establish the relationship between the geometric parameters and process 
variables. 


2.3.1 Knitting Process 

The principal mechanical elements used in knitting are needles, which form the loops to 
interlace the linear insertion yams. Therefore, knitting process can be fully understood if 
the motion of the needles is described. The loops in knitted fabrics are formed essentially 
on a very similar principle. Following Thomas [2.3], the looping process is demonstrated 
for a single latch needle by the consecutive steps shown in Figure 2.7. The general knitting 
action of a latch-needle machine can be found in [2.4], 

Consider the needle which has at its stem a loop already framed during the course of the 
knitting process, as in Figure2.7(a). A thread is then placed under the hook of the needle. 
The loop is restrained in its position whereas the needle is allowed to move through it. As 
the needle moves downward, the existing loop will push the latch and close the hook 
(Figure 2.7c). When the top of the hook reaches the level of the existing loop (Figure 
2.7d), this loop is pulled out of the way by the yam tension. Then as the needle moves 
upward again, the thread in the hook opens up the latch, and it becomes the next 'existing' 
loop. More loops are generated as the process repeats. 
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Figure 2.7 The needle cycle in knit fabrics. 


Generally, the warp knitted loop structur is made of two parts. The first one is the loop 
itself, which is framed by the yarn being wrapped around the needle and drawn through the 
previous loop, as described above. This part of the structure is called an overlap. The 
second part is the length of yam connecting the loops, which is called an underlap. It is 
formed by the shrogging movements of the ends across the needles. Since the underlap is 
connected to the root of the loop, it causes, due to warp tension, an inclination to the loop 
structure. 

Two different lap forms are used in warp knining, depending on the way the yams are 
wrapped around the needles to produce an overlap. When the overlap and the next 
underlap are made in the same direction, an open lap is framed, shown in Figure 2.8(a). 

If. however, the overlap and the following underlap are in the opposition to one another, a 
closed lap is formed, shown in Figure 2.8(b). 
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Figure 2.8 Open and closed lap configurations 


For a chain (pillar ) stitch, it is formed when a needle is being lapped continously by the 
same guide. Since the guide bar does not lap the adjacent needles, there are no sideways 
connections. The chain lapping movement can be open (see Figure 2.9), closed (see 
Figure 2.10) or can be a combination of closed and open laps. The more common open lap 
chain construcuon. shown in Figure 2.1 1. is fronted when the guide laps the needle 
alternately from the right and the left The chain notations, as derived from Figure 2.9, are 
0-1 for the first course and 1-0 for the next. To produce a closed lap chain, the guide has 
to lap the needle continously in the same direction and the chain notations are 0-1 for all 
courses, as illustrated in Figure 2.10. 
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Figure 2.10 A closed-lap chain stitch. 
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Figure 2.1 1 An open-lap chain stitch construction 

For the multi-axial warp knits(MWK), straight warps, straight filling and off-axis lay-ins 
are introduced and the stitch yam is looped over the intersection of filling yam and off-axis 
yams. As shown in Figure 2.5, the unit cell of the MWK includes warp (0°), weft (90°), 
and bias (±9) yams held together by a chain or tricot stitch through the thickness of the 
fabric. The size of the unit cell depends on the orientation of off-axis yams and the 
diameter of the insertion yams. In this report, the chained MWK is mainly focused, and 
the tricot MWK can be studied in a similar way. 


2.3.2 Stitch Loop Model 

As mentioned previously in this section, the stitch loop is formed by needle motion. But, it 
is very difficult to trace the stitch yam path in the 3-D space when it is guided by the 
needle. Unlike the machine motion of 3-D braid, the relationship between machine motion 
and stitch yam path in MWK is imposibly established. Therefore, assumptions have to be 
made in order to describe the stitch yam path. The following points should be considered 
before a loop model is to be built. 
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1. The shape of the loop, loop inclination and fabric weight may change slightly without 
changing the course and wale counts. 

2. The configuration of the yarn within the knitted cell is affected by a great number of 
variables such as yam properties, warp tensions, ta-up tensions, yam lubrication , etc. 

It is proposed that the loop shape of the fabric in the machine state is more likely to be 
determined by the physical pull of the take-up mechanism and the accommodation of the 
root configuration, than by the bending forces. Taking these facts into account, the 
following loop model for the fabric in the machine state, shown in Figure 2.12, is 
suggested. 



Figure 2.12 The 'machine state'' loop modei 


The model can be divided into three parts: 

1. The loop shape which is illustrated in Fig.2.13 can be described as 
the loop's head and the straight arms. The loop head is on 
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horizontal plane, z = 0, while the loop arm is in a 3-D space. 


Rloop = 3r 



r : radius of stitch yam 
Rloop : radius of loop head 


Figure 2.13 The loop pan of the machine state ' loop model. 


loop's head = 2 TCRloop = oTCr 


loops arm = 4 \'(H' - r>- - 4r 2 




(2-25) 






2. The loop-over, which is the loop over insertion yams, can be described in the 
following figure 2.14: 



Rins : radius of t he insertion yam 
r : radius of t he stitch yam 

R’ : radius of the circular arc, which partially wrap around the off-axis insertion yam 

R’ = Rins/cosidS 0 ) = v 2 Rins 

Figure 2.14 The loop over the insertion yams from the loop model 


loop-over $ head = 2 TCR’ = 2 \ 2 rrRins 

H * 

loop-over s armt entering; = Z'\l “Rins + r) 2 - ( - 4r - R') 2 (2-26) 

loop-over's armt leaving) = - A,/ (TRins + r) 2 + ( - 4r - R’) 2 


23 









The second kind of root for the yam leaving the loop-over can be calculated as: 


2jct ( — ) 

K 


where a = tan- If 


7Rins 


- 4r - 42 Rins 


■) 


The amount of yam in the root of the loop is therefore: 


a 9 

loop’s root = 2m ( — ) + 2JCr(l- — 

K 7l 


Hence, the complete model can be described as: 


a 


6 


L = 6rcr - 4 \‘(H' - t)~ - 4r“ - 2V2 TlRins + 2jrr ( — i + 2TCr (1 - — ) 

X K 


+ 2 'y TRins + n 2 - ■ - 4r - R’) 2 (2-28) 

/ tT"* 

+ 2 \j CTRins + r» : - • - 4r - R’) 2 

By applying the loop model stated above, he entire shape of the stitch yam can be 
described provided that processing parameters are given. Thus, the key points of the unit 
ceil of a MWK preform, including stitch yams and insertion yams, can be obtained through 
a geometric analysis. The details of how to acquire the key points mathematically will be 
presented in Chapter three. 
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Chapter -3 Geometric Modeling 


In this chapter, the theoretical framework of computer aided geometric modeling for the 
structural analysis of textile reinforced composite materials is presented. Based on the 
processing model developed in Chapter 2, the illustration of generating key points for each 
yam path is presented for 3-D braid and Multiaxial Warp Knit composites, respectively. 

3.1 Theoretical Background 

The first step in the construction of a geometric model of a textile preform is to model the 
pretorm processing technique. The machinery which forms the textile preform dictates the 
geometry of the yam path. By proper modeling of the process, it is possible to generate the 
basic information associated with the yam pain. The modeling of the process has been 
discussed in detail in Chapter 2. 

The aigonthm is deveioped such that the designer can enter concise information about the 
machine parameters and have the aigonthm output a set of "knots ' which descnbe the yams 
and yam paths within the tabric. The knots are a sequence of points defining the center of 
the yam in such a way that the enure yam can be regenerated from this finite set. 

3.1.1 Yarn Path Modeling 

Given a set of knots irom the modeling ot me textile structure, it is now necessary to 
predict the geometnc behavior of the yam Detween these points. Tne yam path wiil be 
dictated by the matenai properues. the cross-secuonal geometrv or the vam, and the motion 
of the yam relative to the other yams. 

The material properties of the yam are moaeied in terms of the minimization of strain 
energy on the yam element. Strain energy is minimized through the incorporauon of a B- 
spiine funcuon to generate the path of the yam. The basic assumpuon of the B-spiine is to 
minimize the term [3.1] 


L 
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where. 


and. 


L= total length of the yam 



r = g(s) = spline function of the yam path 
s = path length. 


since r" is closely related to k (the curvature ), the minimization of this terra minimizes 
curvature, and hence the bending strain (£ = E k) of the yam. Additionally, the use of the 
Bspline minimizes the twist in the yam [3.1] In this work a constant cross-section has 
been employed, which twists according to the path of the yam as defined by the B-spline. 


Given a set of points describing the center of a yam i, jp j , it is necessary to construct 

additional control points to form the B-spune. These control points are formed in such a 
way as to provide differential continuity for me yam path. To construct the curve around 
point P|, four control points are needed, as illustrated in Figure 3.1: 



Figure 3.1. Construction of Control Points for B-Spline Generation . 


To assist in the construction ot the points, a fifth point is introduced., 
which serves only to make the calculations more readable. These points 
can be constructed in groups. The first two control points are the mid- 
points of the vectors connecting the point to its predecessor and 
successor: 


i Pi +Pj-l 
M) “ 2 


(3-2) 
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r i p i+l jfhl 
= 2 


(3-3) 


The mid-point of the vector connecting these two points is not a control point, but is useful 
for the construction, and is defined as: 


d+d 

r i _ 3 0 

v^4 — n 


(3-4) 


The remaining two control points are formed so as to be parallel to the vector formed 
between C^, and pj. They can be given as: 

C\ = cjj + Pi - C\ (3-5) 

d = d + Pj - C\ (3-6) 

This construction guarantees that the point P; is coilinear with C l , and Ci,, and in fact is the 

midpoint of the corresponding vector. The vaiue of this definition is that it provides 
interpolation of the point Pj in the construction of the B-spiine [3.2]. 

Havins constructed the control points, it is now possible to generate the path of the center 
of the yam using the B-spline. The path of the yam is given as: 

a 

£B};kVt)wlCj 


yB i;k {t)wl 

i=G 

where rtt) = yam spauai earn as a funcuon of arc length 

wl = weighting function 
and Bi.kif) = B-spline of order k 

The B-spline is defined as: 

if, 

tj < t < tj+1 

Bj,k(t) = 1 



else, 


B j.k (t) = ^ %,«) * t^V^ t BwwtO <M> 


Since the B-spiine is constructed using four control points, in order to achieve interpolation 
of the knots, the order of the spline is 3 (k=3). 

The B-spiine, r(t) now describes the path of the center point of the yarn as it moves from 
point to point within the preform. With the construction of an analytical function 
describing the path, it is now possible to identify certain critical parameters associated with 
the mechanical response of the yam. The twisting angle on the yam is given by 

(Ks) = cos'Kvq • v^) (3-9) 

where v o s ,p i+l ’ p i) ® ( p i ' ^i-i- 
V 1 = l ' p i+2 - p i+l) ® fp i-i - p i> 

as illustrated in Figure 3.2. 



The twist on the vam between two points can be given as: 

L 

J0(s) ds 

T = ~ f (3-10) 


where. T = twist per unit length, or torsion. 




Si milar ly, the bending on the yam path can be given by 

cos(<j>) = (P i+ i - Pj) • (Pi - Pj.j) (3-11) 

where <j) = bending angle, as shown in Figure 3.3. 



With these relationships, the basic geometric properties of the yam center-line are known, 
and it is possible to construct the three dimensional form of the structure. 

3.1.2 Solid Modeling of the Yarn 


The solid modeling of the yam is achieved by sweeping a yam cross-section along the 
center-line oi the yam as determined using me algorithm described above. Since the typical 
yam does not have a circular cross-section, it is necessary to identify the orientation of the 
cross-section during the sweeping process. An elliptical cross-section has been chosen to 
represent a typical yam. and has oeen shown to be a reasonable reoresentation. 


The ellipse is actually an n-gon resemoiing an ellipse to an arbitrary cesree. The points on 
the cross section, e;. are calculated as 


a cos 



n i 
. • j»t I 
bsin— ; 
n J 


(3-12) 


In this way, a surface representation ot the physical yam can be constructed. The 
orientation ot the ellipse is normal to the yam path, and the major axis of the ellipse is 
twisted in accordance with the yam path twist. .An initial orientation of the ellipse is 
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defined for the yam at z=0. As twist is identified in the yam path, the major axis is rotated 
at the same amount 

The yam now consists of a surface representation. This surface can be rendered in color or 
shading for informational display. Any of the suitable shading tec hniques work well with 
this structure since it now consists of connected polygons. It is only necessary to identify 
the viewpoint the light source, and the reflectivity of the surface. The result is a color 
shaded rendering of the individual yams within the fabric. 

£2 Application to the Braiding Model 

The processing model developed in Chapter 2 can be used to develop the individual 
motions of the tracks and columns in the braiding machine. Based on the position tensor 
operations, the positions of a bobbin can be identified for each cycle as the braiding process 
progresses. The positions stand for the projection into the x-y plane of the path of one yam 
in a typical braiding machine. 


The yamo(i,j) stands for the initial position of the yam at the i 1 * 1 track and column in the 
loom. The initial position ot the yamti.j) is (Xq(i,j), Yo(i,j), Zo(i,j)), where Zo(i.j) = 0. 
From Equations (2-3) to (2- 1 1), the posiuon of the yam<i,j) is defined by Xk(i,j) and 
Yt(i.j) after k cycles of track/column motions. Thus the sequence of points 
{ (Xo(i,j),Yo(i,j), Zo(i,j)),(Xi(i,j),Y}(i.j), Zi(i.j)), ... ,(X[c(i,j), Yk(i,j), Zk(i,j)) }, where 
Zic(i.j) — k • Az, descnoes the path ot the yamti.j). Az is the distance ot each pick, which 
is determined by the braiding angie and the rauo of track/column movement. The 
sequence, { — iXk.Yk.Z^) } , represents the spatial position of kev points in a yam path, 
and these points can be used as the knots for the three dimensional braid. 

These key points are the basis of generating the yam path using the 3-spiines as discussed 
earlier. The geometric model is complete when we have all the key point information on the 
yams in the braid. The B-spline algorithm is incorporated in a computer code which is 
developed using C language and Sun View graphics on a SUN 3/160. The modeline 
program and its user interfaces is described in Chapter 6. 

The shape of the yam cross-section in a 3-D braid can be determined based on the 
geometric parameters, processing variables and following assumptions: (1). circular yams 
with radius, r, (2). no interactions between yams. Figure 3.4 shows a single yam in space 
with surface angle 0 S after 1/1 track/column motion. 
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9 S = surface braiding angle 
y = braiding angle 


Figure 3.4 The orientation of a yam in space. 


Ay is the distance between two adjacent key points in the y-direction and its value is equal 
to the braiding ratio, r. as mentioned in Chapter 2. The distance of each pick. Az, is 
calculated by: 


Az = Av*cot(0s) = r*cot(0s) = r*tan(\j/)«cos(P) 


the angle 9 can be calculated bv the following reiauon: 


n 1 r COS( 01 , 

9 = tan* 1 [ ; 

tan(0 s ) 


The inclined angle u. or. braiding angle can expressed as: 


” Q . r cost Q) ? 

V = - - 9 = - - tan* 1 [ ] 

tarn 0 S ) 


For circular yams, the yam cross-secdon becomes elliptical on a cut plane parallel to x-y 
plane as illustrated in Figure 3.5: 
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braiding 

yam 


r : radius of a yam 

\\r : braiding angle 

a : major axis of the ellipse 



Figure 3.5 Cross-section of a yam on different planes. 


3.3 Application to Multi axial Warp Knit (MWK) 

Theoretically, a MWK can be made to as many layers of multiaxial yams as needed, but the 
current commercially available machines oniy allow four to six layers of 0°, 90°. +0. and -0 
insertion yams to be stitched together. Based on the processing model presented in Section 
2.3, the MWK structure consists of four basic components: warpiO 1 ! yams, weft (90°) yams, 
bias (±9) yams, and stitch yams through the thickness of the fabric. The dimensions of a unit 
cell of MWK composites can be expresses as H x 2H’ x Hz, where H is the spacing between 
two o° insertion yams. H' is equal to H(cot0), and Hz is the thickness of the unit ceil. The 
geometric relationships of the unit cell is illustrated in Figure 3.6. 
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Figure 3.6(aj The planar dimension of a unit ceil of MWK composites 

As can be seen in Figure 3.6(a), the 0° insertion yams are H apart, while 90° insertion 
yams are H apart. The off-axis insertion yams travel through the diaaonal orientation of 
the H x H rectangie in the through-thicKness projection. As in the thickness direction 
shown in Figure 3.6(b), the z-coordinate ot the center point ot each insertion vam are iisted 
hi Table 3-1. Thererore. the key points tor representing insertion yams can be determined. 
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Hz 1 Y 



Hz = 4<2Rins) + 2(2Rsy; 

Rins : Radius of insertion yam 
Rsy : Radius of stitch yam 


Figure 3.6(b) The through-thickness dimension of a unit cell of MWK composites 


Table 3-1 Z-coordinate of the center point of each insertion yam 


insertion vam 

z-coordinate 

VO 

o 

o 

Rins + Rsv 

0° ! Rins + 3Rsv 

i 

' Rins + 5Rsv 

-0 

Rins + 7Rsv 
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For the stitch yarn, based on the processing model given in Section 13, twenty-seven key 
points are identified in order to represent the stitch yam within a unit celL The location of 
each key point is marked and shown in Figure 3.7. 




Figure a. , Locauon ci each key point in a unit ceil. 


These key points are usually locates at the intersections of two kinos of curves, for 
instance, the intersection of straight line anc arc < loop). The coordinates of each key point 
can be expressed in terms of the processing parameters, such as H. H\ Rinsert. Rloop, 
and Rsv. The origin of the coordinate system is located at the lower-left comer of the unit 
ceiL For example, at point ?fl. it can dc reaciiv calculated that 

u 

x -coordinate : — - Rsv 

OWGSNAL PAGE IS 
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y-coordinate : - Rsy - Rloop 


For another example, at point #2, the geometric relationship is shown in Figure 3.8. 



x' = 2Rsy ( 


4Rsv 

IT ‘ 


H H 

Thus, the x-coordinate is equal to — - x' - Rsy [ = — - Rsy - 4Rsy( ) ] 


The distance z can be calculated in the following trigonometric operation: 

_> _ , w H’ - 4Rsv , _ , Rsv 

z - (- 2 RsyX H , . 3Rsy ) - - -Rsv + H > _ 3Rsy 


' / 





Thus, the coordinates of the point #2 


x - coordinate : — - Rsy - 4Rsy( ) 

y • coordinate : 0 

z - coordinate : • 2Rsv + u , Rs 7 n — 

H - 3Rsy 

For another example. Figure 3.9 shows the geometric relationship between the key point 
#8 and point #9. 



Figure 3.9 Geometric reiancnship for point #8 and #9. 
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Thus, the x -coordinate of key point #8 is 


x = j + x’ = j + 3Rsy(^) 

For the y-coordinate. the enlarged geometric relationship is shown in Figure 3. 10. The 
horizontal distance, d, can be determined from triangular relationship: 

d = 3Rsytan(0) 



Figure 3.10 Enlargea view of geometry around key point #8. 


where 


tan 10 ’) = 


2Rsv 

H’ 


From the geometric reiationship the distance y ' can be calculated bv 


y’ = a sin(0) cost©) = d ( 


2Rsv 


+ dRsyZ)!^ ( (H*- + 4Rsy^)l/2 ; 


H’ 


or. 
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, „„ , 2Rsy w 2Rsy w H’ 

y = 3Rsy( - H^X ( h* 2 + 4Rsy2)l/2> < (jfl 7 4R Sy 2)l/2 ) 

12Rsy3 


H’2 + 4Rsy 2 


12Rsv 3 

Thus, the y-coordinaie of the point #8 is -Rsy - ■ ; 7 . The space coordinates of 


the key point are 


H’ 2 + 4Rsy 2 ' 
(x,y,z) = ( " + 3Rsy(^). -Rsy- H .' 2 f RS 4 ^ sy2 


, -2Rsy) 


From Figure 3.9, the coordinates of the key point #9 can readily be determined as: 


H 


(x,y,z) = f — + Rsy, - Rsy - Rloop. 0) 


For the key point #10. where the stitch yam starts to “wrap” around the off-axis insertion 
yam. its geometric relationship with key point? 1 1 and key point# 12 is described in Figure 
3.11. 

In the case of the orientation of insertion yam :s ±45°, the distance A. on the cut cross- 
Rins . 

section, is equal to — — .or, v 2 Runs, rrom the observation, the coordinates of the 

cos(45 ) 

key point# 10 can be expressed as: 

H H’ 

(x.y.z) = \ 2 Rins. 7Rins i 

Similarly, the coordinates of the key point? 1 i can be found as: 

H H ’ 

'"x.y.z) = i — . — . SRins t Rsv) 
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rigure 3.1 1 Key points in the loop-over portion. 


The location of key roint*13 is a shift of H' of the key point#9 in y-direcdon. The 
coordinates of the key pointffl3 are 

ix.y.z) of the point?# + <0.H\0) = ? ^ Rsv. H’ - Rloop - Rsv. 0 ) 


Similarly, the location of key point# 14 is a shift of H’ of the key point#8 in v-direction. 
The coordinates of the key point?? 14 are 


(x,y,z) of the point#8 + (O.H'.O) 
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For the key points from #15 to #19, their locations are a shift of H’ in y-direction to the key 

points from #7 to #3 (in reverse order), respectively. The key points #20 and #14 are 

H 

symmetrical with respect to the line : x = y . The x -coordinate of key point #20 can be 
found by: 

H H 

y --(x-coordinateofthekeypoint#14 - y ) 


= H - ( x-coordinate of the key point #14) 



- 3Rsy ( 


2Rsy 

H’ ; 


Thus, the coordinates of the key point #20 are expressed by: 


( y - 3Rsy(^), H’-Rsy- 


12Rsv 3 


H’ 2 + 4Rsy 2 


, -2Rsy) 


In the same sense, the key points #21 and #13 are symmetrical with respect to the line : x 

pj 

— . The x-coordinate of key point #21 can be found by: 


H H 

— - ( x-coordinate of the key point if 13 - y ) 

= H - f x-coordinate of the kev point #13) 

H D 

= - - Rsy 


The coordinates of the key point #21 are given by: 
( ^ - Rsy, H’ - Rloop - Rsy. 0 ) 


The locations of the key points #22.#23,#24. #25 and #26 are a shift of H’ in y-direction 
to the key points #10. #1 1. #12. if 21 and #20. respectively. The location of the key point 
#27 is a shift of 2H’ in y-direction to the key points #2. and the coordinates of the key 
point #27 are 


(x. y, z) = (x, y, z) of the key point #2 + (0, 2H\ 0) 
= (S-Rsy-4Rsy(^). 2H\ - 2Rsy + 
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The results of the geometric analysis for the stitch yam will be used as input for solid 
modeling of the MWK. The summary of the coordinates of each key point is shown i 
Table 3-2. 





x - coordinate 


y-Rsy 


S-Rsy-4Rsy<?52L) 


y - Rloop 


f - f V 2 Rsy 




y + ^V2 Rsy 


y - Rloop 


H / 2Rsy ~ 

2 "t" 3 Rsy ( ) 


H „ 
y + Rsy 





y + Rsy 


y + 3Rsy(^) 


y - Rloop 


H 3 

y + JV2 Rsy 


y - coordinate 


Rsy - Rloop 


z - coordinate 




H' - Rsy - Rloop 


H' -Rsy -Rloop(l- -y- ) 


H' - Rsy 


H' -Rsy -Rloop(l- -y-) 


H' - Rsy - Rloop 


-Rsy 12Rsy3 

y H’2 + 4Rsv 2 


- Rsy - Rloop 


H’ 

y - V2 Rins 



y + V2 Rins 


H' - Rsy - Rloop 


H’-Rsy- 


12Rsv 3 



A JLXvJ T ~ 

H’ 2 + 4Rsy 2 


2H' - Rsy - Rloop 


2H' -Rsy -Rloop(l- y— ) 


2H’ - Rsy 


2H' -Rsy -Rloop(l-4y-) 



-2Rsy 


7Rins 


8Rins + Rsy 


7Rins 


-2Rsy 
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19 

y- Rloop 

2H' - Rsy - Rloop 

0 

20 

—EM 


-2Rsy 

21 

H D 
2 + Rsy 

H' - Rsy - Rloop 

0 

22 

H 

2 

^ - V2Rins 

7Rins 

23 

H 

2 

3H' 

2 

8Rins + Rsy 

24 

H 

2 

£ y- + V2Rins 

7Rins 

25 

H D 
2 + Rsy 

2H’ - Rsy - Rloop 

0 

26 



-2Rsy 

27 

Y - Rsy - 4Rsy( ) 

2H’ 
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Chapter - 4 Visualization Verification 


Based on the techniques described in Chapter 2 and Chapter 3, the preform geometries of 
3-D braid and MWK can be generated by computer simulation program. In order to verify 
the simulated geometries, the microsturctures of real composites were e xamine d by photo 
microscopy. The composites were cut, mounted with resin, polished using polishing 
powder. Pictures were taken after each polish and can be used to compare with simulated 
graphics. 

4.1 Material System 

4.1.1 3-D Braid 

Both 3-D braided PEEK/graphite and Epoxy/Mlliken multiaxial warp knit composites were 
used in the present study. The 3-D braided pretorms were made by using commingled 
Victrex® 150G PEEK/AS4 6K yam manufactured by BASF Structural Materials, Inc. The 
commingled yam contains approximately 280 PEEK fibers (fiber Hia m^ r _ 27 
micrometers) and 6.000 graphite fibers ( fiber diameter - 8 micrometers). The graphite 
fiber volume fraction in each commingled yam is 62.7%. In order to obtain a composite 
with corss section of 0.254 mm (0.1 inch) by 1.27 ram (0.5 inch), with a braiding angle of 
±20°. a lxl construction was produced by using 72 yams in a loom consisting of 18 
columns by 4 tracks. After braiding, the 3-D braided preforms were placed in a carbon 
steel mold and consolidated by using a hot press, and held at 400°C and 1.38 MPa (200psi) 
for an hour. After that, the mold was taken out from the hot press and placed in a cold 
press and held at 1.38 MPa for another 30 minutes. 

4.1.2 Multiaxial Warp Knit 

The Epoxy/Milliken multiaxial warp knit composite used in this study was provided by 
NASA. The fiber architecture design of this material is shown in Figure 4. 1. It was 
fabricated with Hercules AS-4 epoxy-sized carbon tows. The 0° and 90° plies were 
produced with 12 K tows, whereas 9K tows were used for the ±45° plies. This tow 
distribution was required to achieve similar liber area for each ply because of the different 
tow count for the off-axis plies. 
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Figure 4. 1 Fabric design of the muluaxial warp knit composites. 
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Only the 4-ply subgroup (-45°, +45°, 0°, 90°) shown in Figure 4.2 was included for this 
investigation. An additional 4-ply subgroup (+45°, -45°, 0°, 90°) would be required to 
produce a symmetric quasi-isotropic fabric preform. Four 4-ply subgroups were stacked to 
form the full 16-ply laminate. Note that the 16-ply stack, shown in Figure 4.2, is 
unsymmetric. The 4-ply subgroup was produced by chain-knitting the carbon tows 
together with a polyester yam. Some of the 16-ply stacks were stitch together with kevlar- 
29 1500 denier stitching yam using a chain stitch. 


MILL! KEN FABRIC 

AS-4 Carbon tow 


Chain 

stitch 


Knit cattern 


-45 y 
-*-15 

0 

90 


,1 ’ 


fl 




j 

4-piy knit 




ill IH 


3-3 


3 © 3 3 

3 © -3 3 3 90 
0 

nn -45 

iLJ. i +45 

© «© 90 


*w« mini -45 

•/J .O + 45 ' 


(16 layer 
unsymrnetric 
layup) 


Figure 4.2 Stacking seguence of the 16 plied-MWK results in an unsymmetric 

conformation. 
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4.1.3 Sample Preparation 

The unit cell geometries of 3-D braided and raultiaxial warp knit composites were observed 
by the optical microscope technique. In this study, the cross-section of the 3-D braided 
composite was cut perpendicular to the braiding axis, whereas the 0°-direction and 90° 
orientation cross sections of the multiaxial warp knit composites were cut for examination, 
as illustrated in Figure 4.3. 


3-D Braid 


MWK 



cutting planes 


ssssssa 

& 3s 3s 3s 3s 3s 3s 3 
3s 3s 3s 3s 3s 3s 3s 3s 3s 3s 3s 3s 3s3s 3 

3*33333333333333 

333333333333333 

333333333333333 

333333333333333 

333333333333333 

333333333333333 

333333333333333 

333333333333333 

333333333333333 

333333333333333 

333333333333333 

333333333333333 


surface 2 


surface 3 \ 

I 


surface 1 


Figure 4.3 Schematic of sectioning 

The cut-offs were mounted in a holder with Lecoset® 7000 cold-curing resin manufactured 
by Leco corporation. After the resin is cured, the samples were polished with 320, 400, 
600 grid silicon carbide polishing paper , as well as 15, 5 micron aluminum oxide powder 
and 1 micron alpha alumina powder, to obtain smooth surfaces for photos. All polished 
samples were examined under an optical microscope with 32x magnification factor. The 
montages are expected to show the variety of the unit cell structure of these samples from 


48 




different cutting planes. Therefore, in order to study the variation of the unit cell structures 
in 3-D space, a series of montage of the unit cell structures of these composites were 
generated by progressive sectioning of the above samples. By doing so, the 3-D fiber 
conformation of the composite can be recorded on a series of planes with certain distance 
apart and thus, the 3-D fiber architecture can be envisioned and described. 

4.2 3-D Braid 

4.2.1 Computer Generated Graphics 

The analysis of textile composites depends directly on fiber architecture of the composites, 
that, in turn, can be accurately characterized by a computer aided geometric models. The 
geometric models for the 3-D braided preforms are given in Chapter Three, in which the 
geometric model considers the relative mouons of the tracks and columns in the braiding 
machine and generates a mathemadcal simuladon of the braiding process. 

In order to compare the computer-generated graphics with experimental observations, a 
virtual loom of 4 tracks and 18 columns is set up on the computer. Six track/column 
movements are simulated, which will constrcut 1.5 unit cells in the braiding axis. The yam 
cross-section is assumed to be ellipse with the ratio of major axis over minor axis to be 0.5. 
The graphics are generated from a Sun workstation. The details of the implementation of 
the graphics model can be found in Chapter 6. 

Figure 4.4(a) - (h) show a series of cut cross-sections of a 3-D braid. The graphs cover a 
track/column movement cycle. The loom state can be seen from each cut cross-section. 
First, the columns move up/down, as shown from Figure 4.4(a) to (c). In the next step, 
the tracks travel back/forth horizontally, as shown from Figure 4.4(d) to (f). Figure 4.4(g) 
and Figure 4.4(h) show the loom state after column movement and track movement, 
respectively. Theretore, the braiding pattern on each cut plane represents a loom state. As 
we go on examining the braiding pattern along the braiding axis, the braiding pattern will 
repeat after a unit cell s length. The braiding pattern generated from computer is slightly 
different from the experimental observations, especially around the boundary. This is due 
to the compression while the composite was consolidated by hot pressing. 

As observed from Figure 4.4(a) to (h), the braiding pattern changes 90° between Figure 
4.4 (a) and (b). The braiding pattern changes 180° between Figure 4.4(c) and (g). Thus, 
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(C) 
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Figure 4.4 Continued. 
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<h) 

Figure 4.4 Continued. 


53 



if the cut planes are the same as the planes on which the photo pictures are taken, the 
computer-generated graphs should resemble the pictures in the central portion of the 
specimen. 

4.2.2 Experimental Verification 

The unit cell structure of 3-D braided PEEK/graphite composite is a diamond-shape 
structure. The diagonal along the braiding axis was measured at 5.1 mm. The schematic of 
the planes which pictures were taken is shown in Figure 4.5. The distances between each 
planes are 0.2mm, 1.6mm and 2.4mm, respectively, which is long enough to represent the 
variation of the structure within a unit cell. The montages showing the cross-section of the 
composite along the braiding axis in these four different layers are shown in Figures 4.6. 
As it can be seen in the figure, the fiber packing pattern changes along the braiding axis. 
The first two pictures, taken from plane A and plane B, respectively, show no significant 
difference because of the short distance separating the planes. The 90° fiber packing 
pattern change between plane B and plane C suggests that the unit cell consists of at least 
two pairs of yams. These two pairs of yams are othogonal to each other. The two yams 
in each pair are interlaced each other. The distance between plane C and plane D is about 
half length of a unit cell. The fiber pattern on plane C is about 180° to the fiber pattern on 
plane D. By induction, it can be predicted that two montages will look the same if they are 
viewed on the planes 5.1 mm apart. 

In conclusion, the unit cell geometry of 3-D braided composites can be described and 
visualized by the computer software. 
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Figure4.5 Schematic of planes on which pictures are taken for 3-D braid. 
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(a) A -plane 



(b) B - plane 



(c) C - plane 



(d) D - plane 

Figure 4.6 Photomicrograph of 3-D braided composites. 
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4.3 Multiaxial Warp Knit 

4.3.1 Computer Generated Geometry 

The geometric model for the MWK preform is given in Chapter Three, in which the 
geometric model relates the key parameter of fiber architecture to the processing variables 
base on the idealized cross-sections of insertion yams and the path of stitch loops. 
Combined with the optical observations detailed in Section 4.3.2, the geometric model can 
be translated into computer solid models for MWK preforms, which facilitates the 
visualization of detailed internal geometries for MWK preforms and the identification of 
their unit cells. 

In order to obtain an overview of the figures from simulation, a schematic of the planes on 
which the computer generated grapgs were taken is shown in Figure 4.7. In the figure, a 
unit cell, which represents the geometry of the whole structure, is used to illustrate the 
planes. Later, the computer code will generate the corresponding graphs based on the 
planes shown. 

Based on the geometric modeling and computer graphics techniques, the MWK preform 
with 8 unit cells is generated on Sun workstation. Figure 4.8 shows the top view, side 
view, back view and (111) view of the MWK preform, respectively. As can be seen, the 
stitch yam loops over insertion yams with a chain configuration. The geometry of the 
stitch yam is simulated according to the geometric model, described in Chapter two and 
three. The actual geometry of the stitch yam can be determined by the yam tension itself, 
yam friction coefficient and yam bending rigidity. However, this is not in the scope of this 
project. Figure 4.9 shows the (1 0 0) cuting plane and cutting section in top view and side 
view, respectively. From a little distance to the previous cutting plane, a section is cut , as 
shown in Figure 4.10. The cutting planes are corresponding to the polishing planes, 
described in the next section. Since the cutting plane does not cut through 0° insertion 
yam, the 90° yam does not show up in the cut section. 

Figure 4. 1 1 shows the (0 1 0) cuting plane and cutting section in top view and side view, 
respectively. From a little distance to the previous cutting plane, a section is cut , as shown 
in Figure 4.12. Figure 4.13 shows the (0 0 1) cuting plane and cutting section in top view 
and side view, respectively. From a little distance to the previous cutting plane, a section is 
cut , as shown in Figure 4. 14. The cutting planes are chosen according to polishing 
planes. The computer generated graphics will be compared with the photomicrographs 
from the polished sample. 
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Figure 4.7 Schematic of simulated planes. 
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(b) side view 

Figure 4.8 Computer rendering of a MWK preform. 
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(c) back view 



(d) (1 1 1) view 


Figure 4.8 Continued. 
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(a) top view and cut section 



(b) side view 

Figure 4.9 (1 0 0) cut plane of a MWK composite. 
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(a) cut section and top view 



Figure 4.10 (1 0 0) cut plane of a MWK composite, (second cut) 
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(b) front view 

Figure 4.1 1 (0 1 0) cut plane of a xMWK composite. 
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(a) cut section and top view 



(b) front view 

Figure 4.12 (0 1 0) cut plane of a MWK composite, (second cut) 
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(b) top view 

Figure 4. 13 (1 0 0) cut plane of a MWK composite. 
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(b) top view 

Figure 4.14 (1 0 0) cut plane of a MWK composite, (second cut) 
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4.3.2 Experimental Verification 

The schematic of viewing planes for multiaxial warp knit composites are illustrated in 
Figure 4 . 15 . The montages from the first polishing on each plane are shown in Figure 
4.16 - 4 . 18 . In Figure 4 . 16 , 90° • insertion yams and stitch "knots" with regular spacing 
are shown as the composite is observed from a top view. Figure 4.17 and Figure 4. 18 
show the front view and side view of the composite, respectively. The insertion yams look 
like a shape of race track due to the compressive pressure in the thickness direction during 
composite processing. By observations, the unit cell dimension from the montages of the 
first layer can be defined as 2.2, 2.2 and 6.2 mm, in length, width and thickness, 
respectively. The montages after the second polishing are shown in Figures 4. 19 - 4.21. 
Comparing Figure 4. 16 with Figure 4. 19, the 90° layer was polished away; the 45° layer 
shows up. From the front view, only knots or cross-sections can be seen; while the stitch 
loop can be observed from side view, as shown in Figure 4 . 21 . Figure 4.22 shows the 
detailed chain stitch structure ( from surface 3). The schematic of the polishing planes is 
shown in Figure 4.7. 

In comparison to the computer generated graphics. Figure 4.9 and Figure 4.10 are 
similar to Figure 4.16 and Figure 4.19, respectively. Figure 4.11 and Figure 4.12 are 
similar to Figure 4.17 and Figure 4.20, respectively. Figure 4.13 and Figure 4.14 are 
similar to Figure 4.18 and Figure 4.21, respectively. As mentioned earlier, the exact 
shapes of the insertion yams are controlled by the the processing conditions. However, the 
relative locations of the insertion yams can be depicted from the photomicrographs. In 
conclusion, the unit cell geometry of MWK composites can be described and visualized by 
the computer software. 
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Figure 4.15 Schematic of cutting planes for a MWK composite. 


68 




Figure 4.16 Photomicrograph of a MWK composite, (surface 1) 
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Figure 4. 18 Photomicrograph of a MWK composite, (surface 3) 








Figure 4.20 Photomicrograph of a MWK composite, (surface 2, second polishing) 
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Figure 4.21 Photomicrograph of a MWK composite, (surface 3. second polishing) 
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Figure 4.22 Photomicrograph of a MWK composite, showing the chain stitch loop. 
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Chapter - 5. Unit c ell characterization 


5.1 Background 

The traditional approach used in modeling of composites reinforced by three-dimensional 
(3-D) preforms is to assume that their fiber volume fraction and fiber orientation are 
known, either obtained from experimental measurement or provided by preform 
manufactures, the relationship between preform fiber architecture and preform processing 
variables is not considered. In this report, we first examine both 3-D braiding and 
Muldaxial Warp Knitting in the light of fiber architecture, followed by the development of 
geometric models for 3-D braided and MWK structures using a unit cell approach. The 
unit cell geometries of these two 3-D fabrics are identified, and the relationship of structural 
parameters such as yam orientation angle and fiber volume fraction with the key processing 
variables is established. The limiting geometry has been computed by establishing the 
point at which yams jam against each other. Using this limiting geometry factor makes it 
possible to identify the complete range of allowable geometric arrangements for 3-D fabric 
preforms. The identified unit cell geometries can then be translated to mechanical models 
which relate the geometrical properties of fabric preforms to the mechanical responses of 
composite systems. 

5.2 3-D Braiding 

The unit cell geometry ot the track-and-column braid has been investigated by many 
researchers since the early 1980 s [5.1-4], A common assumption made in most of the 
analyses is that the braider yams are oriented along the four diagonals in the unit cell. 
However, the fiber volume fraction of 3-D braid is normally over 0.5, so the yams cannot 
be treated as dimensionless lines to cross each other at the center of the unit cell. This unit 
cell geometry is either oversimplified or incorrect. Li, Hammad and El-Shiekh [5.3] 
described a more realistic unit ceil geometry, assuming a cylindrical shape for yams. 
According to their analysis, at the yam jamming point, the yam orientation angle has a 
maximum value of 55° and the yam volume fraction has a maximum value of 0.685. 
Assuming a fiber packing fraction of 0.785, this means that maximum fiber volume 
fraction for the track and column braid is 0.538. At low braiding angles (<20°), the 
geometric model by Li, Hammad and El-Shiekh [5.3] predicts a fiber volume fraction of 
less than 0.328 (also assuming a fiber packing fraction of 0.785). 
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In this report, we propose a raicrogeometric model for the track and column braid based on 
experimental observations and computer simulation. The unit cell geometry has been 
defined to establish the relationship of geometric parameters and processing variables. 

A summary of braiding process is descried in the following prior to the unit cell modeling. 
As described in Chapter 2, Figure 5.1(a) shows a basic loom setups in a rectangular 
configuration. The earners are arranged in tracks and columns to form the required shape 
and additional carriers are added to the outside of the array in alternating locations. Four 
steps of motion are imposed to the tracks and columns during a complete braiding machine 
cycle, resulting in the alternate X and Y displacement of yam carriers, as shown in Figure 
5.1(b)-(e). Since the track and column both move one carrier displacement in each step, 
the braiding pattern is referred to as lxl. 

5.2.1 Assumptions and Nomenclature 

Following assumptions are made to simplify the geometric model of 3-D track-and-column 
braided structures: 

1. No axial yams (0° insertion) are included in the 3-D structure. Axial yams may be 
used to increase composite modulus and strength in longitudinal direction, although 
it is not popular in practice due to the fact that 3-D braids usually have low 
orientation angle as fabricated and hence the strength increase is not phenomenal. 

2. lxl braiding pattern is assumed. Other patterns such as 1x2 and 2x2 can also be 
used but are not popular in practice. The techniques used for lxl can easily be 
extended to other patterns. 

3. Braider yams have circular cross-sections, same linear density and constant fiber 
packing fraction: 

4. Yam tensions during braiding are sufficiently high and hence the yam crimping 
effects are negligible. 
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St e P 3 Step 4 


-igure 5. : - Formation of a rectangular 3-D track 
and column oraid 


78 


The definitions of the symbols used in our analysis are listed below: 

Ac area of braid cross-section vertical to braid axis 
Af area of total fibers in braid cross-section vertical to braid axis 
A y area of braider yam in braid cross-section vertical to braid axis 
b dimension of yam cross-section in x'-z and y'-z planes 
b ' dimension of yam cross-section in x'-y' plane 
d yam diameter 

h z pitch length of braid formed in a machine cycle (four braiding steps) 

Ny number of braider yams 

Vf fiber volume fraction - all fibers to total volume 

r| braid tightness factor 

k yam packing fraction (fiber-to-yam area ratio) 

9 angle of braider yam to braid axis (yam orientation angle) 

Subscripts z, x', y' all refers to co-ordinates. 

5.2.2 Unit Cell Geometry 

The traditional approach used in modeling 3-D braided composites is to artificially define a 
unit cell geometry tor a 3-D braided structure without providing any relationship between 
processing variables and geometric parameters [5.4], All fibers in the unit cell are assumed 
to incline in 4 different diagonal directions, as well as along the longitudinal direction, if 
any. Fiber volume fraction is assumed to be either known or measured. In this work the 
dimension, shape and fiber architecture of the unit cell is based on process and structural 
analysis. Once the unit cell is identified the relationship between processing variables and 
key geometric parameters has been established. 


The key geometric parameters of 3-D braids (which affect reinforcement capability and 
composite processability) include braider orientation, total fiber volume fraction, volume 
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fraction of inter-yam void and axial fiber percentage of total fibers. Although there are only 
two simple process parameters adjustable to control the micro-structure of 3-D braids 
(speed ratio between braiding and take-up and linear density ratio of braider and axial 
yams), the process-structure model of 3-D braid is complicated. 

Normally, yam bundles consisting of numerous continuous filaments are used for fabric 
preforms, thus, the fabric raicrostructure has three levels: geometry of interfiber packing in 
the yam bundle (fiber level), cross-section of yam bundles in the fabric (yam level) and 
orientation and distribution of fibers in the 3-D network (fabric level). The unit-cell 
technique is commonly used to establish the geometric relation. In most of 2-D fabrics a 
unit cell geometry is readily identified, but in complex 3-D fabrics it can be very difficult to 
define. 

The fiber volume fraction of a 3-D fabric depends on the level to which yams pack against 
each other in the structure and the level to which fibers pack against each other in a yam. 
Two basic idealized packing forms can be identified: open-packing, in which the fibers are 
arranged in concentric layers, as illustrated in Figures 5.2(a) and 5.2(c); and close-packing, 
in which the fibers are arranged in a hexagonal pattern as in Figures 5.2(b) and 5.2(d). 
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Figure 5.2 Fiber pacldng in yams, (a) Open-packing in circular 
yams, (b) Close-packing in hexagonal yarns; (c) Open-packing in 
diamond-shaped yam; (d) Close-packing in diamond-shaped yam. 
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In addition to the level of packing fraction, the fibers also establish the yam cross-sectional 
shape, i.e., yam packing in fabrics. This shape plays a very significant role in determining 
how many fibers can be packed into a fabric. One of good examples is the yam packing in 
2-step braided preforms [5.5]. Due to the use of untwisted fiber bundles and high braiding 
tensions, cross-section of axial yams in the 2-step braid is deformed to prismatic shapes, 
giving most compact yam packing within the braided structure. For the track-and-column 
braids, the braiding tensions are lower compared to the 2-step braids and the cross-sections 
of yams actually have a polygonal shape [5.6]. For reasons of simplicity, we idealize the 
polygonal yam cross-section as circular shape. 

In order to understand the braid internal structure and the yam interlacing pattern, analysis 
of braiding carrier motion [5.6], computer graphics simulation [5.7] and computer solid 
modeling [5.8] was performed and experimental observations of cross-sections of a 
CARBON-PEEK braided composite[5.6] were made. Figure 5.3 shows an idealized braid 
cross-section cut longitudinally at a 45° angle to the braid surface. There are four groups of 
yams inclined at angle a with the braid axis (z direction) in different directions; the yams in 
each group are parallel to each other within a specific plane. Two groups of yams are 
parallel to the x'-z plane; the other two are parallel to the y'-z plane. The cutting plane is so 
selected that it cuts through the diameter of a group of yams. 

Figure 5.4(a) shows the unit cell identified from the analysis. The unit cell consists of four 
partial yams being cut by six planes. Clearly, there does not exist such a unit cell which 
only consists of four complete yams. The dimensions of the unit cell are (l/2)h x ’ in x' 
direction. (l/2)hy in y’ direction and (l/2)h z in z direction (braid length). The cross-sections 
of the unit cell at (l/2)h z , (3/8)h z , (l/4)h z , (l/8)h z and 0, are shown in Figures 5.4(b) - (f), 
respectively. 
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Figure 5.3 Braid cross- section cut longitudinally at a 45° angle to 
the braid surface by the x'-z plane ABCD. z is the braid length 
direction. 
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Figure 5.3 Braid cross-section cut longitudinally at a 45° angle to 
the braid surface by the x’-z plane ABCD. z is the braid length 
direction. 
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Figure 5.4 Unit cell geometry of 3-D track and column braid, 
(a) unit cell, (b) unit cell cross-section at z = 1/2 hz. (c) unit 
cell cross-section at z = 3/8 hz.(d) unit cell cross-section at z = 
1/4 hz. (e) unit cell cross-section at z = 1/8 hz. (f) unit cell 
cross-section at z = 0. 
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5.2.3 Unit Call Model 


The portion of O 1 O 2 O 3 O 4 in Figure 5.3 is extracted and enlarged in Figure 5.5 to show 
geometric details. In Figure 5.5, d is the yam diameter, 6 is the braiding angle. 
Dimensions of b, b' (see Figure 5.4) and h z can be obtained from the trigonometric 


relationships: 


b = * 
sin 0 

(5-1) 

b’ = — 

COS0 

(5-2) 

Vl + cos 2 0 
h z — 2 d , 

sin 0 

(5-3) 


h z is actually the pitch length of braid formed in a complete machine cycle (four steps). 
This length is one of the key parameters in controlling the fabric microstructures. Using 
Equation (3), the braiding angle 0 can be expressed in terms of the braid pitch length h z : 


9 ^ W (hz/d) 2 +4 (h z >2d) 

Restriction h z > 2d was applied to the above equation to ensure that O<0 <90°. The physical 
meaning of this restriction is that the braid pitch length must be greater than two diameters 
of the yam to maintain a stable convergence point during braiding process. 

As can be seen in Figure 5.4(b)-(f), each unit cell cross-section consists of four half oval 
cross-sections of yam. Based on the assumption of non-crimp yam paths, total yam area 
in every braid cross-section must be the same. Therefore, a conclusion can be reached that 
any unit cell cross-section contains four half oval cross-sections of yam. The fiber volume 
fraction can then be easily derived as: 


2(k)(|b')7C d 2 

Vr = — K = 9 n k - 

(|hx)(|hy) " hx' hy COS0 


(5-4) 
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Figure 5.5 A portion of braid cross-section (01020304) in 
the x'-z plane consisting of four unit cells. 
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where k is the fiber packing fraction, which normally has a value of about 0.785 for 3-D 
braids after matrix addition and consolidation [5.5]. 

The braid has the tightest structure when each yam is in contact with all its neighboring 
yams, in other words, the yams are jammed against each other. Under this jamming 
condition, the unit cell has a minimum dimension in both x' and y’ directions, i.e.: 

, V 1 + cos 2 0 

hy min = hx min = 2 d (5-5) 

J COS0 

and the braid has a maximum fiber volume fraction, which can be derived by combining 
Equations (4) and (5): 


7C COS0 

Vftnax = 2 K 1+cos 2 0 


(5-6) 


Due to bulky fibers and nonlinear crimp nature, it is difficult to fabricate a braid having the 
tightest structure. In practice, the fiber volume fraction is calculated directly from the area 
of fibers to the area of braided composite in the cross-section perpendicular to braid axis: 



(5-7) 


where Ac is the area of the braided composite and Af is the area of total fibers given by: 


Af = K 


NvAv 

COS0 


(5-8) 


Combining Equations (5-7) and (5-8), we have: 


Vf= K 


NyAv 

A c cos0 


(5-9) 


In order to measure the fiber compactness in the braided structure, we define a braid 
tightness factor, t|, as: 


0 = 


Ny Ay 
A c 


(5-10) 
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According to its definition, the braid tightness factor is actually the ratio of total yam area to 
the braid cross-section area in x'-y’ plane when 0 = 0°. Using Equation (5-10), the fiber 
volume fraction can be simplified to: 


V, = K 


_n_ 

COS0 


(5-11) 


Clearly, the amount of fibers within the braid is limited by the maximum fiber volume 
fraction Vfn.»». Combining Equations (5-6), (5-9) and (5-1 1), we obtain: 


ri it cos0 

Vf = K < -Z K 

COS0 2 l+cos 2 0 


(5-12) 


7C COS 2 0 
^ - 2 l+cos 2 0 


(5-13) 


Based on the geometric analysis given above, braiding angle 0 is determined from yam 
diameter d and braid pitch length h z . The fiber volume fraction, Vf, is controlled by 
braiding angle 0 and braid tightness factor r\. The tightness factor h, in turn, has to be 
properly selected so that the required fiber volume fraction Vf is achieved and over- 
jamming is avoided. Figure 5.6 shows the Vf-e relauonship based on Equation (5-12). 

The fiber packing fraction is assumed to be 0.785 in our calculation. A maximum fiber 
volume fraction of 0.617 is achieved at 0 = 0° (i.e. the "braid" is an assembly of 
longitudinal yams without interlacing). At 0 = 90°, which is impossible to achieve either 
theoretically or in practice, fiber volume fraction approaches zero. There is no such 
maximum braiding angle as suggested by Li, Hammad and El-Shiekh [5.2], wherein 3-D 
track-and-column braids can be made with O<0 <90°. As 0 becomes higher, the yams 
slide away from each other to form a more open structure, resulting in a lower fiber volume 
fraction. In practice, a braiding angle higher than 45° is difficult to achieve because of yam 
slippage creating difficulties in maintaining a constant height of convergence point 
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Figure 5.6 Relationship of fiber volume fraction to braiding 
angle for various tightness factors. Fiber packing fraction k is 
assumed to be 0.785. 


89 



Although a wide range of fiber volume fractions can be achieved for 3-D track-and-column 
braids by using different levels of the braid tightness factor, the actual fiber volume fraction 
is more likely to be closer to the jamming region as shown in Figure 5.6. The reason for 
this is that yarns are usually tightly compacted by their neighboring yams due to the high 
yam tensions applied to overcome inter-yam friction during braid formation. In practice, it 
is possible that the fiber volume fraction achieved will be higher than V fm«« under jamming 
conditions because of non-linear yam paths (yam crimp) in the braided structure. For 
example, at 0 = 20°, fiber volume fraction under jamming conditions is 0.615, while fiber 
volume fractions in composites having a braiding angle of 20° observed from experiments 
normally fall within a range of 0.60 to 0.65. 
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5.3 Multlaxial Warp Knit 

A series of studies on the technology, structure, and properties of the MWK preforms and 
composites have been reported by Ko and his co-workers [5.9-13]. In this section, we 
propose a unit cell model for the four-layer MWK structure. With minor modifications, the 
analysis can be generalized into the MWK system with six or more layers of insertion 
yams. 

The MWK fabric preforms having four directional reinforcements similar to quasi-isotropic 
lay-up can be produced in a single step. The key geometric parameters of the MWK fabric 
preforms, which affect the reinforcement capability and the composite processability, 
include the number of yam axis, the orientation of bias yams, total fiber volume fraction, 
pore size and pore distribution, and percentage of stitch fibers to total fiber volume. The 
processing variables adjustable to control the MWK micro-structure include the type of knit 
stitch, the ratio of stitch-to-insertion yam linear density, the orientation angle of bias yams, 
and the thread count The concept of a unit-cell is used to establish the relationship between 
the geometric parameters and process variables. 

5.3.1 Unit Cell Modeling 

As demonstrated in the development of 3-D braid model, the first step in the unit cell based 
modeling is to determine the unit cell dimension, so that it is the smallest repeating unit of 
the structure. The second step is to assume some idealized cross-section shapes of the yam 
bundles, based on experimental observations. The final and most important step is to 
identify the overall unit cell geometry, from which expressions for the key geometric 
parameters can be derived, and the geometric limits applied to the structure can be defined. 

As can be seen in Figure 5.7, the unit cell for the MWK fabric preform can be defined in 
many ways to meet the definition for a unit cell. However, it would be most reasonable to 
have a unit cell which consists of a complete knitting stitch, and the insertion yams in the 
unit cell are all symmetrical. Figures 5.8(a) and (b) show the unit cells for the MWK 
structure with chain and tricot stitch, respectively. Within the outlined in the figure, the unit 
cell consists of one each of 0° and 90° yams, two each of +8 and -0 yams as well as a knit 
stitch. 

In order to determine the shape of yam cross-sections, a MWK reinforced carbon/epoxy 
composite was cut perpendicularly to its 0° axis, polished and examined using an optical 
microscope. Figure 5.9 shows the photomicrograph of a cross-section of the MWK 
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composite. As can be seen in Figure 5.9, the insertion yams have a rectangular cross- 
sectional shape with rounded comers. For example, the 0° yams have a width-to- thickness 
aspect ratio of 5.7±0.42. 

Based on this observation, the insertion yams are assumed to have a race-track cross 
section having a width-to-thickness aspect ratio / > 1, as illustrated in Figure 5.10(a). The 
stitch yam, which has a much lower linear density than the insertion yams and usually 
contains some twist, has less tendency to spread out in the knitting process. Therefore, a 
circular cross-section with a width-to-thickness aspect ratio / =1 is assumed for the stitch 
yam, as illustrated in Figure 5. 10(b). 
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(b) Tricot stitch 


Figure 5.7 Multiaxial warp knit (MWK) with four 
layers (0°, 90°, and ±9) of insertion yams and (a) chain 
stitch or (b) tricot stitch. 
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(a) Chain stitch 



(b) Tricot stitch 


Figure 5.8 Unit cell which consists of one each of 0° and 
90° yams, two each of +9 and -0 yams as well as a knit 
stitch: (a) with chain stitch; (b) with tricot stitch. 
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Figure 5.9 Photograph of a cross-section of MWK 
reinforced carbon/epoxy composite. 
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(a) / = w/t > 1 


(b) / = 1 


Figure 5.10 Shape and dimension of yam cross-section: (a) 
insertion yam is assumed to have a race-track cross-section 
with the width-to-thickness aspect ratio / > 1; Ob) stitch 
yam is assumed circular in cross-section with / = 1. 
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Figure 5.1 1 shows the idealized unit cell geometry for the MWK structure, including the 
shape, dimension, orientation, and position of all the insertion and stitch yams within the 
fiber 3-D network. The knit stitch is assumed to have the tightest loop construction, and 
the curved loop is idealized to a rectangular shape, as illustrated in Figure 5.1 1(b) and 
Figure 5.1 1(d). The dimensions of the unit cell are X, Y, and Z, respectively, 
corresponding to the 0° axis, 90° axis, and the thickness axis vertical to the 0°-90° plane, as 
shown in Figure 5.11 (a) and Figure 5.11 (c). 



Figure 5.11 Unit cell geometry of MWK structure: (a) with 
insertion and chain stitch yams; (b) chain stitch yams only; (c) 
with insertion and tricot stitch yams; (d) tricot stitch yams only. 


97 





Based on the i dealize d yam shape and unit cell geometry, the geometric model of the MWK 
fabric preform can be established. In the following section, the expressions for yam 
dimensions, unit cell dimensions, yam volumes, yam lengths, fiber volume fractions, and 
criteria for geometric limit will be presented without derivation. 

5.3.2 Nomenclature 

k fiber packing fraction (fiber-to-yam area ratio) 

p fiber density (kg nr 3 ) 

X yam linear density (kg nr 1 ) 

9 angle of bias yams to 0° yam (x direction) 

ds diameter of stitch yam 

w, t width and thickness of insertion yams 
/ aspect ratio of insertion yams (width to thickness) 

S area of yam cross-section (m 2 ) 

L yam length (mi 

1 yam length normalized with 0° yam length (X) 

X, Y, Z unit cell dimensions in height, width, and thickness directions, respectively 
4> thread count (number of insertion yam per unit length, end nr‘) 

n fabric tightness factor (as defined by Eq.t5-44)) 

V volume (m 3 ) 

Vt fiber volume fraction 

v percentage of insertion or stitch fibers to total fiber volume (%) 

Subscripts x. y, Z, and ±9 ail refer to yam orientations in the unit cell frame 
Subscripts c. i, and s refer to unit cell, insertion yam and stitch yam, respectively 
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5.3.3 Dimensions of insertion and stitch yarns. 

The widths of insertion yams: 

w x = /x t x 

Wy = /y ty 
W±Q = /±0 • t±0 

The thicknesses of insertion yams: 




The diameter of stitch yam: 
ds = 

The cross-sectional areas of insertion and sdtch yams: 
Sx = (4 + / x - 1 )• U 2 
Sy = (4 + / y ‘ 1 )■ *y 2 
S±e = (4 + /±e • 1 )■ t±9 2 
S s = J - d s 2 



( 5 - 14 ) 

( 5 - 15 ) 

( 5 - 16 ) 

( 5 - 17 ) 

( 5 - 18 ) 

( 5 - 19 ) 

( 5 - 20 ) 

( 5 - 21 ) 

( 5 - 22 ) 

( 5 - 23 ) 

( 5 - 24 ) 
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5.3.4 Dimensions of unit ceil. 


The dimensions of the unit ceil in x (0°), y (90°) and z (thickness) directions: 
X = <|>x Wy 

Y = <|>yWx 

Z = tx + ty + 2 (t+Q + dg) 

The orientation angle of bias yams: 

0 = tan- 1 

The total volume of the unit cell: 

V c = X Y Z 



5.3.5 Volumes of yarns in the unit cell. 

The volumes of the insertion yams in x, y and =9 directions: 

V ix = S x X 
Vjy = Sy Y 

Vi±e = 2 S ±9 \'X2 + Y2 


The total volume of the insertion yams in the unit cell: 

Vj = Vjx + Vjy t Vj+9 

The volumes of the chain sutch yams in x. >• and ±9 directions: 

V$x = 3 S$ • (X - g ds) 

Vsy — 2 Ss d$ 


(5-25) 

(5-26) 

(5-27) 

(5-28) 

(5-29) 


(5-30) 

(5-31) 

(5-32) 

(5-33) 

(5-34-a) 

(5-35-a) 

(5-36-a) 


V sz = 2 S s (Z - d s ) 

V s ±9 = 0 
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(5-37-a) 



The volumes of the tricot stitch yams in x, y and ±9 directions: 


V sx = 2 S$ • (X - g ds) (5-34-b) 

V sy = 2 Ss • ds (5-35-b) 

V sz = 2 S s (Z»d s ) (5-36-b) 

V s ±e = S s (VX2 + Y2 - d s ) (5-37-b) 

The total volumes of the stitch yams in the unit cell: 

V s = V sx + V sy + V sz + V s±e (5-37) 


5.3.6 Normalized yarn lengths in the unit ceil. 

The lengths of insertion yams normalized with the unit cell dimension X: 
X 

1 ix — y — 


liy = tan 6 


lj+e = 4 V 1 + tan 2 e 


(5-39) 

(5-40) 

(5-41) 


The length of chain stitch yam normalized with the unit cell dimension X: 

„ „ 3 ds -r Z 

lc = 3 + 2 w 


(5-42-a) 


The length of chain stitch yam normalized with the unit cell dimension X: 

. 3 -d s + Z ; — 

l s = 2 + 2 y + \ 1 + tan^Q 


(5-42-b) 


5.3.7 Fiber volume fractions. 

The volume fractions of total insertion and stitch fibers: 


101 



(5-43) 


v « - v c K ' 

Vfs = 


The volume fractions of all fibers in x, y, z, and ±9 directions: 


Vf z = 


V f ±e = 


Kj 

V jx + Ks Vs x 


Vc 

Ki 

' V jy + Ks • V s y 


Vc 

K S 

• V S z 


Vc 

Kj 

• V j±0 + Ks • Vs±9 


Vr 


The overall fiber volume traction of the MWK fabric preforms: 


Vf = Vfj + Vf s , or 
Vf = Vf x + Vfy -p Vfz + Vf+e 


(5-44) 

(5-45) 

(5-46) 

(5-47) 

(5-48) 


(5-49-a) 

(5-49-b) 


Percentages of insertion and stitch fibers to total fiber volume: 


Vfj 


-100% 

1 Vf 

(5-50) 

= . . * • 1 00% 
s Vf 

(5-51) 


5.3.8 Geometric Limitation 

With the geometric model established, one can now proceed to examine the geometric limits 
of the MWK structure. One important factor to be considered in the design and 
manufacturing of fabrics is the yam jamming condition in fabrics. Yam jamming condition 
is a geometric boundary condition under which all the yams are touching each other, 
resulting in the tightest possible structure. The boundary condition which defines the 
limiting geometry of the MWK fabric can be summarized by Eqs (5-52) to (5-54): 
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(5-52) 


X ^ Wy + ds 

Y > w x + 2 d s (5-53) 

X • sine + d s (1 + cose) > w ±0 (5-54) 


Combining with Eqs (25) and (26), Eqs (52) to (54) become: 

f #4 in . * _ W / 1 i o p Q \ 


4» x > max. 

- 1 + w 


d s w ±e - d s (1 + cose) 

+ Wy i Wy sine 


v 

2 • d s 


(5-55) 

(5-56) 


Eqs (55) and (56) serve as the criterion for yam jamming, which not only gives the 
boundary condition for the geometric model, but also defmes the important processing 
limit. In order to measure the degree of the fabric tightness, a processing parameter - the 
tightness factor r\ is introduced: 


d s 

0 

W v 

n = — , *1 


(5-57) 


Clearly, the higher is the tightness factor, the more compact the insertion yams are within 
the fabric structure. When n = 1, the unit cell has a minimum dimension in X direction, and 
the fabric has the tightest structure at an given angle of bias yams. After determining the 
value of the tightness factor r|, the orientation angie of the bias yams (0), and the ratio of 
stitch-to-insertion yam linear density OW-M) must be so selected that both Eqs (5-55) and 
(5-56) are satisfied. 


5.3.9 Results and Discussion 

Using the equations developed in the analysis, we have established the key relations 
between process variables and fabric geometry. The results of the geometric modeling can 
be represented in many ways, and the selection of relations to be evaluated depends on the 
specific technical need. We have focused on the overall fiber volume fraction (Vf), the 
percentage of stitch fibers (v s ), and the geometric limit to describe the MWK structure. The 
key process variables include the orientation angle of bias yams (9), the ratio of stitch-to- 
insertion yam linear density (W^)> and the fabric tightness factor (n). 
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The re l ati o nship between these geometric parameters and process variables is shown in 
Figures 5.12 -5.14, using Eqs (5-14) to (5-57). In our calculations, we assume that tricot 
stitch is used, fiber material and fiber packing fraction are the same for all insertion and 
stitch yams, i.e., k = iq = k$ ^d p = pj = ps; and all the insertion yams have the same 
linear density and width-to- thickness aspect ratio, i.e., Xj = X« = *iy = Xj±e and f \ = f \ x = 

f\y — f i±fl- 


5.3.10 Effect of yarn linear density ratio on fiber volume fraction. 

The fiber volume fraction relation in Figure 5. 12 shows that for the fixed parameters 
selected, only a limited window exists for the MWK fabric construction. The window is 
bounded by two factors: yam jamming and the point of 90° bias yam angle. Fabric 
constructions corresponding to the curve marked "jamming” are at their tightest allowable 
point, and constructions at the 0 = 90° curve have the most open structure. When 0< 30°, 
jamming occurs in the whole range of yam linear density ratio from zero to infinite. When 
9 is in the range of 30° to 40°. the fiber volume fraction decrease with the increase in yam 
linear density ratio until jamming occurs. When 0 = 45°, the fiber volume fraction 
decreases with the increase in yam linear density ratio to a minimum at about l</lj =1, and 
starts to increase until jamming occurs. When 9 > 60°, the fiber volume fraction has the 
same trend as when 0 = 45°, but yam jamming never occurs. The fiber packing in the 
yams, taken as 0.75, limits the maximum fiber volume fraction in the fabric. 

5.3.11 Effect of fiber orientation on fiber volume fraction. 

The fiber volume fraction is plotted in Figure 5.13 against the orientation angle of bias 
yams for seven levels of the fabric tightness factor. It can be seen that for a given tightness 
factor, the fiber volume fraction decreases as the bias yam angle increases. At t| > 1, yam 
jamming always occurs. At n = 1 and 0 = 50°. the fabric has the tightest structure 
(maximum fiber volume fraction) for the parameters indicated in the caption. At lower 
levels of t|, smaller bias yam angles are possible, but the fiber volume fractions also 
become smaller. 

5.3.12 Effect of yarn linear density ratio on fractional stitch fiber volume. 

Figure 5.14 shows the percentage of tricot stitch fibers as a function of the ratio of stitch- 
to-insertion yam linear density. The possible range of bias yam angle is from 30° to less 
than 90°. As can be seen in Figure 5.15, the stitch fiber percentage increases as the yam 
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linear density ratio increases, but decreases as the bias yam angle increases. The process 
window of v s is bounded by two curves as shown in Figure 5.14. The upper bound is 
given by the yam jamming limit, and the lower bound is reached at 0 =90°. This provides a 
guide to determine the yam linear density ratio to achieve a required stitch fiber percentage 
for a given bias yam angle. For example, if we want to fabricate a MWK fabric preform 
which has a 10% of stitch fibers and 45° of bias yam angle with other parameter being the 
same as given in the caption of Fig.5.2. the ratio of stitch- to-insertion yam linear density 
required is about 0.15 as can be read direcdy from Figure 5.14. 
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Figure 5.12 Fiber volume fraction versus ratio of 
stitch-to-insertion yam linear density (tricot stitch, 
k = 0.75, p = 2.5 kg m-3, /i = 5, and r| = 0.5). 
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Figure 5.13 Fiber volume fraction versus bias yarn 
orientation (tricot stitch, k = 0.75, p = 2.5 kg m-3, 
fi = 5, and ls/li = 0.1). 


107 







Figure 5.14 Percentage of stitch fibers versus ratio of 
stitch-to-insemon yam linear density (tricot stitch, k = 
0.75, p = 2.5 kg m-3, fi = 5, and r| = 0.5). 
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Chapter - 6 Applications 

6.1 Implementation of the Geometric Model for pre-form design 

The geometric model is implemented on a Sun 3/160 work station using C language. The 
computer code has a modular structure, each module performs a clearly defined function. 
New modules can be added to increase the capabilities of the code. Current capacity of the 
code include the computer simulation of the 3D braided structure (3DB) and the Multiaxial 
Warp Knit (MWK) structure. 

The simulation is divided into stages. In stage 1 , a group of modules, which were compiled 
into an executable file named preform generator, is executed. The processing and/or 
geometric parameters are input from keyboard. Based on the logic described in earlier 
chapters, the yam path is determined and a device independent, neutral graphics file is 
written. In stage two a model depictor converts this neutral data files into screen images. 
The Sun graphics library calls are used to draw the graphics primitives. The logic flow of 
different modules are given below. The stage 1 is device independent The data files after 
stage one are written in ASCII text so that the files are transferable from one computer 
system to the other and all the modules need not to be run on the same machine. 
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6.1.1 The Preform Generator 


Current capabilities for preform generation include the 3D braid and the MWK. The 
computer codes to perform the preform generation were developed in the C computer 
language. C is a powerful language and is available in practically all unix based 
workstations. The codes are portable with minimal effort and hence are hardware 
independent. The basic principle of preform generation is to identify a number of key points 
along the yam trace and then to pass a yam through the keypoints. 

Preform gen eration for 3D braid 

The keypoints for the 3D braided preforms are generated from the machine motion. The 
yam trace is determined by taking keypoints along the machine motion and then passing a 
B-spline along the keypoints. The yam path is simulated such that it assumes a shape that 
minimizes the energy of the system. To accomplish this additional keypoints along the yam 
center line is generated. These keypoints through which the B-spline passes are called knots 
in spline vernacular. The details of the yam path modeling are given in Chapter 3. 

Once the yam trace is established the yam cross-section is taken as a input and it is swept 
along the trace to generate the solid yam. The twist of the yam at the turning point is 
calculated. A neutral, device independent graphics file is generated for later rendering by a 
depiction program. 

The kevpoint generation for 3-D braided preform is accomplished by 3dbraid. The default 
name for the output the intermediate data file is movementdat but it can always be 
renamed. Once the yam-path is generated the solid yam is built along the path by calling the 
module build_yam. The graphics output of this stage is written in neutral graphics file. 
The default name for the device independent graphics file is n.dat but is usually renamed to 
represent the geometry it renders. 

Preform ge neration for MWK 

The kevpoint trace for the MWK structure follows the same principle as 3-D braids, 
although the actual implementation involves a few more steps. The added complexity arises 
from the fact that the MWK structure consists of both insertion and stitch yams that have to 
be considered separately. The keypoints for the stitch yams are generated by mwk. The 
keypoints for the insertion yams are selected by observation. The body around the trace of 
insertion yam is generated by the program xybuild, the body around the trace of the stitch 
yam is built by build_yarn (same program as used for the 3D braid). The theory of 
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selecting the keypoints is given in Chapter 3 and the details of the implementation is given in 
the user manual section of this report 

6.1.3 The Model Depictor 

The neutral graphics data file generated by the model generators 3dbraid (for 3DB) or 
mwk ( for MWK) in conjunction with build .yarn and/or xybuild, is rendered on the 
computer screen at this stage. The depiction program used is called render. The depiction 
program is more than just a image renderer, it can also manipulate the image by rotating it 
for different views and produce section diagrams by passing a cutting plane through the 
structure. 

The rendering program reads in the neutral graphics file and depicts it on the screen. While 
the rendering logic is similar for all devices the actual subroutine calls are particular to the 
hardware used. In the current project the hardware used to render the model is a Sun 3/160 
workstation. The workstation is equipped with a Graphics Processor (GP). Sun 3/160 is 
equipped with a 8 bit color display. The value of each pixel serves as an index to a 
colormap table. This techniques allows the mapping of up to 16 million colors. At anyone 
time the 8 bit monitor can however display up to 256 colors simultaneously. 

The rendering of the neutral graphics file is accomplished by fully utilizing the SUN 
systems Sun View GUT (graphics user interface). SunView is Sun Microsystems 
proprietary system. It supports interactive graphics based applications running within 
window environment. The driving engine for the renderer is render. The use of the 
modules 3dbraid. mwk. build .yarn and render are described in section 6.1.4. The 
code is listed in Appendix A. 

6.1.4 Installation of Software on a SUN workstation. 

The software is available on a SUN QIC 24 9 track data cartridge in tar format To install 
the software on to a SUN workstation the cartridge is to be placed in the drive and the files 
to be extracted using the following command: 

Car -xvf /rst0/nasa3d. tar 

rstO is the drive name and may vary from system to system. This will automatically create 
the necessary subdirectories and install the code. 
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The program is supplied as executables that will run on a SUN 3/160 workstation. To 
obtain copies of the program contact FMRC at the address on the cover page. 


6.1.5 Using the modeling software - User Instructions 
Overview 

The Simulation takes the following three steps to complete: 

Step 1: Depending on which preform structure (either the 3d braid or the Multiaxial warp 
knit) is to be depicted, the modules 3dbraid or mwk is invoked to generate yam 
keypoints. The default name of the intermediate output after this step is the ASCII text file 
movement.dat for 3dbraid and outfil for MWK. 

Step 2: The solid yam is built around the yam trace using the module build_yam and/or 
xybuiid. The default name of the output file from this step is n.dat. 

Step 3: The neutral graphics file generated is then rendered on the screen by the module 

render. 

In the rest of the section step by step user instruction to generate yam keypoints for a simple 
3d braided structure and a simple Multiaxial warp knit structure is given. In the following 
section all the user commands are described in detail for use as a reference. 

Preform Generation for 3d braided structure - User Instruction 
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1. After the software is installed, the model generation is started by typing 
3dbraid <Retnm> 

The user will be prompted for the number of tracks and number of column information. The 
queries and input are shown below. The user input is bolded. 

How many tracks? 4 <Retum> 

How many columns? 5 <Return> 

Once the track and column information is input the user is prompted for the number of 
movements to be simulated. For 4-step braiding, 4 movements complete a cycle. The final 
prompt is for takeup distance. 

How many movement?2 <Retum> 

Distance? 1 <Enter or Retum> 

The input is checked for loom size. Currently the maximum allowable loom size is 50x50. If 
either track or column exceeds the allowable value the input is rejected and the user is 
prompted again. 

Loom too big! # of track or # of column should be less than 50. 

Try again! 

How many tracks? 

The data file generated by this step is named movementdat. The data file for the example 
inputs shown is appended in Appendix A2. The trace of the yam path is recorded in 

movement.dat . 

2. Once the movementdat is generated , the program build yam is used to place a 
solid yam along the trace. Proper yam twist at changeover points is also incorporated at this 
point The task is started by typing 

build _yam <retum> 

The user is prompted for the geometry of the yam. The cross-section is modeled as an 
ellipse and the user is input for the major axis and the minor axis. For a circular cross- 
section, the user needs to input the radius twice. .An example session is shown below. 
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buikLyarn moveraenLdat<retura> 

Yam radial 1? 0.5<retum> 

Yam radial 2? 0.25<return> 

The output is a neutral graphics file in text format and is automatically named n.dat 
(included in Appendix A2). 

It is important to note that steps 1 and 2 do not necessarily has to be carried on the Sun. The 
user can run the program at his computer of choice and transfer the output n.dat to the Sun 
for graphical rendering. 

Preform Generation for MWK - User Instruction 

The preform generation for the MWK follows similar principles as the 3Dbraid. 

1. The first step is to generate the keypoints for the stitch yam. This is done by typing 
mwk. On starting the mwk the user is taken through a series of queries for the processing 
parameter. A sample query is shown below. 

mwk <retum> 

Enter # of insertion yam/unit length 
2 <retum> 

Enter the angle (degree) of bias yams 
40 <retum> 

Enter the radius of the insertion yam 
.4 <retum> 

Enter the radius of the stitch yam 
.125 <retum> 

Enter the radius of the stitch loop 
.2 <retum> 

color of the stitch yam: 1, 2, 3, 4 
1 <return> 

2. Once all the necessary input is given to the computer an intermediate outfil is generated. 
This outfil records all the keypoints along the trace of the stitch yams. 
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3. Next step is to build the trace of the insertion yams. The insertion yams are all arranged 
along straight lines and hence the two end points of each yam are taken as the keypoints. 
The user needs to know the number of insertion yams he or she needs to plot and then 
calculate the coordinates of the end points of the trace. The user then edits a small file, 
outfiI2, with all the keypoints. A sample outfil2, file is shown below. 

9 120 
1 0 

0.0 0.5 0.225 

1.0 0.5 0.225 
1 0 

0.0 0.5 0.225 

1.0 0.5 0.225 
1 0 

0.0 1.5 0.225 

1.0 1.5 0.225 
2 90 

0.0 0.0 0.625 
0.0 2.0 0.625 

2 90 

1.0 0.0 0.625 

1.0 2.0 0.625 

3 45 

0.0 0.0 1.025 

1.0 1.0 1.025 

3 45 

0.0 1.0 1.025 

1.0 2.0 1.025 

4 -45 

1.0 0.0 1.425 
0.0 1.0 1.425 
4 -45 

1.0 1.0 1.425 
0.0 2.0 1.425 

The data is entered in free format The first line has two elements. The first element denotes 
the total number of yams to be plotted the second number denotes the number of sections 
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each yam is to have. The more the number of section the smoother is the yam but at the 
expense of longer computer time. After the first line there is a three line block of data for 
each line. The first line in each block denotes the yam color and the second number the 
orientation of the yam. The remaining two lines give the starting and end point coordinates. 
In the sample file the data is prepared for nine insertion yams. 

4. The body is built now around the yam traces. The program build_yarn is used to build 
around the stitch yams. The build_yarn is the same program as used to build the 3dbraid 
yams. A modified version of build_yarn, named the xybuild is used to build body 
around the trace of the insertion yams. 

The programs are started by typing 

build_yarn outfil <retum> 
mv n.dat m.dat 
xybuild outfil2 <retum> 

The user dialog is very similar to that described for build_yam in the section for 3Dbraid 
and the user is referred to that section. The default output from the code is n.dat . So care 
should be taken to rename the n.dat generated by build_yam before xybuild is run or the 
original n.dat will be overwritten. 

Once the two data files m.dat and n.dat are generated they are combined to one file. Each 
of the files have two header lines which contain important information that is needed for 
successful combination. A sample header for n.dat is shown below. 

9 6462 6480 

1.2000 - 0.2000 2.1414 - 0.1414 1.5982 0.0518 
718 720 

0.000000 0.700000 0.225000 

Each file contains one block of data for each yam. Each block can be identified by its 
header line which contains only two integers. In this example file the block header file starts 
with the integers 718 720. The first block from m.dat contains only one block of data. The 
first block of n.dat is replaced by the data from n.dat. The second line of n.dat is now 
compared to second line of m.dat and modified as explained below. The second line 
contains the information of the three-dimensional parallelepiped that encompasses the MWK 
structure. It is defined by inputting the range of x, y and z as shown: 
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imai xmin ymax ymin zmax zmin 
1.2000 - 0.2000 2.1414 - 0.1414 1.5982 0.0518 


The max. xmin , ymax, ymin, zmax and zmin of the two files are compared and the 
maximum value for max, ymax and zmax and the minimum values for min, ymin and 
zmin are chosen. These values are used to modify the second line of n.dat The neutral 
graphics file is now ready for rendering. 

Graphical Rendering 

The graphical rendering program is developed for a sun workstation environment The 
program is highly interactive. Iniual input is by Keyboard but once the image is drawn on 
the screen, further manipulation is done by the Sun Mouse. 

The optimum display is obtained when the Sun View window covers the whole screen. To 
ensure this the user has two choices. 

(a) The user can click at the corner of the active screen and select the full screen window 
option of the popup menu. Some useful Sun View window management techniques are 
discussed in Section 5. 

(b) The user can run the rendering program without starting SunView. In this case the 
rendering program starts SunView and makes the active window span the whole screen. 

To start the processing of n.dat neutral graphics file, type 


model n.dat <retum> 

input cut plane normal 0.1 0.1 0.1 <retum> 

At this point the control buttons for image manipulation pops up on the right side of the 
active window. The screen is shown in Fig. 6.3. 

After the program starts the image is manipulated by the control strip. Different image 
manipulation functions are performed by placing the cursor at the appropriate button and 
depressing the middle mouse button. A schematic of the Sun Mouse is shown in Fig. 6.4 
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Fig. 6.3 Screen Display during Simulation 



Fig. 6.4 - The sun three button mouse. 

In the rendering program only the left button is used. The right button is used by 
Sun View to change the size of the windows and to exit from Sun View. A schematic of the 
control strip is shown in Fig. 6.5. As can be seen the control strip offers numerous options. 
All these options are explained in the next section. In this session the minimal steps 
necessary to render the neutral graphics file is described. 
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Image 


Line 


X-Frame] 


Out Box 


Solid 


Reset 


View In; 


Larger 

m 


Smaller 


[Menu Hid 


Sw Direct 


Backwarc 


Cutting 


Forward 


Get Rast 


[ Move Rasj 
Load RasI 


Restore 


Cut Fil 


EXIT 


1. The first step is to select a bounding parallelepiped inside 
which the model will be drawn. The user points and clicks 
at the option Outbox. A wire frame parallelepiped is drawn 
on the screen. 

2. The user then controls the size of the parallelopiped by 
clicking the viewing buttons. The R, L, U and D 
buttons move the viewpoint to right, left, up and down 
respectively. The corresponding lowercase buttons perform 
the same function but at a slower rate and hence can be used 
to fine tune the position of the box. The Larger, Smaller, 
1 and s buttons are used to increase the size of the box. 
Larger and Smaller buttons increases or decreases the size at 
a higher rate than the 1 and s buttons respectively. 

3. After the box is placed the yam can be drawn by using 
one of the image buttons. The X-Frame button draws a 
wire frame image of the yam. The Solid option draws a 
hidden line solid model. 

4. The image can now be saved as a raster file by using the 
Move Rast button. 

5. To view a cross-section along a cutting plane the user 
uses the cutting in conjunction with the Forward and 
Backward buttons. While the cutting plane is 
predetermined, it can be moved along its normal by the 
Forward and Backward buttons. The position of the 
plane can be seen on the display. Once the cutting plane is 
positioned, the cross section is displayed by clicking on the 
cutting button. 


Fig. 6.5 

Schematic of the 
Control Strip 
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6.1.5 Command reference of the image manipulation controls 


The image manipulation commands are cataloged below. To select (or use) any of the 
commands the user needs to position the cursor in the appropriate box and press the left 
mouse button. 


Imagii 


Line 


X-Frame) 
Out Box 
Solid 
Reset 


The Image Buttons 

These buttons Control the Type of Image Generated by the 
code. 

Line 

Selecting the Line button draws a wire frame mesh of the 
model within the bounding box. 


X-frame 

This button is not in use at this time. 


Outbox 

Draws the bounding box for the model. 

Solid 

Draws a hidden plot model of the yam. 

Reset 

Resets the viewing parameters to default and erases the 
screen. 
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The Viewing Buttons 

The viewing buttons are used to position the outbox to look 
at the unit cell from different aspects. 

Randr 

The R button is used to rotate the outbox counter 
clockwise. The r button is used to rotate the outbox in the 
same direction but at a slower rate. 


Land 1 

The L button is used to rotate the outbox clockwise. The 1 
button is used to rotate the outbox in the same direction but 
at a slower rate. 


U and u 

The U button is used to move the outbox up in the positive 
z direction. The u button is used to move the outbox in the 
same direction but at a slower rate. 



The Size Buttons 

The size buttons control the size of the outbox 
Larger 

The Larger button increases the screen size of the outbox. 
The effect is to zoom in along the viewing line. 

Smaller 

The Smaller button decreases the screen size of the 
outbox. The effect is to zoom out along the viewing line. 

lands 

The 1 and s buttons are the slow speed counterparts of the 
Larger and Smaller buttons. 
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Menu Hid 

The Menu Hid button hides the menu. This is useful for 
talcin g snapshot of the screen. 


Sw Direct 

This button is not in use 


Backwarc 

Cutting 

Forward 


The Cutting Buttons 

Backward 

The Backward button moves the cutting plane along a 
normal to the plane in the negative direction. 

Forward 

The Forward button moves the cutting plane along a 
normal to the plane in the positive direction. 


Get Rast 


Move Rast 
Load Rast 

Restore 
CutFil 


The File Buttons 

The file buttons manipulates the sun raster files. 

Get Rast 

The Get Rast button reads are raster file and displays it 
on the screen. 


Move Rast 

The Move Rast button creates a raster file with the 
current screen image. 

Load Rast 

The Load Rast button displays the contents of a raster file 
Restore 

The Restore button restores the position of the cutting 
plane 
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( EXIT ) The Exit Button 

Selecting the exit button quits the program and closes the 
model window. 


6.1.6 Starting the Rendering program inside SunView 

The rendering code when started from outside SunView automatically invokes SunView. 
That is the easiest way to start the code. On the other hand starting the code from inside the 
code has the advantage that the user can predetermine the size of the window in which the 
model will be drawn. The user can also have other windows open that he can use in 
conjunction with running the program. In this section limited instructions are given to 
manipulate the windows for such a purpose. For exhaustive discussion on the SunView 
environment the user is referred to [6. 1.2 ]. 

To start the rendering program place the cursor in a command or shell window. Fig. 6 
depicts a typical SunView desktop with a open shell window. The frequently used 
techniques are described below 

Opening a New Shell 

New Shell Windows are opened from the SunView menu. SunView menu is opened by 
clicking the right mouse button when none of the windows are active. From the SunView 
menu selecting the item Shells pops up a new Shell. 

Resizing the Shell Window 

The size of the shell window determines the portion of the screen that will display the 
model. To change the size of the shell window the cursor is placed on the top window 
border. Pressing the right mouse button pops up a Frame Menu. The Resize menu item is 
used to control the window size. The Resize is a pull right menu. The choices offered by the 
menu are constrained and unconstrained resizing, zoom and Full Screen. The cursor is then 
placed near the comer or frame border. Pressing the middle mouse button the bounding box 
can be adjusted to requirement. Zoom expands the window vertically to maximum and Full 
Screen expands the window to encompass the whole screen. 

Exiting SunView 
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To exit SunView the SunView menu is popped up and the item Exit SunView is selected. 


6.2 Application for Engineering Design 

The visualization of the architecture on the computer allows us to look at the model without 
going through the time consuming and expensive operation of actually making a prototype. 

Once a textile pre-form is decided upon it can be readily discretized into a finite structure. 

The software developed can allow us to connect it to a structural analysis program. 

6.2.1 Automesh Generation 

The graphical rendering of any shape requires the object to be discretized into a finite number of 
polygons. The larger the number of polygons is, the realistic the depiction is. This discretization 
can directly be used to obtain a finite element model of the unit cell. The polygons that describe the 
surfaces are taken as a single element and the triangular elements Each triangle is systematically 
broken up into four smaller triangles. Adjacent triangles are then combined to form as many 
quadrilateral elements as possible. 

6.2.2 A Case Study 

As an example of application of the modeling techniques in engineering design we present a work 
carried out at the fibrous materials research center. The center was called upon to design and 
fabricate flanged tube for a helicopter. The flanged tube was designed to transmit torque from 
engine to the hub of a helicopter blade assembly. The tube was to be made of composite material. 
It was suggested by the client to use IM7-W 12k graphite fiber and to use PR500 resin as the 
matrix. The matrix and fiber properties are given in Table 6.2 


Properties 

IM7-W 

PR500 

Tensile Strength(Ksi') 

769 

8.3 

Tensile ModulusfMsi) 

44 

.507 

Tensile Strain(%) 

1.79 

1.9 

Densitv(g/cc) 

1.78 

1.245 

Tg ( 0 O 


205 


126 




The flanged tube was subjected to static bending and fatigue loading. The static loading comprised 
of 771,400 in-lb torque and 38,750 in-lb moment; the fatigue loading was 175,000±175,000 in-lb 
torque and G±27,680 in-lb moment The torque and bending moment occured at the same time. 

The fabrication process included first layer 2-D braiding on tubular flange, thickness built-up with 
woven fabrics and then overbraiding on tube and flange for final process. The braiding angle of 
the flange tube is ±45°. The fiber volume fraction is about 50%. 

Geometric Modeling 


In braiding operations, all braiding yams move with their carriers and wrap around the surface of 
the mandrel. One group of braiding yams moving clockwise forms an angle of +0 with the 
mandrel axis; the other group of braiding yams, moving anticlockwise, has an angle of -0. These 
two groups of braiding yams (±0) are interlocked, forming a biaxial fabric on the mandrel. 

Braids can be formed with different yam interlacing patterns in a manner similar to that of woven 
fabrics simply by changing relative position of carriers on the track ring. If one bias yam 
continuously passes over and under one yam of the opposing group, the pattern is designated as 
1/1 (or diamond) braid. Other simple interlacing patterns in common use include 2/2, 3/3, 2/1 and 
3/1 braids. After examing each pattern and considering the resulting volume fraction and effective 
moduli the 2/2 biaxial braid was chosen for the flanged tube. The effective properties were then 
used in a finite element model to complete the design process. 

Finite Element Modelling 

The finite element analysis(FEA) is performed in order to examine the stress distribution in the 
flanged tube when it is under the operation environment. The FEA procedures include the 
definition of the geometry, the boundary conditions, the material properties and element type in the 
analysis. In this project, 2-D shell element is chosen for the finite element analysis. In the analysis 
the total number of nodes is 260, and the total number of elements is 240. The created geometry 
and the element meshes are shown in Figure 6.6 from different views. 

The element type used for the finite element analysis is nonlinear composite quadrilateral 
shell element This is a 4-node muti-layer quadrilateral shell element with bending 
consideration, which is capable of analyzing three-dimensional structures. Six degrees of freedom 
(three translations and three rotations) per node are considered. 
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For the boundary conditions, fixed dispaceraents at the end of tube are imposed, while the torque 
and bending moment are applied on the edge of the flange according to the provided loading 
environment The loadong conditions are illustrated in Figure 6.7. Both the torque and the 
bending moment can be equated by a set of force components uniformly distributed over the nodes 
on the edge of the flange. For the purpose of simplicity, only 9 inches of the tube in length is 
considered in the analysis. 



Figure 6.7 The loading boundary conditions of the flanged tube. 


The stiffness matrix of the flanged composite tube obtained using the Fabric Geometry Model is 
given as follows: 
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6.15 4.88 0.27 0 0 

4.88 6.15 0.27 0 0 

0.27 0.27 1.33 0 0 

0 0 0 0.63 0 

0 0 0 0 0.63 

- 0 0 0 0 0 


0 -I 
0 
0 
0 
0 

5.13 


Msi 


The stiffness matrix can be used as input into the COSMOS finite element package to compute the 
stress and displacement distributions of the flanged composite tube. The results, such as a x , <Jy, 
a z , T X y, T X z, Tyz, von-Mises stress, Uy,Uz and resultant displacement Utot, are used to see if any 
material or design limit is exceeded. As an example of the von Mises stress fringes are presented in 
Figures 6.8. From the figures, the maximum stress occurs in the tube portion and the value of the 
maximum stress does not exceed the material limit. 
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Uontl i ses 

■ 6438 1 
MB 60586 
I 1 5679 1 

■ 52996 
HW 4920 1 
HI 45406 
OH 41611 
H 378 16 
M 34022 
H 30227 
H 26432 
H 22637 
H 13842 
H 15047 
EH 1 1252 

I I 7456 . 6 


Figure 6.8 The contour plot of von Mises stress in the flanged tube. 
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Appendices 




Model Generators 

3 D BRAID - MOVEMENT GENERATOR 

/* Movement Generator */ 

# include <math.h> 

#inciude <usercore.h> 

^include <stdio.h> 

#define SIZE 50 
#define FASLE 0 

static short x_matrix[SIZE][SIZE][2] > y_matrixiSIZEISIZEI2]4nitial[SIZEl[SIZE]; 

static int n_track,n_coiumn; 

static int movement^ightjUpjtrack, column; 

static int yarnjndex; 

static float dist,yam JocadSIZE][3]; 

int yam_color,n_pair,n_knot; 

short yam_matri?dSIZE][SIZE][10l; 

struct { 

float x; 
float y; 
float z; 

)knotxyz[100]; 

struct { 

short indexO; 
short indexl; 
short i; 
short j; 
short ii; 
short jj; 
short flag; 

}knot[100]; 

mainO 

{ 

int i j^ii; 

FILE *fptr, 

again: printfTNnHow many tracks? M ); 

scanfT%d",&n_track); 

printfTNnHow many columns? M ); 

scanfT%d’\&n_coluinn); 

printfTNnHow many movement? 1 '); 

scanfT%d’\&movement); 

printf( M \nDi stance?"); 

scanfT%f\&dist); 

yam_index=n_knot; 

if (n_track>SIZE I l n_column>SIZE ) 

printfTNnLoom too big! # of track or # of column should be less than 50. \n Try again!”); 
goto again; 

} 


1 


buil d_xy _m a t rix( ); 
print jcyjnatrixO; 


yam_index*0; 
for (i*0d<n_track4++) 

for(j=0;j<n_column;j++) 

if (y_matrix[i][j][0] !*0 I I x_matrix[i][j][0] !=0) 

{ 

yamjndex++; 

initial{ilj]=l; 

} 

for (i*0;i<n_track;i++) 

{ 

for(j=03<n_column;j++) 

printff"%d H > initial[i][j]); 

printfTW); 

} 

fptr » fopenCmovement.dat’VV’); 

/*fprintfTfptr,"%d %d\n",yarn_i ndex, movement/2+1); 

*/ 

fprintflfptr,"%d %d\n",yamJndex,movement+l); 

movement++; 
yam_color = 3; 
yam_index=0; 
for (i=Od<n_trackd++) 

for (j=0*j<n_column*j++) 

{ 

if (initial[i][j] !=OY* if initial[i][j]=l, it is yarn! */ 

{ 

for ( k=0;k<movement;k++) 

{ 

if (k=0) 

{ 

track=yarnJocat[k][2]=i+l; 

yam_locat[k][l]=0;/* y coord, is ZERO */ 

column=yam_locatfk]tO]=j+l; 

right=0; 

up=0; 


else 


if (k%2 !=oy* move track first */ 
move_track(k); 

else 

move_column(k); 

} 

yamJocatfk][l]= k*dist ; 

J /* end of movement */ 

printfTyam # ^dXn^yamJndex); 

fprintftfptr/ , %d\n",yarn_color); 
yam_color= (yam_color==3) ? 4 : 3; 
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for (k*0 , je<movement^++) 

/* if(k%2==0)*/ 

{ 

printf("%2.4f %2.4f %2.4f\n",yam_locattk][0],yamjocat(kll]oram Jocatlk][2]); 
fprintf(fptr,"»2.1f %2.1f 

%2.1f\n" jam locat(k][0],yam_locat[k][l],yamjocat(k][2]*0.7); 

} 

fprintflfptr,"\n"); 

yam_index++; 

} 


} /* end of yam location */ 
initially am_matrix(); 
iaO; 

for (k»0;k<movement-l;k++){ 

if(k%2) shift_track(i,i+l); 
else shift_column(i,i+l); 
if(k%2 && k !=0 M 
n _pair * 0; 
knot_search(i-l,i+l); 
pickJcnot(n_pair,k*dist); 


i++; 

} 

for (i=K)p<n_knotd++) 

printfT\n%2.2f %2.2f %2.2f ^knot3cyzfi]jc,knotxyzQ]^^ 

fprintflfptr, M %d\n , \n.knot); 
for (i=0d<n Jcnot;i++) 

fprintflfrtr,"\n%2.2f %2.2f %2.2f ”, 

knotxyz[i]jc,knotxyz[i].y,knotxyz[i].z); 
close (fptr); 

} /* end of main program *1 

move_track(k) int k; 

{ 

if (right==0) 

{/* move right */ 

column*yarnjocatfk][0]=column+x_matrix[track-ll[column-l][right]; 

yam _locat(k][2]=yam Jocat[k-l][2]*/* keep Z direction */ 


rights 1; 

} 

else/ 1 *' move left */ 

{column=yam_locat[k][0]=column+x_matrix[track-l][column-l][right]; 
yam_locat[k][2]=yam Jocat[k-l][2]*/* keep Z direction */ 

right=0; 

1 


} f* end of move_track */ 


move_column(k) int k; 

{ 

if (up=FALSE) 
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else 


{ 

track-yarn Jocat0c][2]»track+y_matrix[track-Hcolumn-lIup]; 
yam _locat[k][0]=yam Jocat[k-l ][0]; 


l 

track-yam Jocat(k][2]-track+y_matrix{track-lIcolumn-l][up]; 
yam Jocat[k][0]=yamJocat[k-l][0]; 


up-O; 

} 

} /* end of move_column *1 


print_xy_matrix() 

{ 

intijjc; 

printfTXn"); 
for (k=0;k<2;k++) 

for(i=0*4<n_trackp++) 

{ 

if(k=0) 

{ 

for(j=0;j<n_coluinn;j++) 

printfT%2d "pc_matrix[i][j][0]); 
printfT\n M ); 

} 

else 

{ 

if (k=l && i=-0) printfTW); 
for(j=0;j<n_column;j++) 

printfT%2d ”pc_matrixli]lj][l]); 
printfCW); 

) 

} 

printfT\n M ); 
for (k-0;k<2;k++) 

ford- 04 <n_track;i++) 

( 

if (k==0) 

1 

for(j=0;j<n_column;]++) 

printfT%2d ”,y_matrix[i][j][0]); 
printfTNn"); 

1 

else 

{ 

if (k=l && i==0) printfTW); 
for(j=Oy<n_coluinn*j++) 

printfT%2d ”,y_matrix(i][j][l]); 
printfTW); 

) 

} 

} /* end of print x_y-matrix */ 
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build^jcy.matrixO 

{ 

intijjt; 

for (i»0;i<n_track;i++) 

for(j=OJ<n_column;j++) 

for(k*0;k<2;k++) 

x_matrix[i][j][k]=y_matrix[i][jMairutialti]Ij]»0; 


for(k=0;k<2;k++) 

for(i=0;i<n_track;i++) 

{ 

if (k=0) 

{ 

ifli !=0 && i !=n_track-l && i % 2 !=0 ) 
for(j=0;j<n_column-l;j++) 

x_matnx(iljl[0]= 1; 
if (i !=0 && i !=n_track-l && i % 2 !=1 ) 
for (j=ly<n_column;j++) 

x_matrix[i][j][0]= -1; 


else 


ifft !=0 && i !=n_track-l && i % 2 !=0 ) 
for(j = l*j <n_colum n 3 ++ ) 

x_matrix[i][j]Il]= -1; 
if (i !=0 && i !=n_track-l && i % 2 !=1 ) 
for (j=0;j<n_coluinn-l;j++) 
x_matnx(i]]j][l]= 1; 


) 

} f* end of build x.matrix */ 


for(k=0;k<2;k++) 

for(i=Od<n_column;i++) 

{ 

if (k=0) 

{ 

iffi !=0 && i !=n_column-l && i % 2 !=0 ) 
for(j= 0 *j<n_track-l*j-H-) 

y_matrixfj][i][0]= 1; 

if (i !=0 && i !=n_coiumn-l && i% 2 1 = 1 ) 
for (j=l;j<n_track;j++) 

y_matrixij][i][0]= -1; 


else 

{ 

ifli !=0 && i !=n_column-l && i % 2 !=0 ) 
for(j = 1 *j <n_tra ck*j ++) 

y_matrix(j][i][l]= -1; 
if (i !=0 && i !=n_column-l && i % 2 !=1 ) 
for (j=0*j<n_track-l;j++) 

y_matrix(j][il{l]= 1; 


) 

} /* end of build y_matrix */ 
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/* end of build xy_matrix 


1 


*/ 


initially am_matrix() 

{int ij^yamjndex; 
yamjndex » 1; 
for (i=*0^<n_trackp-H-) 

for (j=Oy<n_columny++) 

yam.matrudi][j][0] = (initialB][jl«l) ? yamjndex++ :0; 
} /* end of initial_yam_matrix *1 

shift_track( indexO, index Dint indexO, index 1; 

{int i j; 

for (i=0;i<n_track;i++) 

for (j^O;j<n_column;j++) yarn_matrix[i][j][indexl] * 
yarnjnatrixli][j][index0] ; 


for (i» 1 d <n_tr ack- 1 d ++)1 

if (yarn_matrix[i][0][indexl]>0){/* RIGHT */ 
for (j=l-j<n_column*j++) 

yarnjnatrix[i][n_column-j][indexl] = 

y arn_m atrix[i][n_column-j- 1 ][index 1 ]; 
yam_matrixli][0][indexl] ■ 0; 


else {/* LEFT V 

for (j=0J<n_column-l*j++) 

yam_matrix{iltjlindexl] = 

y arn_matrixti][j+ l][indexl]; 
yarn_matrix[i][n_column-l][indexl]= 0; 

) 


}/* the next track */ 

} /* end of shift track */ 

shift_column( indexO, indexDint indexO, indexl; 
lint ij; 

for (i=0p<n_track*4-H-) 

for (j= 0 *j<n_column 3 ++) yam_matrix{i][j][indexl] = 

yam_matrix[i](j][indexO] ; 


for (i=l;i<n_column-l;i++){ 

if (yam_matrix(0][i][indexl]>0){/* DOWN */ 
for (j=l J<n_track;j++) 

yam_matrix[n_track-j][i][indexl] = 

yarnjnatrix[n_track-j-l][i][indexl]; 
yarn_matrix{0][i][indexl] = 0; 

1 

else {/* UP */ 

for (j=0*j<n_track-ly++) 

yam_matrix[j][i][indexl] = 

yam_matrixlj+l][ilindexl]; 
yarn_matrix[nj;rack-l][i][indexl]= 0; 

} 

1/* the next column */ 

) /* end of shift column */ 
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knot_search(indexO,indexl)int indexO,indexl; 


{ 

intiJ4ijj; 

for (i»Od<n_trackp++) 

for (jaO*j<n_columny++){ 

if (yam_matrix[illj][indexl] !=0 ){ 
for (ii*0;ii<n_track;ii++) 

for (jj=0*jj<n_columni++){ 

if (yarn^matrixfilljllindexllMyarn^matrixliilijlindexO] && 

yarn_matrix[ii]|jj]Iindexl]**yarn_matrix[i]|j][indoxO] && 
yam_matrix[ii][ij][indexl] !» 0){ 

/*printflr\nPAIR %d 

%d^yarn_matrix[i][j][indexl], 

yarn_matrix[ii][ju][indexl]);*/ 

knot[n_pair].indexO = yarnjnatrix[i][j][indexl]; 

knotfn_pair].indexl = yam_matrix(ii]Ijj][indexl]; 

knot[n_pair].i =i; 

knotfn_pair].j = j; 

knot[n_pair].ii = ii; 

knot(n_p air]. flag = -1; 

knot(n_pair++].jj =jj; 

} 

) 


} 

}/* end of knot search */ 


pidc_knot(n_pair,y)int n_pair;float y; 

{int ij^k; 

for (i=0;i<n_pair;i++){ 

if (knot[i].flag == -1){ 

knot[i].flag= 1; 
for(j=0’j<n_pair*j++) 

if ((knot(i].indexO =knot{j].indexO && knot(j].flag != 1 && 
knot(i].mdexl=knot{j].indexl) 1 1 
(knot[i].indexO = knot(j].indexl && knot(j].flag != 1 && 
knot{i]. index 1 = knot(j]indexO)) knot(j]ilag = 0; 

} 


for (i*0;i<n_paird++) 

if (knot[i].flag = 1){ 

knotxyz[n_knot].z = (knot(i]i + knotji] ji+2)/2.0; 
knotxyz[nJcnot]jc = (knot{i].j + knodi].ij+2)/2.0; 
knotxyz[n_knot++].y = y; 

1 

} 


MWK - KEYPOINT GENERATOR 


/* 


mwk unit cell builder 
history: 

original development in Fortran -> Dr. Charles Lei 
C conversion Anisur Rahman 4^16/93 


*/ 
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#indude <stdio Ji> 

#indude <atdlib Ji> 

#include <mathJi> 

idefine LOOP 10 

#define POINTS J>ER_STITCH 33 
fdefine DIMENSIONS 3 

int ncolor; 

int indl, ind2, ind3; 

float din; 

double pi, dx,dy,sqrt2, theta; 
double rstitch, rinsert, rloop; 

double coor[LOOP][POINTS_PERJSTITCH][DIMENSIONS]; 

FILE *fp; 

main(){ 

puts (" Enter # of insertion yam/unit length"); 

scanfT%f,&din); 

dx=l-/din; 

puts (" Enter the angle (degree) of bias yams "); 

scanfT%lg’\&theta); 

pi *4.* atan(l.); 

theta»theta/180.*pi; 

dy = dx * tan(pi/2. - theta); 

puts( M Enter the radius of the insertion yam”); 

scanfC%lg M ,&rinsert); 

putaC Enter the radius of the stitch yam”); 

scanf( H %lg H ,&rsti tch); 

putsC Enter the radius of the stitch loop”); 

scanfT%lg",&rloop); 

/* 

Compute the geometry within an unit cell 
*/ 

for (indl = 0;indl != LOOP; indl++) 
for (ind2 =0; ind2 != POINTS_PER_STITCH; ind2++) 
for (ind3=0; ind3 != DIMENSIONS; ind3++) 
coor[indl][ind2][ind3]=0; 
coor[0][0][0]=dx/2.-rstitch; 
coor{0][0][l]=-rstitch - rloop; 

coor[0][l][0]=dx/2.-4.*rstitch*(2.*rstitch/dy)-rstitch; 

coor[0][l][2]=-2.^stitch-(2^stitch/dy)*4.*rstitch; 

coor{0][2][0]=dx/2.-rloop; 

coor[0][2][l]=dy-rloop*rstitch; 

sqrt2=sqrt(2.); 

coori0][3][0]=dx/2-rloop/sqrt2; 
coor{0][3][l]=coor{0][2][l]+rloop/sqrt2; 
coor[0][4][0]=dx/2; 
coor[0][4][l]=dy-rstitch; 
coor[0][5][0]=dx/2.+rloop^sqrt2; 
coor{0][5][ l]=coor[ 0][3][1] ; 
coor[0][6][0]=dx/2. +rloop; 
coor(0][6][l]=coor[0][4][l]-rloop; 
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cooK0l7][0]»clx/2.+r»titch+2,*r8titch*(2.*r8titch/dy); 

coorf0l7Il]*-ratitch; 

coorf0l7I2]*-2 *rstitch; 

coorf0l8l0]*coor[0][4l0]+r8titch; 

coorf0l8Il]*-rloop-rstitch; 

coorf0][9l0]*dx/2.; 

coor(0][9Il]ady/2.-rinsert*8qrt2; 

coorfOI9][2]=7.*rinsert; 

coor(0][10][0]«dx/2.; 

coor[0ll0][l]=dy/2.; 

coorf0][10I2]=8.*rinsert+r8titch; 

coorfO][ 1 l][0]=dx/2, ; 

coor[0][ 1 1 ][ l]ady/2.+sqrt2*rins ert ; 

coorfOl 1 1I2>» 7.*rinsert; 

coor{0]I12][0]*coor{0][8][0]; 

coor[0][12][l]=coor[0]l6][l]; 

coor{0][13][0]=coort0][7][0]; 

cooK0][13][l]=coor[0][7][l]+dy; 

coor(0][13][2]=-2.*rstitch; 

coorf0][14l0]=coor{0][6][0]; 

coor(0][ 14][ l]=coor[0][6][ l]+dy; 

coor{0][14][2]=coor[0][6][2]; 

coor[0][15][0]=coor[0][5][0]; 

coor[0][15][l]=coor[0][5][l]+dy; 

coor{01[16l0]=coor(0]r4][0]; 

coor[0][16][l]=coor[0][4][l]+dy; 

coor[0][17][0]=coor(0][3][0]; 

coor[0][17][l]=coor(0][3][l]+dy; 

coor{0][18][0]=coor(0][2][0]; 

coor[0][18][l]=:coor(0][2][l]+dy; 

coor[0][19][0]=dx/2.-rstitch-2. # rstitch*(2.*rstitch/dy); 

coor{0][19][l]=dy-rstitch; 

coor[0][19][2]=-2.*rstitch; 

coor[0][20][0]*dx/2.-rstitch; 

coor[0][20][ l]=dy-rstitch-rloop; 

coor(0][21][0]==coor(0][9][0]; 

coor[0][21][l]=coor[0][9][l]+dy; 

coor[0][21]{2]=coor{0][9][2]; 

coor[0][22][0]=coor[0][ 10][ 0]; 

coor[0][22][ l]=coor[0][ 10][l]+dy; 

coor[0][22][2]=coor{0][10][2]; 

coorf 0][23][0]=coor[0][ 1 1 ][ 0] ; 

coor[0][23][ l]=coor[0][ 1 1 ][ 1 ]+dy ; 

coor[0][23][2]=coor[0][ll][2]; 

coor(0][24][0]=coor(0][20][0]; 

coor(0][24][l]=coort0][20][l]+dy; 

coor[0][24][2]=coori0][20][2]; 

coor[0][25][0]=coor{0][ 19][0]; 

coor{0][25][l]=coor[0][19][l]+dy; 

coor[0][25][2]=coor{0][19][2]; 

coor{0][26][0]=coor[0][l][0]; 

coor[0][26][l]=coor[0][l][l]+dy*2; 

coor[0][26][2]=coor[0][l][2]; 

/* 

Write Sun File 

*/ 

fp = fopen("outfilY w ); 

f^rintfttp,"! 27\n"); 
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putaT color of the stitch yam: 1 , 2 , 3, 4 \n"); 
8canfi’%d",&ncolor); 
fprintfffp > "%d\n”,ncolor); 
foKindl - 0; indl !« 27; indl-H-) 


fprintf(fp,"%10.6f %10.6f %10.6f\n", coortOlfindllO], 
coor[01ind l][l],coor(OI[ind 112]); 


3D BRAID YARN BUILDER 

#include <math.h> 
findude <usercore.h> 
finclude <stdio.h> 

#define NC 50 
#defme NSET 6 
fdefine MAXLINE 400 
#define MAXVLINE 5000 
#define PI 3.141596 

float cp(NC][3]; 
float cm£NC*2][3]; 

float bbox(3][2]; 
intpi,pj^LN; 

double a{3][3]; 

float radial l,radial2,u,v,x,y,z; 
float pp[60000][3]; 
int index(60000][6],ind; 

int count ,nvert,npoly,nzset,nvpzset,countl,total_vert,old_count; 
int polyJndex,finish__yarn, movement, yarn_color; 
int knotkjmotn^ic; 
float linejvert(MAXVLINE][3]; 
int Une_index[MAXLINE][8],nline,nlvert; 
static float cube{8][3]=:{{0.0, 0.0, 0.01, 

{ 0 . 0 , 0 . 0 , 1 . 01 , 

{ 1 . 0 , 0 . 0 , 1 . 01 , 

{ 1 . 0 , 0 . 0 , 0 . 01 , 

{0.0,1.5,0.01, 

{0.0,1.5,1.01, 

{1.0,1.5,1.01, 

{1.0,1.5,0.011; 

static float pmatrix[4l4]={{-0.5,1.5,-1.5,0.5(, 

{1.0, -2.5,2.0, -0.51, 
{-0.5,0.0,0.5,0.01, 
{ 0 . 0 , 1 . 0 , 0 . 0 , 0 . 011 ; 


FILE *fptr; 
struct v_p_count{ 

short v_count; 
short p_count; 

) v_p_count[600]; 

main(argc, argv) 
int argc; 
char *argv(]; 

{ 
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intij,ii,ncube; 


int kAtmp[7],color, 
int n_poly_pobj,n_vert_pobj 
printfl"\nYam radial 1? "); 
scan«*%r,&radiall); 
printff’XnYam radial 2? "); 
scanf("%f , &radial2); 


,vert_count,poly_count; 


poly.index ■> total_vert = old_count = 0; 

finish_yam = 0; 

fptr = fopen(argv[l],"r H ); 

fscanKfptr ,"%d",&movement); 

fscanfffptr,"%d",&nc); 


for (finish_yam =0;finish_yarn<movement;finish_yam++) 

( 

fscanf(fptr,“%d",&yam_color); 

for (i=0; i<ncp++) 

fscanflfptr, "%f %f %f',&cp[i][0],&cp[i][l],&cp(i][2]); 


ind=0; 

/* printfi"#%d",finish_yam); */ 

/*build_control_mesh(nc);*/ 

v_p_count[fmish_yam].v_count=v_p_count(finish_yam].p_count 

build_yarn(); 


1 

/* 

cellO; 

*/ 

close (fptr); 

fptr = fopen (”n.dat",'V); 


for (i=Op<total_vert;i++) 


if (i == 0) 

{ 

bbox{0][0]=bboxI0][l]=pp[i][0]; 

bbox[l][0]=bbox[l][l]=pp{i][l]; 

bbox£21[0]=bbox£2][ll=ppti][2]; 

1 

else 

{ 

if Cbbox£0][0]<pp[i][0]) bboxt0][0] = pp[i][0]; 
if (bbox£01[ l]>ppti][0]) bbox(0][l] = pp[i][0]; 
if (bbox{l][0]<pp[i][l]) bbox{l][0] = pp[i][ll; 
if Cbboxf 1][ ll>ppfilt 1]) bboxll][l] = pptiltl]; 
if (bboxt 2][0]<pp[i][2]) bboxt2][0] = pp[i][2]; 
if (bbox{21[l]>pp[i][2]) bbox[2][l] = pp[i][2]; 
} 

1 
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%4.4f 


% 4 . 4 f 


fprintfl[fptr,"%d %d %d\n", movement, poly jndex,total_vert); 
fprintf(fptr,"%4.4f %4.4f %4.4f 

%4.4f\n H ,bboxt0][0],bbox[0ll],bboxa][0],bboxn][l],bboxl2][0],bbo3d2Il]); 

poly_count = vert_count ■ 0; 

n_poly_pobj * poly_index/movement; 

n_vert_pobj » total_vort/movement; 

for (i=*0?< movementd++){ 

n_poly_pobj=v_p_count(i].p_count; 
n_vert_pobj=v_p_count(i].v_count; 
fprintf(fptr,"%d %d \n’ - ,n_poly_pobj,n_vert_pobj); 
for ( j=0a<n_vert_pobjy++){ 
fprintRfptr,"%f %f %f\n",pp[vert_count][0], 

pp[vert_count][l],pp[vert_count][2]); 

vert_count++; 

} 


for ( j=03<n_poly_pobj3++){ 

fprintfi:fptr,’'%d %d %d %d %d %d\n",index[poly_count][0], 
index[poly_count][ 1 ],index[poly_count][2],indextpoly_count][3] 
,index[poly_count][4],index[poly_count][5]); 
poly_count++; 

} 

} 

i =j = 0; 

fprintflfptr,"%d %d \n",i j); 
close (fptr ); 

} /* end of main */ 


cell() 

{ 

int ij,k4i,ncube, color; 
int tmp{7]; 
float x,y,z; 


fscanfTfptr,''%d%d",&ncube,&color); 
if (8*ncube>MAXVLINE) 


printfT"\n Too many vertices on lines!!"); 
exit(l); 

} 

nlvert=0; 

for (i=0p<ncubed++) 

{ 

fscanflfptr,"%f %f %f*,&x,&y ( &z); 
for (j=0-j<8-j++) 

1 

Iine_vert(nlvert][0]=cube|j][0]+x; 

line_vertfnlvert][l]=cube{j][l]+y; 

line_vertfnlvert][2]=cube[j][2]+z; 
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nlvert++; 

} 

} 

tmp(0]=5;tmp{ l]=color, 
for (i=0;i<ncube;i++) 

( 

k=i*8+5; 
for (j=la<43++) 

{ 

tmp[2]=j+i*8; 

tmp(3]=k; 

tmp[4]=k+l; 

tmp[5]=tmp[2]+l; 

tmp[6]=tmp[2]; 

k++; 

for (ii=0;ii<7;ii+-t-) 

line_indextnline][ii]=tmp[ii]; 

nline++; 

} 

tmp[2]=i*8+4; 

tmp[3]=tmp[2]+4; 

tmp[4]=tmp[3]-3; 

tmp[6]=tmp[2]-3; 

tmp[6]=tmp[2]; 

for(ii=Opi<7;ii-H-) 

line_indextnline][ii]=tmp(ii]; 

nline-M-; 

if (nline>MAXLINE) 

{ 

printfTXn Too many lines!!"); 
exi«2); 

} 

1 

} t* end of cell */ 

build_yarn() 

( 

int kJvj4i; 
int tmp[6]; 
int p_count; 
float ppp[10000][3]; 

float tangent[3],vect01[3],vectl2[3],b_norm[3],vectpp[3],vect_temp[3],bt(3], 

vect_temp0[3],vect_templ[3]; 

float dot_product(); 

float bn[3]; 

float bend, alpha, p[3],anglel,f,pv[3], twist; 
int flag,num,go_back; 

float vtt[3],w[3],vp[3],vect00[3],vt[3],twist_angle,new_bt[3],new_bn[3]; 
float old_b_norm[3],old_bt[3],twist_angle_i; 

struct path 

{ 

float bt[3]; 
int node; 
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int twist _flag; 

} path! 200]; 
k«3 ; 

n*ind-l; 
p_count»0; 
flag - TRUE; 

printfT\nOLD POINTS\n M ); 
for (i=K)?<=anp++) 

printfT%4.4f %4.4f %4.4f\n\cm[i]l0],cni[i][l],cm[i][2]); 

for (j=0g<=nc-43++){ 

for (u=0.0;u<=.9001;u=u+0.2) 

{ 

cards p(uj); 

PPPfP- coun t][0]=x;ppp[p_count][l]j=y^)pp[p_coiint][2] 

p_count++; 

} 

} 

printfTXnNEW POINTSW *); 


/* 


*/ 


for (i=0;i<p_count;i++) 

{ 

path{i].twist_flag=FALSE; 
if (fabs(cm[j][l]-ppp[ill])<=0.05) 
{ 

path! i] .node =TRUE; 

j++d++; 

) 

else 

path[i].node=FALSE; 

} 


for (i=0;i<p_countp-M-) 

printfT%4.4f %4.4f %4.4f\n' , ,ppp[i][0],ppp[i][l] f ppp[i][2]); 


count =0; 

anglel=0.0; 
flag = 1; 

for (i=0^<p_count;i++) 

{ 

printfl" %d\n l *,i); 

if (i=0 ) f* no bending and twisting at the first points */ 


bend=0.0; 
twist=0.0; 
twist_angle * 0.0; 
for (j=0*j<3*j++) 


vectOO[j]=ppp[i+l][j]-ppp[i]Ij]; 

vect01Ij]«ppp(i+2][j]-pppfi][j]; 

} 

cross_product(&bt[0],&vect00[0],<Stvect01[0]); 

u=dot_product(&vect00[0],&vect01[0]); 
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printf("%2.4f \n\u); 

crosi^producU&b^normfOl^&btfOl.&vectOOtO]); 
normalize(&bt[0]); 
normalize &b_norm[0]); 
for (j«03<NSET3++) 

( 

alpha* twist_angle + 6.283*yNSET ; 

pp(total_vert][0]=.ppp[0][0]+radiall*cos(alpha)*bt(0]+radial2*8in(alpha)*b_normt0]; 

ppt total _vert][ l]=ppp[0][l]+radial l*cos(alpha)*bt( l]+radial2*sin(alpha)*b_norml 1]; 

pp[total_vert][2]=ppp[0][2]+radiall*cos(alpha)*bt(2]+radial2*sin(alpha)*b_nonnC2]; 

total_vert++; 

count++; 

} 

1 /*endofi=0 */ 

else 

if ( i = p_count-l) /* the last one as same as i-1 point */ 

{ 

for (j=0a<3y++) 

I 

vect00[j]=pppli-l][j]-ppp{i-2][j]; 

vect01[j]=ppp[i][j]-ppp[i*l][j]; 

1 

/*cross_product(&bt[0],&vect00[0],&vect01[0]); 
u=dot_product(& vect00[0], &vectO 1[ 0]) ; 
printfl"%2.4f \n“,u); 

cros8_product(&b_norm[0],&btl0],&vect01[0]); 

normal! ze( &bt(0]); 

*/ 

bt[0] = old_bt[0];bt[ 1] = old_btUl;bt[2] = old_bt(2]; 

cross_product(&b_norm[0] ( &btl0],&vect01[0]); 

normalize(&b_norm[0]); 

for (j=0g<NSET*j++) 

{ 

alpha= twist_angle + 6.283*j/NSET ; 

pp(total_vert][0]=ppp[i][0]+radiall*cos(alpha)*bt(0]+radial2*sin(alpha)*b_nonn[0]; 

pp[totaLvert][l]=ppp[i][l]+radiall*cos(alpha)*bt(l]+radial2*sin(alpha)*b_norm[l]; 

ppttotal_vert][2]=ppp[i][2]+radiall*cos(alpha)*bt[2]+radial2*sin(alpha)*b_norm[2]; 

total_vert++; 

count++; 

1 

} /* end of i=p_count-l */ 

else 

if ( i==l) 

1 

for (js0-j<3a++) 

{ 

vectOO[j]=ppp[i][j]-ppp[i- l][j]; 
vect01[j]=ppp[i+l][j]-ppp[i][j]; 
tangent(j]=vect01[j]+vect00[j]; 

1 

cross_product(&bt[0] ( &vect00[0],&vect01[0]); 
u=dot_product(&vect00[0],&vect01[0]); 
printfT%2.4f \n”,u); 

cross_product(&b_norm[0],&bt[0] ( &tangenU0]); 
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normalize(&bt[0]); 
normalize(&b_nonnlO]); 
for (j» 03 <NSET 3 ++) 

{ 

alphas twi strangle + 6.283*yNSET ; 

pp( total jrert][0]sppp[i][0]+radiall*co8(alpha)*bt[0]+radia]2*sin(alpha)*b_norm[0]; 
pp[total_vert][l]sppp[i][l]+radiall*cos(alpha)*bt(l]+radial2*sin(alpha)*b_nonn[l]; 
pp(total_vert][2]sppp[i][2]+radiall*cos(alpha)*bt(2]+radial2*8in(alpha)*b_norm[2]; 

total _vert-H-; 
count++; 

} 

} /* end ofi=l case *1 

else 

if ( i»»p country* start p-2 case */ 

1 

for (j=0;j<3;j++) 

{ 

vectOO(j]=pppti][j]-ppp{i- l][j]; 
vect01{j]=ppp[i+l][j]-ppp[i][j]; 
tangent(j] = vect 0 1 [j] + ve ctO Olj ] ; 

1 

/*cross_product(&bt(0],&vect00[0],&vect01[0]);*/ 
btfOl = old_bt(0];bt{l] = old_btUl;bt{2] = old_btf2]; 

Usdot_product(&vect00[0],&vect01[0]); 
printfT%2.4f \n'\u); 

cross_product(&b_norm[0],&bt[0],&tangent{0]); 

normalize(&bt{0]); 

normalize(&b_norm[0]); 

for (jsO-j<NSET*o++) 

{ 

alphas twist_angle + 6.283*j/NSET ; 

pp[totaljvert][0]=ppp[i][0]+radiall*cos(alpha)*bt[0]+radial2*sin(alpha)*b_norm[0]; 

pp(totaljvert][l]sppp[i][l]+radiall*cos(alpha)*bt[l]+radial2*sin(alpha)*b_norm[l]; 

pp[total.vert][2]=ppp[i][2]+radiall*cos(alpha)*bt{2]+radial2 ,,t sin(alpha)*b_norm[2]; 

total_vert++; 

count++; 

1 

} f* end of i=p-2 case */ 

else 

{/* start the other case */ 
for (j=0;j<3;j-»-+) 

{ 

vect00[j]=ppp[i-l][j]-ppp[i-2][j]; 

vect01[j]sppp[i][j]-ppp[i-l][j]; 

vectl2(j]sppp[i+l][j]-ppp[i]tj]; 

tangent{j]=vect01[j]+vectl2Ij]; 

) 

cross jDroductf&vect.tempOtO^&vectOltO^&vectOOtO]); 
crosa_product(&vect_templ[0],&vectl2[0],&vect01[0]); 
twist= dotj?roduct(&vect_temp0[0],&vectJ;empl[0]); 
if (twist >s 0.99 ) twist = 1.0; 
if (twist <= -0.99 ) twist = -1.0; 
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if ( twist—— 1.0 1 1 twist — -1.0) /* no twisting */ 

{ 

bt[0] ■ old_btlO];bttl] - old_bt(l]Jbtt2] = old_btl2]; 
cross_product(&b _norm[0],&bt[0],&tangent(0]); 
normalize! &b_norm(0]); 
for (j«03<NSETa++) 

( 

alphas twist_angle + 6.283*j/NSET ; 

pp(total_vert][0]»ppp[i][0]+radiall*cos(alpha)*bt[0]+radial2*sin(alpha)*b_norm{0]; 
pp[ total vert][l]=ppp[i][l]+radiall*cos(alpha)*bt(l]+radial2*sin(alpha)*b_norm[l]; 
ppttotalIvert][2]=ppp[i][2]+radiall*cos(alpha)*bt{2]+radial2*sin(alpha)*b_norm[2l; 

total_vert++; 

count++; 

1 

)/* end of no twisting */ 

else 

{ /* handle the twisting first */ 

twist = acos( twist); 

/*if (twist > (PI/2.0) && twist <= PI ) twist -= PI/2.0; 
if (twist > PI && twist <= ( PI *3 /2.0)) twist -= PI; 
if (twist > (PI *3/2.0) && twist <= ( PI *2.0)) twist = 2*PI - twist; 
twists 0.0;*/ 

/*twist_angie» twist_angle - twist;*/ 

printf (" twist_angle %f %f \n",twist_angle, twist); 

cross_product(&bt(0],&vect01[0],&vectl2[0]); 

u=dot_product(&vect00[0],&vect01[0]); 

normalize/ &bt[0]); 

if (fabs( twist) > (PI/6.0)){ 
bt[0] = -btfO]; 
bt[l] = -bt[l]; 
btI2] s -bt[21; 

1 

printfl"%2.4f\n'’,u); 

bt[0] = old_btt0];btll] = old_bt(l];btl2] = old_btl2]; 

cross_product(&b_nonn[0] > &bt(0],&tangent(0]); 

normalize(&b_norm[0]); 

for (j=0 - j<NSET - j++) 

{ 

alphas twist_angle + 6.283*j/NSET ; 

pp[total_vert][0]=ppp[i][0]+radiall*cos(alpha)*btl0]+radial2*sin(alpha)*b_norm[0j; 

pp(total_vert][l]=ppp[i][l]+radiall*cos(alpha)*bt{l]+radial2*sin(alpha)*b_norm[l]; 

pp[total_vert][2]=ppp[i][2]+radiall*cos(alpha)*bt(2]+radial2*8in(alpha)*b_nonn[2]; 

total_vert++; 

count++; 

1 

)/* the end of twisting */ 

}/* end of ( i !=1 && i != p_count-l ) */ 

for (j=03<3'j++)( 

old_bttj]=b0j]; 
old_b_norm[j]= b_norm[j]; 
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} 


/•end of i loop */ 


nvert * count; 

v_p_countffinish_yam].v_count » count; 

printfTXn # of V %d M , count); 

npoly * count-NSET; 

nzsetacount/NSET; 

nvpzsetaNSET; 

tmp[0]*4;tmp[l]=yarn_color; 

/*************** bottum ****************/ 

k *nvpzset+ old_count; 

tmp[0]»4; 

for (i=lp<nvpzset/2*4++) 

{ 

tmp{2]=i+ old.count; 
tmpt3]=tmp[2]+l; 
tmp[4]=k-l; 
tmp[5>=k; 

k-; 

for (ii=0;ii<6;ii++) 

index[polyJndex][ii] = tmpfii]; 
polyJndex++; 

v_p_count[finish_yarn].p_count++; 


y*************** gj^g ****************/ 

for (i=l;i<nzset;i++) 

( 

k=i*nvpzset+l+old_count; 
for (j=l'j<nvpzset*j++) 

( 

tmp[2]=j+(i-l)*nvpzset+old_count; 

tmp(3]=k; 

tmp(4]=k+l; 

tmp[5]=tmp[2]+l; 

k++; 

for (ii=0;ii<6;ii-H*) 

index(polyJndex][ii]=tmp[ii]; 

poly_index++; 

v_p_count[fmish_yarn].p_count++; 

) 

tmp[2]=( i - 1 )*nvp zset+j +old_count; 
tmp[3]=nvpzset+i*nvpzset+old_count; 
tmp[4]=i*nvpzset+l+old_count; 
tmp[5]=l+<i-l)*nvpzset+old_count; 

for (ii=0;ii<6;ii++) 

index(polyJndex][ii]=tmp[ii]; 

poly_index++; 

v_p_count[ finish __yam].p_count++; 

) 
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****** ******** top **** 

k » nvpzset*nzset+old_count; 
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for (i=l;i<nvpzset/2;i++) 

{ 

tmp(2] = i+nvpzset*(nzset-l) + old_count; 
tmp[3] = k; 
tmp[4] = k-1; 
tmpf5] » tmp[2]+l; 

k-; 

for (ii=0;ii<6;ii++) 

index[polyJndex][ii] = tmp[ii]; 
polyJndex++; 

v_p_count[finish_yam].p_count++; 

} 

old_count + = count; 


} I* end of build yarn */ 


cardsp (u,n) 

float u; 
intn; 


{ 

float t[4],pc[4]; 
t[0]=u*u*u; 
tt2]=u; 

pcf0l.tio]*pmatrix[0][0]+t[l]*pmatrixll][0]+t[2]*pmatrixt2][0]+pmatrix[3l0]; 

pcf l]=tt0]*pmatrix[0][ l]+t[ l]*pmatrixll][ l]+t[2]*pmatrix{2][ l]+pmatrix[3][l]; 

pct2]=tl0]*pmatraloj[2]+tll]*pmatrix[l][2]+tf2]*pmatm[2]12]+pmatru43][2]; 

pc{3]=t(0]*pmatrix£0][3]+tt l]*pmatrixt l][3]+t[2]*pmatrix[2][3]+pmatrix[3I3]; 

x = cp[n][0]*pc[0]+cp[n+l][0]*pc{l]+cp[n+2][0]*pc(2]+cp(n+3][0]*pc{3]; 
y = cp[n][l]*pc(0]+cp[n+l][l]*pc{l]+cp[n+2][l]*pc[2]+cpln+3][l]*pcl3j; 
z = cp[n][21*pc{0]+cp[n+l][2]*pc{l]-H:p[n+2][2]*pc{2]+cp[n+3][2]*pcl3]; 

} /* end of cardsp*/ 


/* Cross_product,dot_product and vector normalization */ 


Dot_product ******************/ 
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/* Input the points of vector A and B 
/* Return value : cos( theta) 


*/ 


*/ 


#include <mathii> 

float dot_product( pva,pvb ) float *pva,*pvb; 

{ 

float cos_theta; 

normalize(pva); 

normalize(pvb); 

cos_theta=(*(pva)**(pvb) + *(pva+l)**(pvb+l) + *(pva+2)**(pvb+2)); 
retum(cos_theta); 

) /* end of dot_product */ 

/m «*m^****«m Normalize******************/ 

/* Input the point of the vector which will be normalized */ 

normalize(pv) float *pv; 

{ 

float m; 

m= sqrt(*(pv)**(pv) + *(pv+l)**(pv+l) + *(pv+2)**(pv+2)); 

*(pv) / =m; 

*(pv+l) A=m; 

*(pv+2) /=m; 

} /* end of vector_normalize */ 

/***»***« »**•* Cross_product **************/ 
croBS_product(pvc,pva,pvb) 
float *pva,*pvb,*pvc; 

{ 

normalize(pva); 

normalize(pvb); 

*pvc » *(pva+l) * *(pvb+2) - *(pva+2) * *(pvb+l); 

*(pvc+l) = *(pva+2) * *pvb - *pva * *(pvb+2); 

*(pvc+2) = *pva * *(pvb+l) - *(pva+l) * *pvb; 
if (*pvc >= -0.001 && *pvc <= 0.001 ) *pvc = 0.0; 
if (*(pvc+l) >= -0.001 && *(pvc+l) <= 0.001 ) *(pvc+l) = 0.0; 
if (*(pvc+2) >= -0.001 && *(pvc+2) <= 0.001 ) *(pvc+2) = 0.0; 

if ((fabs(*pvc) <=0.01 ) && (fabs(*(pvc+l))<=0.01) && (fabs(*(pvc+2))<=0.01)) 
if (fabs(*(pva))<=0.01) 

{ 

*pvc = -1.0; 

*(pvc+l)= 0.0; 

*(pvc+2)=0.0; 


{ 

if (fabs(*(pva+l)) <= 0.01) 

{ 

*pvc = 0.0; 
*(pvc+l)= -1.0; 
*(pvc+2)=0.0; 


} 

normalize(pvc); 
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} /* end of crosa_product *1 

Model Depictor 

♦include <usercore.h> 

♦include <sun/fbioJi> 

♦include <math Ji> 

♦include <stdio Ji> 

♦include 7home/Cai/demolib.h" 
♦include "model. h" 


int nverMipoly; 

float bboxf3][2]; 

float planeq[MAXPOLYI4]; 

float vertices [MAXVERTI3]; 

float n_vert[MAXVERT][3]; 

float normal[MAXVERT][3]; 

short my color l,cindex[MAXVERT]; 

int npvertfMAXPOLY]; 

int *pvertptr[MAXPOLY]; 

short planJnfo[MAXPOLY]; 

int pvertfMAXFVERT]; 

float dxlistt 100] > dylist[ 100],dzlist( 100]; 

int indxlist(100l,nuinber; 

float red[256],grnI256],blu[256],dred[256],dgm[256] t dbluf 256]; 

int rmm ahade color.num shade Jevel; 
float forward; 

int vertex_index[6][4]=Ul,2,6,5), 

{2,3, 7, 6}, 

{3,4,8,71, 

{4,1,5,81, 

{1,4,3,21, 

{5,6,7,811; 


f* line variables */ 

int nline,nlvert4ine_index_list{ MAXLVERT] ,npline[MAXLINE], 

short clinefMAXLINEl; 

float line_vert(MAXLVERT][3]; 

float xmax,xmin,ym ax, ynrin,zmax,zmin,xcent,ycent,zcent,length,emin,emax, scale; 

float xlistt 100] ,y lis t[ 1 0 0] ,zlist( 100]; 


float menu_x,menu_y,menu_h,menu_w,menu_f,window_fact; 

int n_button,face_remove; 

float T3,L,R,dot, alphabets, r, 

int nobj,npoly,nvert; 

int free_vert,free_poly,free_obj,free_pvert; 

int ol divert, old_poly,old_p vert; 

int raster Jd; 

struct menu_table_format menu_table[30]; 

struct view_parametersjbrmat my_view_parameters; 

struct image_status_ format image_status; 
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struct obj *objlst; 
struct poly *polylst; 
struct vert *vertlst; 
struct pvert *pvertlst; 

VECTOR cut_normal; 

POINT plane Jocate; 

POINT cut_plane[6]; 

main(argc, argv) 
int argc; 
char *argvO; 

{ 

char stringt81]; 
int i,endj4^ 
int buttonjength =0; 
float gema,w,cosb,x,y,z; 

POINT cut_point; 

LINE line; 

FILE *fptr, 

printfTInput cut plane normal "); 

scanfT%f %f %f\&cut_nomal4i,&cut_normal.v,&cut_normal.w); 
if ( load_file(argvtl]) !* NULL ) exit(l); 


get_view_surface{our__surface > argv); 

start_up_core(); 

init_menu(); 


I* initialize */ 


r=5000.0; 
alpha=beta=0.0; 
face_remove=FALSE; 
window_fact» 1.0; 

num_shade_color=4,num_shade_level=52; 
image_s tatus .1 ine =i m age_s tatus . wire_.fr am e = 
image_status.out__box=image_status.solid= 

image J3tatus.special=image_status. cutting =image_status.move = FALSE; 

forward = 12.0; 
rasterjd = 900; 

planejocate jc =* 0.0;plane_locate.y = 0.0;planejocate.z = 0.0; 

/* 

image_status.out_box=TRUE; 

create_retained_segment( OUTBOX); 
draw_box(); 

close_retained_segment(OUTBOX); 

*/ 

new_view_point( alpha, b eta, r); 
for(n) 

{ 

switch(call_menu(MENU,TRUE)) 

{ 
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caM 1: /* Right +10 *f 

if (image_check()) break; 
alpha +-10.0; 

new_view _point(alpha,beta,r); 
image_8witch(); 

break; 

case 20: /* Right + 2 *f 

if (image_check()) break; 
alpha +*2.0; 

new_view_point( alpha, beta, r); 
image_s witch( ); 

break; 

case 2: /* Left -10 
if (image_check()) break; 

alpha - *10.0; 

new_view_point( alpha, be ta,r); 
imagejawitchO; 

break; 

case 21: f* Left -2 
if (image_check()) break; 
alpha - =2.0; 

new_view_point( alpha, beta,r); 

image_switch(); 

break; 


case 3:/* down +10 

if (image_check()) break; 
if (beta <90.0) beta +=10.0; 

new_view_point( alpha, beta,r); 
image_switch( ); 
break; 

case 23: f* down +2 

if (image_check()) break; 

if (beta <90.0) beta +=2.0; 

new_view_point(alpha,beta,r); 

image_s witchO; 


break; 

case 4: /* up +10 */ 

if (image_check()) break; 

if (beta > -90.0 ) beta -=10.0; 
new_vi ew_point( alpha, beta,r); 
imagers witchO; 


break; 
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case 22: /* up +2 */ 

if (image_check()) break; 

if (beta > -80.0 ) beta -—2*0; 

new_view_point(alpha,beta,r); 

image_switch(); 

break; 

case 5: 

if (image_check()) break; 

r +*2500.0; /* zoom out *f 

new_view_point( alpha, beta, r); 
image_switch(); 

break; 

case 6: 

if (image_check()) break; 

r -*2500.0; /* zoom in */ 

if (r< 100.0) r=100.0; 
newjview__point(alpha,beta,r); 
image_switch(); 

break; 

case 7: I* reset */ 

if (image_check()) break; 
r =5000.0; 
al pha=beta=0 .0; 
windo w_fact= 1 .0; 
setvwpo(x,y,z,bbox); 
new^view^poinUalpha^eta^); 
image_awitch(); 

break; 

case g. /************* solid **********************/ 

if (face_remove I I image_status.move) { 
newJrameO; 

face_remove = image_status.move = FALSE; 

} 

newjview.pointfalphajDeta^); 

setjshading_paraineters(.01,.9,.0,0.0,7.,Q,l); 

set_primitive_attributes(&PRIMATTS); 

set_polygon_interior_style( SHADED ); 

create_temporary_segment(); 

drawobjlO; 

face_remove=TRUE; 

close_temporary_segment(); 

break; 
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cu t ft********************/ 


case 9: /*****♦♦****♦•*•**♦* 

if (face_remove )( 

newJrameO; 
face_remove * FALSE; 

} 

if (image_statufl.wire_frame){ 

image_8tatus.wire_frame*FALSE; 
deleU_retained_segment( WIREFRAME); 

1 

if (image_8tatusiine){ 

image_8tatus.line=FALSE; 

delete_retained_segment(DRAWLINE); 

) 

if (image.status.out _box)l 

image_status.out_box=FALSE; 

delete_retained_segment(OUTBOX); 

} 

if (image_8tatu8.move I I image_status.cutting){ 
if (image_status .move) { 

delete_retained_segment(SCREENTEXT); 
image_statua.move = FALSE; 

} 

if (image_statm.cutting)image_status.cutting =* FALSE; 
new_frame(); 

delete j-etained_segmentf PLANE); 

) 


cut_obj(cut_normal,plane_locate); 

create_temporary_segment(); 

draw_box(); 

/*draw_wireframe();*/ 
close_temporary_segment(); 
crea tejr etai n ed_s egm e n t( PLANE ) ; 

new_cut_plane(cut_normal,plane_locate,&cut_plane[0].x > &i); 

draw_cut_plane(i); 

close_retained_segment(PLANE); 

rebuild_planeq(); 
image_status.cutting = TRUE; 

section_draw(); 

/*inquirejviewing_parameters(&myjview_parameters); 

setjview_reference_point(0.0,0.0>0.0);*/ 


break; 


case 10: 

shut_down_corel(); 

exitO; 

case 11: 
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if (image status. wire_frame) 

l 

image_status.wireJrame=FALSE; 

del ete_re tained_s egm ent( WIREFRAME); 

} 

else 

{ 

image_status.wire_frame=TRUE; 
new_view_point( alpha, beta, r); 
create_retained_segment( WIREFRAME); 

draw_wireframe(); 

close_retained_segment( WIREFRAME); 

) 

break; 

case 12: 

if (image_status.out_box) 

{ 

image_status.out_box=FALSE; 

delete_retained_segment(OUTBOX); 

) 

else 

l 

image_status.out_box=TRUE; 
new_view_point( alpha, beta, r); 
create_retained_segment( OUTBOX); 
draw_box(); 

close_retained_segment( OUTBOX); 

} 

break; 

case 13: 

if (image_status.iine) /* line */ 

{ 

image_status.line=FALSE; 

delete_retained_segment(DRAWLINE); 

} 

else 

{ 

i m age_s tatus . 1 i n e = TRUE ; 
new_view_point(alpha,beta,r); 
create_retained_segment(DRAWLINE); 
drawl ine(); 

cl os e_retained_segment( DRAWLINE); 

} 

break; 

case 14: /* larger*2 */ 

if (window_fact>0.1 ) windowjact -»0.1; else if (window_fact>0.01) windowjact -=0.01; 
new_view_point( alpha, be ta,r); 
image_switch(); 
break; 


case 24: /* larger*! */ 

if (window_fact>0.01) window_fact -=0.01; 
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new_view_point( alpha, beta,r); 
image_switch(); 
break; 


case 15: /* smaller*2 *1 

windowjact +*0.1; 
new_view_point( alpha, beta,r); 
imagers witchO; 

break; 


case 25: /* smaller* 1 */ 

window_fact +*0.01; 
new_view_point( alpha, be ta,r); 
image_switchO; 
break; 


case is: /* Switch cutting direction 

cut_normal.u = -cut_normal.u; 
cut_normal.v = -cut_normal.v; 
cut_normal.w = *cut_normal.w; 

break; 

case 19: /* Cutting plane move forward */ 

if (bbox[0][l]<=piane_locatejc && 

bbox[l][l]<=plane_locate.y && 
bbox[2][l]<=planeJocate.z ) break; 
move_cut_plane( 1); 
break; 

case 26: /* Cutting plane move backward */ 

if (bbox(0][0]>=planejocatejc && 

bbox{l][0]>=planejocate.y && 
bbox(2][0]>=plane Jocate.z ) break; 
move_cut_plane(-l); 
break; 

case 27: f* Raster file */ 

raster^fileO; 
break; 

case 28: f* Load Raster file •/ 

load_raster(); 
break; 

case 29: /* Move Raster */ 

move_raster(); 
break; 

case 30: /* Delete Raster */ 


re_store(); 

delete_raster(); 

break; 

case 31: /* Menu hiding */ 

set_segment_visibility( MENU, FALSE); 
set_echo( LOCATOR, 1 ,1); 
for (;;) 

{ 

do 


27 



await_anyjbutton_getjocator - 2(20000000 > l»&button,4x,fty); 

while(button ** 0); 

iflbutton—1 1 1 button =2 I I button 3) break; 

set_segment_visibili ty(MENU,TRUE); 
set_echo(LOCATOR,l,0); 

break; 

case 32: /* cutting file *1 

if (face_remove ){ 

newJrameO; 
face_remove * FALSE; 

} 

if (image_status.wire JrameH 

image_status.wire_frame=FALSE; 

del ete - retained_segment( WIREFRAME); 

} 

if (image_status.line){ 

imagers tatus .1 in e=F ALSE ; 

delete_retained_segment(DRAWLINE); 

1 

if (image_s tatus. out_box)( 

imagers tatus. out_box=FALSE; 
delete_retained_segment(OUTBOX); 

} 


if (image_status.move I I image jstatus. cutting)! 

if (image_status.move) image_status.move = FALSE; 
if (image_status.<nitting)imagej5tatus.cutting « FALSE; 
new_frame(); 

delete_retained_segment(PLANE); 

} 

if C(fptr = fopenCcutjfile", "r")) == NULL) ( 

printfTCan’t open cut_file fileW); 
exit (1); 

} 


fscanflfptr,"%d",&k); 

ere ate_r etain ed_s egm ent( PLANE) ; 
draw_box(); 

for(j=0*j<k3-M*H . _ . . 

fscanflfptr/ , 9fcf%f%f\&cut_normal.u,&cutjnorinai.v,&cut_normal.w); 

fscanflfptr, M %f%f%r,&plane_locate.x,&planeJocate.y,&planeJocate.z) 

plane Jocatejc =plane_locate.x*scale; 
planejocate.y =planeJocate.y*scale; 
planejocate.z =planeJocate.z*scale; 
cut_obj(cut_normal, plane Jocate); 

new_cut_plane(cutjionn al, plane jocate,&cut_plane[0].x,&i), 

draw_cut_plane(i); 

rebuild_planeq(); 

} 

close(fptr); 
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cl os e_r etained_segmen t( PLANE); 
image_status .cutting » TRUE; 

/*Bection_draw();*/ 


break*/* end of cut files */ 

default: 

break; 


} /* end of switch */ 

} /* end of waiting for */ 

} f* end of main *f 

int image_check() 

if (image_status.move I I image_s tatus. cutting) ( 
delete_retained_segment(PLANE); 
newJrameO; 

image_status.cutting = image_status.move = FALSE; 

1 


if (image_status.cutting) { 

newJrameO; 

image_status.cutting = FALSE; 

if ( image_s tatus .line FALSE && image_status.wire_frame 

— FALSE && image_status.out_box== FALSE && image_status.solid 
== FALSE && image_status. cutting == FALSE) retum(TRUE); 

else 

retum(FALSE); 


} 


new_view_point( alpha, b eta, r) 
float alpha,beta,r, 


float cosb,x,y,z; 

cosb= fabs(cos(beta*PI/180.0)); 
x = r*cosb*sin(alpha*PI/180.0); 
y = r*sin(beta*PI/180.0); 

2 = r*cosb*cos(alpha*PI/180.0); 
setvwpo<x,y,z,bbox); 

I* set_single__window; */ 

} /* end of new_view_point *1 


move_cut_plane( direction) int direction; 
{inti; 

char buffI120]; 


if (image_status.cutting){ 

delete_retained_segment(PLANE); 
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image_8tatus.cutting = FALSE; 

} 

if (image_Btatus.move»* FALSE && 

facejremove — FALSE && image jjtatus.wire_frame 
create_temporary_segment(); 

draw_box(); 

/*draw_wireframe();*/ 


FALSE ){ 


close_temporary_segment(); 

create_retained_segment(PLANE); 

} 

if (image_Btatus.move** FALSE && 

(facejremove I limage_status.wire_frame )){ 
create_temporary_segment(); 

draw_box(); 

close_temporary_segment(); 
create_retained_segment( PLANE); 

plane Jocatejc +* forward*cut_normal.u * direction ; 
planejocate.y +« forward* cut_normal.v * direction ; 
plane Jocate.z +» forward*cut_normai.w * direction ; 

new_cut_plane(cut_normal,planeJocate,&cut_plane[0].x,&i); 

if (image_status.move ){ 

delete_retained_segment(PLANE); 

create_retained_segment(PLANE); 

delete_retained_segment(SCREENTEXT); 

1 

draw_cut_plane(i); 

close_retained_9egTnent(PLANE); 

inquire_viewing_parameters(&my_view_parameters); 

set_view_reference_point(0.0,0.0,0.Q); 
set_view_plane_normal(0.0,0.0, -1.0); 
set_view_plane_distance(0.0); 

s et_proj e ction( PARALLEL , 0 .0 , 0 . 0 , 1 . 0 ) ; 
set_view_up_3(0.0,1.0,0.0); 
set_window(0.0, 100. 0,0. 0,767.0); 
setjview_depth(0.0,1.0); 
set_window_clipping(FALSE); 
set_viewport_3(0.0, 0.9 15, 0.0, .75,0.0,1.0); 

create ^retain ed_segmen t( SCREEN TEXT); 

move_abs_3(10.0, 10.0,0.0); 
textf’Cutting Plane Location: "); 

sprintfibuff;%3.3f %3.3f %3.3f\ 4 , 

planejocate.x/scale, plane Jocate.y/scale, plane Jocate.z/scale); 

move_abs__3(35.0, 10.0,0.0); 
text(buff); 

close_retained_segment(SCREENTEXT); 


set_viewingj?arameters(&my_view_parameters); 


image_status .move* TRUE; 
}/* end of move_cut_plane */ 
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iinclude "mycut" 

raster_file{) 

{ 

intrasfidjj; 

int butnum, button; 

float xmin, ymin, x, y,z, xmax, ymax.zmax.mx.my; . , 

float wx, wy, wz.wwx, wwy,wwz,kx,ky,xmin_ndc,xmax_ndc,ymin_ndc,ymax_ndc; 

char string[80]; int length; 
struct suncore.raster my .raster; 
extern struct vwsurf *our_surface; 
struct { int type, nbytes; 

char *data; ) map; 

FILE *fptr, 
char *path_name; 

path.name * "/home/Cai/load_fileJib ; 

inquire_viewing_parameters(&my_view_parameters); 

set_view_reference_pointf0.0,0.0,0.0); 


set_view.pl ane_normal(0. 0,0.0, -1.0); 
set_view_plane_distance{0.0); 
set_projection(PARALLEL,0.0,0.0,1.0); 
set_view_up_3(0. 0,1. 0,0.0); 
set_window(0.0, 100.0, 0.0, 767.0); 
set_view_depth(0.0,1.0); 
set_window_clipping(FALSE); 
set_viewport_3(0. 0,0.915, 0.0, .75, 0.0, 1.0); 


set_echo(LOCATOR, 1 ,1); 
for (;;) 

{ 

await_any_button_get_locator_2(1000000,l,&button,&mx,&my); 

while(button == 0); 
break; 


} 

set_echo_position( LOCATOR, 1, mx,my); 

print£C\nfirst point %f %f\mx,my); 
map_ndc_to_world_3( x, y, 0.0, &wx, &wy,&wz); 
printfTxyz %f %f %f\wx, wy, wz); 
set_echo( LOCATOR,!, 6); 

do{ 0 , 

await_any_button_get_locator_2( 1000000,1, &butnum, &xmax, &ymax); 

} while (butnum != 3); 
printfCXnfirst point %f ,xmax, ymax); 
printf("\n%f %f %r,bboxt0][0],bbox[ll[0],bbox[2][0]); 
printfT\n%f %f %F,bbox[0][l],bbox[l][l],bbox[2][l]); 

set echo position KEYBOARD, 1, 0.25, 0.1); /* move to start positn */ 
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/* set user char buf size */ 


setJteyboarcK 1, 80, 0); 

Tmin ndc a (mxcaxmax ) ? mx : xmax; 
xmaxjndc * (xmax>» mx) ? xmax : mx; 
ymin.ndc ■ (my<*ymax ) ? my : ymax; 
ymax_ndc a (ymax>* my) ? ymax : my; 

size^astertour_surface,xmin_ndc,xmax w ndc,ymin_ndc,ymax_ndc,&my_raster); 

allocate_raater(&my_raster); 

get_raster(our_surface,xmin_ndc,xmax_ndc,ymin_ndc,yinax_ndc,0,0,&niy_raster); 

create_r etained_se gment( S CREENTEXT) ; 

move_abs_2(10.0, 10.0,0.0); 
text* "Raster File Name: "); 

set_echo( KEYBOARD, 1,1); f* echo the text *f 

await_keyboard( 10000000,1, string, ^length); 

if flength && stringflength-l] = \n ) stringOength-1] = W; 


move_abs_2(25.0, 10.0,0.0); 
text ( string); 
move_abs_2( 50. 0,50.0); 

close_retained_segment(SCREENTEXT); 

set jmage_transformation_type( XLATE2); 
create jretained_8egment(99); 
setj?rimitive_attributes( &PRIMATTS); 
setjpick_id( 1); 
put_raster(&myjraster); 
set_segment_detectability( 99, 5); 
if( (rasfid = open( "try", 0)) == -1) ( 

rasfid a creaK "try", 0755); 

) 

if (rasfid != -1) { 

map. type a 1; map mbytes = 768; map. data = 85 ; 
file _to_raster( my.raster, &map, rasfid ); 
clos e( rasfid); 

) 

free_raster(&my_raster); 

close_retained_segment(99); 
strcat(path_name, &string{0]); 
save_segment(99, string); 

delete_retained_segment(99); 

delete_retained_segment(SCREENTEXT); 

set_viewing_parameters(&my_view_parameters); 

)/* end of raster_file */ 

load_raster() 

{ 

intrasfidjj; 
int butnum, button; 

float xmin, ymin, x, y,z, xmax, ymax,zmax,mx,my; 

float wx, wy, wz,wwx, wwy,wwz,kx,ky,xmin_ndc,xmax_ndc,ymin_ndc,ymax_ndc; 
char string[80]; int length; 
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struct auncore_raster my_raster, 
extern struct vwsurf *our_surface; 
struct { int type, nbytes; 

char *data ; ) map; 
FDLE*fptr; 


int segnam, pickid, segtype; 
float sxO, syO, angO, txO, tyO; 
float sx, sy, ang, tx, ty; 
float px, py; 

char *path_name; 


inquire_viewing - parameters(&my_view_parameters); 

set_view_reference_point(0.0,0.0,0.0); 

set_view_plane_normal(0. 0,0.0, -1.0); 
s e t_vi e w_p i ane_dis tan c e( 0 .0 ); 

se t_proj ection( PARALLEL , 0 . 0 , 0 . 0 , 1 .0 ) ; 

set_view_up_3(0.0,1.0,0.0); 

set_window(0.0, 100.0, 0.0, 767 .0); 
set_view_depth(0.0,1.0); 

set_window_clipping( FALSE); 
set_viewport_3(0.0, 0.915, 0.0, .75, 0.0, 1.0); 

set_echo_position( KEYBOARD,!, 0.25, 0.1), 
set_keyboard( 1, 80, 0); 

create_retained_segment(SCREENTEXT); 

move_abs_3(10.0, 10.0,0.0); 
texK "Raster File Name: H ); 

await_keyboard( 10000000,1, string, ^length); 
set_pick_id( 1); 

if Qength && string(length-l] = An) string{length-l] = \0 ; 

move_absJ3(25.0, 10.0,0.0); 

textf string); 

close_retained_segraent(SCREENTEXT); 

if ((fptr =* fopen( string, V)) == NULL) { 

move_abs_3(25.0, 20.0,0.0); 

text( "Can’t open raster file "); 
for (i=0d<10000;i++); 

dele te_retained_segment( SCREEN TEXT); 


restore.segmentfraster^id, string); 


set_echo(LOCATOR, 1 ,1); 
for (”) 

I 

await_any_button_get_locator__2( 1000000, l,&button,&mx,&my) 
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while(button =* 0); 


break; 

l 

set_echo( LOCATORS 0); 

do { 

await_any_button_get_locatorJ2( 0,1, &butnum, &x, &y); 
tx = x-mx; ty = y-my; 


set_segment_image_t rans l a kG-2(raster_id , tx,ty); 

} while (butnum != 3); 
raster Jd-M-; 

set_viewing_parameters(&my_view_parameters); 

delete^r e tai n ed_s e gm e nt( SC RE EN TEXT); 

}/* end of load raster */ 

move_rasterO 

{ 

int rasfiddj; 
int butnum, button; 

float xmin, ymin, x, y,z, xmax, ymax,zmax,mx,my; 

float wx, wy, wz,wwx, wwy,wwz,kx,ky,xmin_ndc^niax_ndc,ymin_ndc,yniax_ndc; 

char string[80]; int length; 

int segnam, pickid, segtype; 
float sxO, syO, angO, txO, tyO; 
float sx, sy, ang, tx, ty; 
float px, py; 

/* inquire jviewing_parameters(&my_view_parameters); 
set_view_reference_point(0.0,0.0,0.0); 


set_view_plane_normal(0.0,0.0, -1.0); 
set_view_plane_distance(0.0); 
s et_proj e ction( PARALLEL , 0 . 0 , 0 . 0 , 1 . 0 ) ; 
set_view_up_3(0. 0,1.0, 0.0); 
set_window(0.0, 100.0, 0.0, 767.0); 
set_view_depth(0. 0,1.0); 
set_window_cli ppi ng( FALSE) ; 
set_viewportJ3(0.0, 0.915, 0.0, .75, 0.0, 1.0); 

create_retained_segment( SCREEN TEXT) ; 
move_abs_2(10.0, 10.0); 
textTFind your image and move ’); 

close_retained_segment(SCREENTEXT); 

await_pick( 1000000000, 1, &segnam, &pickid); 
await_any_button_get_locator_2( 0,1, &butnum, &px, &py), 

set_echo( LOCATOR, 1, 0); 
do { 

await_any_button - get_locator_2( 0,1, &butnum, &x, &y); 

tx = x-px; ty = y-py; 


set_segment_image_translate_2(segnam , tx,ty); 
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} while (butnum !■ 3); 

set_vie wing_p aram e ter8( &my_vi e w_param e ters ) ; 

delete_retained_segment(SCREENTEXT); 

*/ 


delete_raster() 

{ 

new_frame(); 

} 


section draw() 

{ 

intijjc; 
int index^flag; 

float vecl0[3] f vecfcl[3],scale f plane[8][3],area,cut - area,pointl400][3],u,v*w; 

FILE *fptr; 

float x,y,z,xc,yc,zc,sina,cosa,sinb,cosb, minx, maxx,minz,maxz > fact_x,fact_z, fact; 

if ((fptr = fopenCnc.dat", V')) == NULL) { 

printfTCan't open nc.dat fileW); 
exit(l); 

) 

fscanf(fptr > "%f%f%f\&vect0[0],&vect0tl],&vect0[2]); 

fscanflfptr, M %f,&scale); 

fscanflf^tr, n %d M ,&k); 

u=vect0[0]; 

v=vectO[l]; 

w=vectG[2]; 

sinb» (u— 0.0 && v=0.0) ? 0.0 : -u/sqrtfu*ii+v*v) ; 
cosb* (u==0.0 && v==0.0) ? 1.0 : v/sqrt(u*u+v*v); 

CO sa= (u==0.0 && v==0.0 && w==0.0) ? 1.0 : sqrt(u*u+v*v)/sqrt(u*u+v*v+w^w); 
sina* (u«»0.0 && v==0.0 && w==0.0) ? 0.0 : w/sqrt(u*u+v*v+w*w); 
maxx=minx=maxz=miriz=0.0; 

for (i=0;i<k;i++) 

fscanflfptr/^f^f^r^planeliltOl^planetiltll^&plan^nW]); 
plane{i][0] = plane[i][0Vscale; 
plane(i][l] = plane(i][l]/scaJe; 
plane{i][2] = plane[i][2]/scale; 

) 


for (i=0?<kp++) 

{ 


if(i— 0M 

xc=plane[i][0]; 

yc=plane[i][l]; 

zc=plane(i][2]; 
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} 

xc=yc=zc=0.0; 

plan e(i][ 0 }= plan e{i][0]-xc; 

plane(i][l]=plane{i][l]-yc; 

plane(il 2 ]=plane(i][ 2 ]-zc; 

xaplane[i]I 0 ]*co 8 b+plane[i]ll]*sinb; 

y«-plane(i]I 0 ]*sinb+plane[i][l]*co 8 b; 

z= plane(iI2]; 

planeMO] = x^lanetill] = y;plane{il[2] = z; 
x=plane(ilO]; 

y=plane(i][l]*cosa+planeli][2]*8ina; 

z=-plane{i][l]*sina+plane[i][2]*cosa; 

plane(i][0] = x;plane{ill] = y^lane{i]I2] = z; 
plane(i][0]=planetiIO]+xc; 
plane[i][l]=plane{il[l]+yc; 
plane(i][2]=planeti][2]+zc; 

printfT %3.4f %3.4f %3.4f\n'\plane[i][0],plane[i][l],planeUJl2J), 

ifli==0)( 

maxx=minx=plane[i][0]; 

maxz=minz=plane[il[2]; 


maxx= (maxx<plane[i][0]) ? plane[i][0] : maxx; 
minx= (minx>plane[i][0]) ? plane[i][0] : minx; 
maxz= (maxz<plane[iI21) ? plane[i][2] : maxz; 
minz= (minz>plane[i][2]) ? plane{i][2] : minz; 

} 


flag = (fabs(maxx-minx)>=fabs(max 2 -minz)) ? TRUE : FALSE; 
f ac t= (flag) ? fabs(5.4/(maxx-minx)) : fabs(5.4/(maxz-minz)); 

inquire_viewing_parameters(&my_view_parameters); 

set_view_reference_pointf0.0,0.0,0.0); 


set_view_planejiormal(0.0,0.0, -1.0); 
set_view_plane_distance(0.0); 

set_projection(PARALLEL,0.0,0.0,1.0); 

setjvdew_up_3(0.0,1.0,0.0); 

set_window<0.0,6.0,0.0,6.0); 

set_view_depth(0.0,1.0); 

set_window_clipping(TRUE); 


/*set_viewport_3(0.0, 0.375, 0.375, 0.75, 0.0, 1.0);*/ 

set_viewport_3(0.0, 0.5, 0.25, 0.75, 0.0, 1.0), 
create_temporary_segment(); 

/*xlist{ 0 ]= 0 . 0 ; 2 iist[ 0 ]= 0 . 0 ; 

xlist(l]=6.0*4ist[l]=0.0; 

xlist(2]=6.0^1ist(2]=6.0; 

xlistf3]=0.0;zlistf3]=6.0; 

xlist(4]=0.0;zlistt4]=0.G; 

move_abs_2(xlist(0],zlist(0]); 

polyline_abs_2(xlist,zlist,5);*/ 

setjinewidth(0.5); 

setJine_index(MENU_TEXT_COLOR); 
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for (j«Oa<kj++) 


l 

if (flag){ 

xlistjj] a (planeU][0]-minx)*fact+(6.0-(maxx-minx)*fact)*0.5; 
zlist(j] * (plane{j][2]-minz)*facM6.0^maxz-minz)*fact)*0.5; 

1 

else { 

zlistfj] = (plane(j][0]-minx)*fact+(6.0-<maxx-minx)*fact)*0.5; 
xlistlj] » (plane{j][2]-niinz)*facU<6.0-(maxz-minz)*fact)*0.5; 

1 

} 

xlistjj] » xlistfO]; 
zlistjj] a zlistf 0]; 
move_abs_2(xlist(0],zlist[0]); 
polyline_abs_2(xlist,zlist,k+ 1); 


cut_area«0.0; 


while ( fscanfTfptr/ , %d , *,&k)!=EOF) { 

fscanftfptr/'^f%f%r f &vect0[0] ( &vect0[l],&vect0[2]); 

fscanf(f^tr, , '%d",&k); 

for (i= 03 <lp++) 

{ 

fscamTfptr, M %f%f%r,&point[i][0],&point[i][l],&point(i]I2]); 

point{iI0] = point(i][0yscale; 

pointOIl] a point(i][iyscale; 

point[il2] = point[i][2yscale; 

xc=yc=zc=0.0; 

point(i][0]=point{i][0]-xc; 

point(iIl]apoint[i][l]-yc; 

pointfil2]=point(i][2]-zc; 

x=point(i][0]*cosb+pointf i][ l]*sinb; 

y= -point[i][0]*sinb+point[i][ l]*cosb; 

z= pointf iI2]; 

pointfilO] = x^point[i][l] = y;point{il2] = z; 
x=point[i][0]; 

y=pomtjj][l]*cosa+pointfi][2]*sina; 
z= -pointfi][l]*sina+point[i][2]*cosa; 

pointftlO] = x^Doint[i][l] = y;point[iI2] = z; 
pointjil0]=point[i][0]+xc; 
point[i I l]=point[i][ 1 ]+y c; 
pointfiI2]=point{i][2]+zc; 

/*printfT %d %3.4f %3.4f %3.4An*\i,point[i][0],point[i][l],point[i][2]);*/ 
if (flag){ 

xlistji] = (point[i][0]-minx)*fact+(6.0-(maxx-minx)*fact)*0.5; 
zlistji] = (point[i][2]-minz)*fact+<6.0-(maxz-minz)*fact)*0.5; 

1 

else { 

zlistfi] = (point[i][0]-minx)*fact+(6.0*(maxx-minx)*fact) ,# 0.5; 
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xlistti] * (poinUi][2]-minz)*fact+(6.0-(max2-min2)*fact)*0.5; 


} 

} 

xlistti] * xlistfO]; 
zlistfi] m zlisttO]; 
move_abB_2(xlist{0],zlist(0]); 
polyline_aba_2(xlist > zlist,k+ 1); 


}/* end of cut file */ 
setJinewidth(O.O); 
close_temporary_segment(); 

se t_vi e win g_p aram e te rs ( & my _vi e w_p ar am e t e rs ) ; 
close (fptr ); 

}f* end of section draw *1 


re_store() 

{ 

int ij,k4ndex,numjvert,polyJndex,new_vert,new_index,new_nobj; 
float x,y,z,xO,yO,zO, length, bboxl[3][2]; 
short *ptr,*vert_indexjst; 
struct vert *vertlstptr; 

FILE *fptr, 


new_nobj=new__vert=0; 

vert_index_lst = (short *)malloc(free_vert * si zeoft short)); 
for (i=0;i<free_vertd++) 

vert_index_lst[i]= -1; 
fptr = fopenCnn.dat", "w"); 

for (i=0d<nobjd++) if (objlst{i],npoiy !=0) new_nobj++; 
npoly=0; 

for (i=sO;i<nobjd++) npoly=npoly +objlstli].npoly; 

for (i=0;i<nobjd++)t 

poly_index = objlstfi].index; 
for (j= 03 <objlstIi].npoly-j++){ 

linking<&poly_index); 
ptr = polylst(poly_index].pvert.ptr; 
num_vert = polylst[poly_index++].npvert; 
for (k=0;k<num_vert;k++){ 

vert_index_lst[*ptr++]=-2; 

)/* mark all the vertexes of one polygon */ 

)/* end of j poly loop *f 
}/* end of i loop *f 
for (k=0;k<free_vert;k++) 

if (vertJndexJst[k]»-2) new_vert++; 


bboxl[0][0] = bboxl[0][l] = bboxl[l][0] = 
bboximtl] = bboxl[2][0] = bboxl[2][l] = 0.0; 


38 



for (i-0;i<free_vert;i++){ 

if(vertjndexjst(i]=« -2)1 

if (bboxl[0][0] >- vertlst(i].vertexjc) bboxl[0][0] 

= vertlstfil.vertexjc; 

if (bboxl[0][l] < vertlst(i].vertexjc) bboxlfOll] 

= vertlstli]. vertex .x; 

if (bboxl[l][0] >= vertlstfi].vertex.y) bboxl[110] 

= vertlst(i].vertex.y; 

if (bboxlUUl] < vertlst(i].vertex.y) bboxlUIl] 

= vertlstli].vertex.y; 

if(bboxl[2][0] >= vertlst{i].vertex.z) bboxl[2][0] 

= vertlstti].vertex-z; 

if (bboxl[2][l] < vertlst(i].vertex.z) bboxl[2][l] 

= vertlst(i].vertexj; 

) 

vert_index_lst{i]= -1; 

) 

fprintRfptr,"%d %d %d\n",new_nobj,npoly,new_vert); 

for (i=0;i<3d++) for(j=0'j<2y++) bboxl[i][j] * = 1.2; 

fprintflfptr, "%f %f %f %f %f bboxl[0][0]/scale, bboxl[0][lVscale, 

bboxl[l][01/scaie ( bboxl[l][l]/scale,bboxl[2][0]/scale,bboxlt2][lV8caie); 

new_vert=0; 
for (i= 0 p<nobj 3 ++)i 

if (objlsttil.npoly ==0) goto next_obj; 
poly_index = objlst(i].index; 
for (k=0;k<free_vert;k++) 

vert_index_lst[k]= -1; 

for (j=Oj<objlstti].npolyj++){ 

linking(&poly_index); 
ptr = polylsdpoly_index].pvert.ptr; 
num_vert = polylst[poly_index++].npvert; 
for (k=0;k<num_vert;k++){ 

vert_index_lst[*ptr++]=-2; 

)/* mark all the vertexes of one polygon */ 

}/* end of j poly loop *1 

num_vert=0; 

for (k=0*4c<free_vert;k-H-) 

if (vertJndex_lstflc]=*-2) num_vert++; 


fprintfffptr/’%d %d\n^objlsUi].npoly,num_vert); 
for (k=0;k<free_vert;k++){ 

if (vert Jndex_lst[k]*«-2){ 

x=vertlstfk].vertexjc/scale; 

y=vertlst[k].vertex.y/scale; 

z=vertlst[k].vertex.z/scale; 

fprintflfptr, ”%f %f %f\n'\x,y,z); 

vert_i n dex J s t[k] *ne w_v ert; 

new_vert++; 

} 

*/ 


}f* end of sending 
polyjndex = objlst[i].index; 
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for (j=0g<objl8t(i].npolyj++){ 

linking(&polyJndex); 
ptr ■ polylst{polyJndex].pvert.ptr; 
num_vert * polylst[poly_index].npvert; 
fprintflfptr,"%d ",num_vert); 
fprintfffptr,"%d ",polylst(polyJndex++].info); 
for (k=0;k<num_vert;k-M-){ 

newJndex*vert_indexJstI*ptr++]+l; 
fprintflfptr,"%d ",new_index); 

}/ * send all the vertexes index of one polygon *f 
fprintflfptr/An"); 

}/* end of j poly loop */ 
next_obj:k=0; /• no poly !! *f 
}f* end of i obj loop */ 

free(vertjndexjst); 
fprintflfptr > M %d %d\n*\0,0); 
fclose(fptr); 

} /* end of restore */ 

Subroutines for Model Depiction 

#include <sun/fbio.h> 

#include <usercore.h> 

#include <math.h> 
finclude <stdio.h> 

#include " model. h" 


/* Externals */ 

extern nobj,npoly,nvert,vertex_index[6][4]; 
extern free_vert,free_poly ,free_obj ,fr ee_p vert; 
extern ol divert, old_poly,old_pvert; 
extern float bbox[3][2]; 

f* Externals */ 

extern struct obj *objlst; 
extern struct poly *polylst; 
extern struct vert *vertlst; 
extern struct pvert *pvertlst; 


int plane j:ut_line(normal, locate, line,cut_point) 
VECTOR normal; 

POINT locate; 

LINE line; 

POINT *cut_point; 

{ 

float px,py,pz,d,dLd2, cosaO,cosal, delta; 
float a[3][3]; 

VECTOR v0 9 vl,l,vtemp; 
int endO_in,endl_in; 

float dot_product(); 

vO.u = line.pO jc - locates; 
vO.v = line.pO.y - locate.y; 
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vO.w m line.pO.z - locates; 
vLu * line.pl jc - locates; 
vl.v ■ line.pl.y - locate.y; 
vl.w ■ line.pl.z - locates; 

coeaO * dot_product(& normal, &vO); 
co sal * dot_product(&nonnal,&vl); 


if (cosaO> PLANETHICK && cosal>PLANETHICK ) retum(NOCUT); 

endO in * (cosaO <= PLANETHICK && cosaO >= -PLANETHICK) ? TRUE :FALSE; 

endl in - (cosal <- PLANETHICK &.&. cosal >= -PLANETHICK) ? TRUE JALSE; 


if (endl Jn && 

endOJn) return! ONPLANE); 
if ((co«aO< -PLANETHICK && cosal< -PLANETHICK ) 
I l( endOJn && cosal< -PLANETHICK) 

I I ( endl jn && cosaO< -PLANETHICK ) ) 
return! ALLCUT) ; 

if ((endOJn && cosal >PLANETHICK) 

I I (endl Jn && cosaO>PLANETHICK)){ 
if (endOJnM /* end at pO */ 

cut_point->x = line.pOjc; 
cut_point->y = line.pO.y; 
cut_point->z = line.p0.z; 
retum(ENDCUTO); 

} 

if (endl Jn){ /* end at pi */ 

cutjx>int->x m line.pljc; 
cut_point->y = line.pl.y; 
cut_point->z = line.pl.z; 
return! ENDCUTl); 

} 


} 


l.u = line.pl jc - line.pOjc; 
l.v = line.pl.y - line.pO.y; 
l.w = line.pl x - line.pO.z; 

delta * ((line.pl jc + 10.0 la line.pOjc) && (line.pljc + 10.0 != line.pOjc) 
&& (line.pla + 10.0 != line.pOjc)) ? 10.0 : 20.0; 
vtemp.u = line.pljc + delta - line.pOjc; 
vtemp.v = line.pl.y + delta - line.pO.y; 
vtemp.w * line.pl.z + delta - line.pO.z; 


cross_product(&vO,&vtemp,&l); 
cross jproduct(&vl,&vO,&l); 

a[0][0] = normal.u; a[0]U] = normal.v; a[0][2] = normal.w; 
a[l][0] = vO.u; a[l][l] » vO.v; all][2] = vO.w; 
a[2][0] * vl.u; a[2][l] = vl.v; a[2][2] = vl.w; 

d = normal.u*locate.x+normal.v*locate.y+normal.w*locate.z; 

dl = v0.u*line.plJC+v0.v*line.pLy+v0.w*line.pl.z; 
d2 = vl.u*line.pljc+vl.v*line.pl.y+vl.w*line.pl.z; 

matinv(&aI0][0]); 
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cut_poinfc->x » a[0][0]*d + a{0][l]*dl + a[0][2]*d2; 
cut_point->y * a[l][0]*d + a[l][l]*dl + a[l][2]*d2; 
cut_point->z ■ at2][0]*d + a(2][l]*dl + a[2][2]*d2; 
if (co»aO<- PLANETHICK) return(PAKTCUTl); 
if (cosaK- PLANETHICK) retum(PARTCUTO); 


rebuild_planeq() 

{ 

struct count{ 

short n 

}; 

struct count ^normal count; 
float x t y,zjength; 

int ij,vl p v2,v3 > vtmp,poiyJndexjj,nuni_vert; 

short *ptr; 

if (normalcount != NULL ) free(normalcount); 
normalcount= 

(struct count *)malloc(free_vert * sizeoflstruct count)); 
for (i=0;i<free_vert;i++) ( 

normal co unt(i].n = 0; 
vertlst{i].normal jc = vertlst(i].normal.y 
= vertlstfi] .normal j . =0.0; 

) 

/* rebuild planeq A3 ? C,D and normal at ehch vertex */ 

for (i=0' f i<nobj4++){ 
polyjndex = objlst(i].index; 
for (j=Oy<objlst(i].npoly*j++){ 

linking(&poiyJndex); 

ptr = polylsttpoly Jndexj.pvert.ptr; 
numjvert = polylst[polyJndex].npvert; 
polylst(poly_index].planq.a = poiylst[polyJndex].planq.b 
= polylst{poly_index].planq.c = polylst[poly_index].planq.d = 0.0 
v3 = *ptr++; 
v2 = *ptr++; 
vl = *ptr, 

for (jj = 0; jj < 3; jj++) 

( 

polylstJpoly_index].planq.a += vertlst{vl].vertex.y* 

(vertlst(v2].vertex.z - vertlst[v3].vertex.z); 
polylst{polyJndex].planq.b += vertlst[vl].vertex.x * 

(vertlst[v3].vertex.z - vertlst[v2].vertex.z); 
polylst{poly_index].planq.c += vertlst[vl]. vertex jc * 

(vertlstiv2].vertex.y - vertlst(v3].vertex.y); 
polylst{polyJndex].planq.d += vertist[vl].vertex.x * 

((vertlst{v3].vertex.y * vertlst{v2].vertex.z) - 
(vertlst(v2].vertex.y * vertlst(v3].vertex.z)); 
vtmp = vl; vl = v2; v2 = v3; v3 = vtmp; 

} 

x = polylstfpoiy_index].planq.a; 
y = polylst[poly_index].planq.b; 
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z m polylst[poly_index].planq.c; 
length ■ sqrK x*x + y*y + z*z>, 
ptr =» polylstlpoly JndexJ.pvert.ptr; 
for (jj * 0; jj < polylst(polyJndexl.npverti++) 
( 

vtmp* *ptr++; 

vertlst(vtmp] Jiormal jc +- x/length; 
vertlst(vtmp].normal.y +=» y/length; 
vertlst(vtmp].normal.z +- z/length; 
normalcount(vtmp].n++; 

} 

poly_index++; 

) /* end of poly, j loop */ 

} l* end of obj. i loop *1 


for (i = 0; i < free_vert; i++) 

{ 

if (normalcountfij.n != 0) { 

vertlst[i].normaljc/= normalcount(i].n; 
vertlst{i].normal.y /= normalcount[i].n; 
vertlst[i].normal.z /= normalcount(i].n; 
1 


/* Free the temporary storage */ 
free( normal count); 

} /* end of rebuild„planeq */ 


new_cut_plane(cut_normal, plane Jocate,new_planejptr,n_vert_ptr) 

VECTOR cut_normal; 

POINT planejocate; 

POINT *new_plane_ptr; 
int *n_vert_ptr; 

< 

int i jX n - vert >n_cut, search, plane_index,kl f flag; 

intkkjj; 

POINT box(8],new_point; 

LINE cutjine; 

short tempJndexO,temp_indexl,e_index,pass; 
struct { 

float x; 
float y; 
float z; 

short e_index; 
short vJndexO; 
short vjndexl; 
short pjndex; 
short flag; 
leut _point[12]; 

VECTOR vtemp; 
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POINT new_plane[6]; 
float dot_product(); 
float coaa[4]; 


box[0]jc * bboxtOlO]; boxtOj.y a bbox£l][0]; box[0].z = bboxt2I0]; 
box(l]jc = bbox[0ll]; baxtlj.y » bbox(l][0]; box[l].z » bboxt2I0]; 
box{2]jc ■ bbox[0][l]; boxt2].y - bboxtl][0]; boxl2].z - bboxf2ll]; 
box(3]jc * bboxf0][0]; box{3].y a bboxtl][0]; box(3].z = bbox[2ll]; 

boxl4].x » bboxlOKO]; box{4].y ■ bbox£l][l]; box[4].z = bbox{210]; 
box{5]jt * bbox[0][l]; box[5].y a bboxllUl]; box[5].z = bboxt2I0]; 
box(6]jc » bbox[0][l]; box[6].y * bbox[l][l]; boxt6].z = bboxf2Xl]; 
bo3tf7]jc » bboxf0][0]; box[7].y a bboxf box{7].z = bbox[2ll]; 

n_cut a 0; 

printf C\n%3.3f %3.3f %3.3F, planejocate.x, plane Jocate.y,plane_locate.z); 


for (ia0p<6;i++){ 

for (j=0;j<4;j++){ 

k a vertex Jndex[i][j]; 

vtemp.u = boxt-k] jc - planejocatejc; 

vtemp.v = box[k].y - plane Jocate.y; 

vtemp.w * box[k].z - planejocate.z; 

cosa[j] a dot_product(&cut_nonnal,&vtemp); 

) 

if (cosat 0] <=0.0 && cosa( 1] <=0.0 && cosa{2] <=0.0 && cosat3] <=0.0 ) 
goto jump; 


for (j=0*j<4*j++){ 
k = vertex_index[i][j]; 

kl a ( j =s 3 ) ? vertex Jndexti][0] : vertex jndexti][j+l]; 

cutJine.pO jc = boxt~k] jc; 

cutJine.pO.y = boxfkl.y; 

cut_line.p0.z = boxflc].z; 

cut_line.pl jc = boxt-kl]jc; 

cutjine.pl.y = boxtkl].y; 

cutjine.pl.z = boxtklj.z; 


flag = plane_cut_line(cut_nonnal,planejocate,cut_line,&new_point); 
if (flag == PARTCUT1 I I flag == ENDCUTl 

I I flag PARTCUT0 I I flag == ENDCUT0 I I flag =* ONPLANE ){ 
cut_pointtn_cut]jc = new_point.x; 
cut_point{n_cut].y = new_pointy; 
cut_pointtn_cut].z = new_point.z; 
cut_point[n_cut].v_indexO = k; 
cut_pointtn_cut].v_indexl = kl; 
cut_pointfn_cut].flag = UNPICKED; 

if (flag =- ENDCUTl I I flag == ENDCUT0 I I flag =» ONPLANE) 
cut_point(n_cut].e_index = (flag == ENDCUTl) ? kl:k; 

else 

cut_pointjn_cut].e_index = NOTEND; 
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cut_po in t[n_cut++].p Jndex = i; 

if (flag ONPLANE )1 

cutj?oint[n.cut]jc = new_pointJc; 
cut_pointIn_cut].y » new^point-y; 
cut_point[n_cut].z = new_pointz; 
cut_pointln_cut].vjndex0 » k; 
cut_point(n_cut].vJndexl * kl; 
cut_pointln_cut].flag = UNPICKED; 
cut_point{n_cut].eJndex = kl; 
cut j?ointIn_cut++],p Jndex * i; 

} 

if (n_cut>12)( 

printfTn.cut more than 12 !D; 

shut.down.corelO; 

exit(l); 


} 

}/* end of j loop (line loop) *1 

jump: j=0; 

)/* end of i loop (plane loop) *f 


cut_point{0].flag = PICKED; 
tempJndexO = cut.pointfOLvJndexO; 
tempjndexl = cut_point(0].v Jndex 1; 
e .index = cut_point{0].ejndex; 
planejndex = cut.point(0].pjndex; 
new_plane(0].x = cut_point[0].x; 
new_plane[0].y = cut_point[0].y; 
new_plane[0].z = cut_point[0].z; 


if (cut_point[0]. ejndex == NOTEND){ 
for (i=ld<n.cutd++) 

if ((tempJndexO = cut_point[i].vjndex0 
&& tempjndexl =* cut_point!i].vJndexl) I I 
(tempJndexO =« cut_point(i].v Jndex 1 

&& tempjndexl = cut.point[i].vjndex0)) cut_point[i].flag = END; 


else { 

for (i=ld<n_cutd++) 

if (ejndex == cut_point[i].eJndex 

&& plane Jndex != cut _point[i].p Jndex) cut.pointlil.flag = END; 

1 ; 


n.vert = 1; 
pass =0; 

search = BEGIN; 
kk= 1; 

while (search == BEGIN ){ 

if (ejndex != NOTEND ) { 
for (i=l;i<n_cut;i++) 

if (plane Jndex ** cut.point[i].p Jndex 
&& ejndex— cut_point[i].e Jndex) 

cut_point(i].flag = PICKED; 

for (i=14<n_cut;i++) 
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if (plane Jndex »■ cut_point[i].p Jndex 
&& cut _j)oint[i] .flag *■ UNPICKED )( 
cut_point(i].flag ■ PICKED; 
if (e Jndex l* cut_pointfi].e_indexH 
temp JndexO ■ cut __point{i].v JndexO; 
tempjndexl * cut j3oint[i].v Jndex 1; 
e Jndex « cut_point[i].e Jndex; 
new_plane[n_vert]jc * cut_point[i] jc; 
new_plane[n_vert].y - cut_point(i].y; 
new plane[n vert++].z * cut_point(i].z; 
1 

1 

1 

else ( 

M; 

while(i<n_cut){ 

if (plane Jndex == cut_point(i].p Jndex && 
cut_point(i].flag == UNPICKED ){ 
cut_point(i],flag = PICKED; 
temp JndexO = cut_point{i].v JndexO; 
tempjndexl = cut_point{ i].v Jndex 1; 
ejndex = cut_point[i].e Jndex; 
ne w_pl ane[n_vert] jc = cut_point[i] .x; 
new_plane[n_vert].y = cut_point[i].y; 
new_plane[n_vert++].z = cut_pointIi].z; 
break; 

) 

i++; 

} 

} 

jj-i; 

while(jj<n_cut){ 

if (ejndex »» NOTEND ){ 

if(cut_point[jjLflag != PICKED 
&&((temp JndexO == cut_point(jj].v JndexO 
&& tempjndexl == cut_point[jj].vJndexl) 

I 1 (temp JndexO == cutj?oint(jj].vJndexl 
&& tempjndexl == cut_point(j[j].vJndexO))){ 
cut_point[jj].flag = PICKED; 
planejndex = cut_point[jj].pJndex; 
break; 

) 

if (cut_point[jj].flag == END 

&&((tempJndexO == cutj?oint[jj].v JndexO 
&& tempjndexl == cut_point(jj]. v Jndex 1) 

I KtempJndexO = cut^pointfjj]. v Jndex 1 
&& tempjndexl == cut_point[jj].v JndexO ))){ 
search = END; 
break; 


if (ejndex != NOTEND ){ 

if (ejndex == cut_point[jj]. ejndex 

&& cut_point[jj].flag == UNPICKED)! 
cut_point{jj].flag = PICKED; 
planejndex = cut_point[jj].pJndex; 
break; 

} 


46 



if (e_index »» cut_point[jj].e_index 

&& cut_point(jj].flag mm ENDM 

search * END; 

break; 

} 

} 

ij++; 

}/*endofjj loop */ 
if ( pass++ > 6 ) break; 

}/* end of search loop */ 

for (i=0;i<n_vert;i++)( 

new_plane_ptr->x = new_plane[i]jc; 
new_plane_ptr->y = new_plane[i].y; 
new_plane_ptr->z = new_plane[i].z; 
new_plane_ptr++; 

} 


*n_vert_ptr = n_vert; 

}/* end of new_cut_plane */ 


#include <usercore.h> 
iinclude <sun/fbio.h> 
tinclude <stdio.h> 
#include <math.h> 
#include "model.h" 


extern 

extern 

extern 

extern 

extern 

extern 

extern 


POINT cutjplaneG; 

int nvert,npoly,face_remove,num_shadeJevel; 

int npvertQ f *pvertptr[],pvert(]4ndxiistQ, number; 

float planeqQ[4],vertices[][3],normalO[3]; 
short mycolorl,cindex[],pian_infoG,cline[]; 
float xmax^anm,ymax,ymin,zmax f zmin r xcent 
,ycent,zcent,length,emin,emax; 
float dxlistQ.dybstO^dzlistO 

r xhstQ,ylist[],zlistQ > n - vert[][3]Jine_vertO[3]; 


/* line variables */ 

extern int nline,nlvertjinejndexjist{],nplinen; 

/* Externals */ 
extern nobj,npoly,nvert; 

extern free_vert,free_poly,free_obj,free_pvert; 
extern old_vert,old_poly,old_pvert; 
extern float bbox0[2]; 

extern struct image_status_format imagers tat us; 

extern struct obj *objlst; 

extern struct poly *polylst; 

extern struct vert *vertlst; 

extern struct pvert *pvertlst; 
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defin_color() 

I 

int 5 

float x,y,z,x0,y0,z04ength; 

mapjidc_to - worldJ3(-348.,348.,-870.,&x t &y,&z); 
map_ndc_to_world_3(0.0,0.0,0.°>& x 0>&y0,&z0); 
x -*x0; y -»yO; z -*zO; 
length = sqrt(x*x +y*y + z*z); 
if (length != 0.0) 

{ 

x /^lengthy /=length;z /=length; 

} 

for (i=0; i< nvert; i++) 

I 

cindexli]=* 

fabs(normal[i][0]*x+normal[i][l]*y+normal[i][2]*z) * 255.0; 
if (cindex[i]<1.0) cindex(i]=1.0; 
if (cindex{i]> 255.0 ) cindex(i]=255.0; 

I 


} /* end of defin_color *1 

backBurface checkO 

{ 

intijXp; 

int *ptr; 

float x[6] > y[6] f z[6] > nx,ny,nz,c,ai,aj,bi,bj,ci,cj,di,dj;; 
for (p=0;p<npoly;p++) 

( 

ptr=pvertptr{p];k=npvert(p]; 
for (i=04<kd++) 

{ 

j= *ptr++; 

xlist[i]=vertices[j][0]; 
ylistf i]=vertices[j][ 1]; 
zlist{i]=vertices(j][2]; 

} 

for (i=0;i<npvert{p];i++) 

map_world_to_ndc_3(xlist(i] > ylist(i],zlist(i],&xfi] > &y[i] > &z[i]); 

c=0.0; 

ai=x[l]-x[0];aj=ytl]-y[0]; 

bi=x(2]-xil];bj=y(2]-y[l]; 

ci=x[3]-x[2];cj=yf3]-y[2]; 

di=x[0]-x[3];dj=y(0]-yi3]; 

c=ai*bj-aj*bi+bi*cj-bj*ri+ci*dj-cj*di+di*aj-dj*ai; 

if (c <0.0 && plan _info[p]<0 I I c >0.0 && plan Jnfo[p]>0 ) 
plan_info[p]= -plan_info[p]; 

I 

} /* end of backsurface_check */ 
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ndc_to_world() 

{ 

intiik; 

float xy f z; 

for (i*Op<nverti++) 

{ 

n_vert(i][0] *=1020; 
n_vert(i][l] *=1020; 
n_vert[i][2] *=1020; 
if (i**0) 

I 

xmin=xmax=n_vert(i][0]; 
ymm=ymax=n_vert(il 1]; 
zrnin=zrnax=n_vert(i][2]; 

} 

else 

{ 

xmin= (xmin>n_vert[i][0])? n_vert(i][0]:ximn; 
xmax= (xmax<n_vert[i][0])? n_vert[i][0]:xmax; 
yrain= (ymin>n_vert[i][l])? n_vert[i][l]:ymin; 
ymax= (ymax<n_vert(i][l])? n_vert{i][l]:yinax; 
zmin= (zmin>n_vert[i][2])? n_vert[i][2]:zmin; 
zmax= (zmax<n_vert(i][2])? n_vert(i][2]:zmax; 
emin= (emin> cindexti]) ? cindex[i]:emin; 
emax= (emax< cindexli]) ? cindex[i]:emax; 

} 

} 

xcent=xmin+(xmax-xmin)/2; 
ycent=ymin+<ymax-ymin)/2; 
zcent=zmin+(zmax-zinm)/2; 
for (i=0d<nvertd++) 

{ 

if (cindex[i]>l) 

cindex[i]=(int) (cindex[i]-emin)*nuin_shade_level/(emax-emin); 

} 

} /* end of ndc_to_world */ 


drawline() 

{ 

int i j,k,np, color, vjndex; 
int *ptr; 

ptr= &line Jndex_list{0]; 
for (i=0d<nline*4++) 

{ 

np=npiine[i]; 
color = clinefi]; 
if (color < 0) 

ptr +=np; 

else 
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{ 

setjinewidth(0.3); 
setJine_index(color); 
for (j» 03 <npa++) 

{ 

vjndex* *ptr; 

xlistlj] » line_vert{vjndex][0]; 

ylistlj] » Une_vert[v_index][l]; 

zlist4j] = line_vert[y_index][2]; 

ptr++; 

) 

polyline_abs_3(xlist,ylist,zlist,np); 

} 

I 

setJinewidth(O.O); 

} /* end of drawline */ 

draw_axis() 

{ 

int color; 

/♦ 

set_line_index( 100); 

move_abs_3(0.0,0.0,0.0); 

line_rel_3(200.0,0.0,0.0); 

set_line_index(60); 

move_abs_3(0.0,0.0,0.0); 

line_rel_3<0.0, 200.0,0.0); 

setjine_index(80); 

move_absJ3(0.0,0.0,0.0); 

line_reL3(0.0, 0.0,200.0); 

set_linejndex(4); 

*/ 

) /* end of draw axis */ 


defme_color() 

{ 

int i; 

float x,y,z,x0,y0,z0, length; 

map_ndc_to_worldJ3(-348. ,348.-870. ,&x,&y,&z); 
map_ndc_to_world_3(0. 0,0.0, 0.0, &x0,&y0,&z0); 
x -=x0; y -=yO; z -=z0; 
length=sqrt(x*x+y*y+z*z); 
ifGength !=0.0) 
l 

x /=length;y /=length; z /=length; 

} 

for (i=0;i<nvert;i++) 

{ 

cindex[i]=(int)fabs(normal[i][0]^+normal[i][l]^+normal[i][2]*z)* 254.0; 
if (cindex[i]<1.0) cindex[i]=1.0; 
if (cindex[i]>255.0) cindex[i]=255.0; 

1 

} /* end of define_color */ 
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world_to_ndc() 

{ 

intiik; 

for (i»Od<nvertd++) 

map_world_to_ndc_3(verticeB[i][0],vertices[i][l],vertice«[i][2], 

&n_vert(i][0],&n_vert(i][l],&n_vert{il[2]); 

} 

} /* end world_to _ndc */ 


draw_box() 

{ 

inti; 

float xl5]j[5],z[5]; 

set_linewidth(0.3); 

set_lineJndex(BOX_COLOR); 

xi0]=x{4]=x[3]=bboxi0][0]; 

x(l]=x{2]=bbox[0][l]; 

z(0]=ztl]=zl4]=bbori2I0]; 

z[2]=z(3H>boxt2][l]; 

ytO]=yt l]=yC2]=yC3]=yf 4]=bbox{ 1][0]; 

polyline_abs_3(x,y,z,5); 

yl0]=ytl]=y[2]=y(3]=yt4]=bboxt 1][ 1]; 

polyline_abs_3(x,y,z,5); 

y{0]=y(4]*y[3]=bbox[ 1][0]; 

yll]=yt2]=bboxIl][l]; 

zfO)=zf I]=z(4]=bboxl2l0]; 

z{2]=z[3]=bboxI2][l]; 

x[0]=x[ I]=xt2]=xt3]=x[4]=bbox{0][ 1]; 

polyline_abs_3(x,y,z,5); 

x[0]=x( l]=x[2]=x{3]=x(4]=bbox[0][l]; 

polyline_abs_3(x,y,z,5); 


x[0]=x[4]=x[3]=bbox(0][0]; 

x[l]=x{2]=bboxt0][l]; 

y[0]=ytl]=yl4]=bbox[l][0]; 

yt2]=y(3]=bbox[l][0]; 

yl2]=yt3]=bbox[l][l]; 

z[ 0]=z { l]=z£21=z£3j=z£4j=bbox£2][0]; 

polyline_abs_3(x,y,z,5); 

z[0]=z[ l]=z(2]=z(3]=z(4}=bbox{2][ 1]; 

polyline_abs_3(x,y,z,5); 

draw_axis(); 

} /* end of draw_box */ 

drawface(p) 

intp; 

( 


intijjc; 
int *ptr; 

float x(6] ) y[6],z(6],nx,ny > rLZ,c,ai,aj,bi ) bj,ci,cj ) di,dj;; 
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ptr*pvertptr[p]; 

for (i«On<npvertf p];i++) 

{ 

j» *ptr++; 

xlist[i]*vertices(jIO]; 

yliatli]avertice«ljll]; 

ziistti]=vertices(jI2]; 

1 

for (i»0*4<npvert(p];i++) 

map_worid_to_ndcj3(xlist[i] > ylistIi] > zlistIi] f ixli],4y[i]*«*lii); 

c=0.0; 

ai»x[l]-xIO];aj»y[ l]-y[0]; 
bi=*x[2]-x[l];bj=y(2]-ytl]; 
ci=x£3]-x£2];cj=>i3]-y[2]; 
di=x[0]-x[3];dj=y[0]-yf3]; 

c=ai*bj-aj*bi+bi*cj-bj ,(, ci+ci*dj-cj*di+di*aj-dj*ai; 

xliatl4}=xlist(0]ylistt4]=ylist[0]* > 2list[4]=zlisU0]; 
if (c <0.0 ) 

polyline_abs_3(xlist,ylist,zlist f npvert[p]+l); 


} I* end of drawface */ 


drawobjO 

{ 

inti; 

float xf5l,y[5],z[5]; 
set_line_index(LINE_COLOR); 
if (image_status. special = TRUE) 

{for (i=618*4<npolyd++) 

drawface{i); 

} 


else 

for (i=0;i<npoly4++) 
drawface(i); 


draw_axis(); 

} /* end of drawobj *f 

draw_wir efr am e( ) 

{ 

short *ptr; 

struct vert *vertlstptr; 

int i j,k4ndex,num_vert,polyJndex; 


set_line_index(LINE_COLOR); 
for (i=04<nobj4++){ 

poly Jndex = objlstjil.index; 
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for (j= 03 <objlst(i].npol 5 rj-H-){ 

linking(&polyJndex); 
ptr * polylsttpolyJndex].pvert.ptr, 
num_vert ■ polylst[poly_index++].npvert; 
for (k»0 Jt<num - vert;k++M 

vertlstptr = vertlst + *ptr++; 
xlistfk] 3 vertlstptr->vertexjc; 
ylistfk] = vertlstptr->vertex.y; 
zlistfk] = vertlstptr->vertex.z; 


} 

xlistfk] =» xlist(0];ylistfk] = yiistJO]; zlist[k] = zlistfO]; 

polyline.abs.Sfxlistjlist^zlis^num^vert+l); 

}/* end of j poly loop */ 

}/* end of i obj loop *f 

} /♦ end of draw_wirefram */ 


drawobjlO 

{ 

int i j.kAndex.num^vert^poly Jndex; 
float x*y,z^tO,yO,zO, length; 
short *ptr,*colorlst, color Jevel; 
struct vert *vertlstptr; 

map_ndc_to_world_3(-348. 0,348.0, *870.0, &x,&y,&z); 
map_ndc_tojworld_3(0.0,0.0,0.0,&x0,&y0,&z0); 
x -= x0;y -= yO; z -= zO; 
length * sqrt(x*x+y*y+z*z); 
if (length != 0.0) 

{ 

x /=length;y /^length; z /^length; 

1 


colorist = (short *)malloc(free_vert * sizeoflshort)); 
for (i=0;i<free_vert;i++){ 
colorlst{i]= fabs 

(vertlst[i].normaljc*x+vertlst[i].nonnal.y^y 

+vertistii].normal.z*z)*1000.; 

if (i==0) 

emin=emax=colorlst[0]; 

else 


{ 

emin= (emin> colorlstf i]) ? colorlst(i]:emin; 
emax= (emax< colorlsUi]) ? colorlst(i]:emax; 
1 


for (i=0;i<free_vert;i++) 

colorlstfi] = (colorlst[i]>l) ? 

((intXcoiorlst[i]-emin)*num_shade_level/(emax-emin)):l; 


for (i=0;i<nobj;i++){ 

poly_index = objlst(i].index; 

for (j=0*j<objlst[i].npoly'j++){ 

/♦ 

again3: while (polylst[polyJndex].info = DELETE) poly_index++; 
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if (polylstf polyJndexLinfo < 0 ) 

polyjndex = polylstIpoly_index].pvert.index; 
if (polylst(polyjndex].info =* DELETE) goto again3; 

*! 

linking( &poly Jndex); 

ptr ■ polylsttpolyjndexj.pvert.ptr; 

num_vert = polylst[polyJndex].npvert; 

color Jevel ■ (polylst[poly - index++].info-l)*num_shadejevel+l; 
for (k=0;k<num_vert;k++){ 

indxlisttkl^colorlstl^ptrj+color Jevcl ; 
vertlstptr = vertlst + *ptr++; 
xlistfk] = vertlstptr->vertexjc; 
ylistfk] = vertistptr->vertex.y; 
zlistfk] = vertlstptr->vertex.z; 

I 

set_vertexjndices(indxlist,num_vert); 

polygon_abs_3(xlist,ylist,zlist,num_vert); 

}/* end of j poly loop */ 

}/* end of i obj loop */ 

free(colorlst); 

} /* end of drawobjl *t 


image switchO 

{ 

if (face_remove) 

{ 

new_frame(); 
face_remove= FALSE; 


else 

{ 

if (image_status.wire_frame) 

{ 

delete_retained_segmentf WIREFRAME); 
create_retained_segment(WIREFRAME); 
draw_wireframe(); 

close_retained segment WIRE FRAME); 

) 

if (image_status.out_box) 

{ 

delete_retained_segment( OUTBOX); 
create_retained_segment(OUTBOX); 
draw_box(); 

close_retained_segment(OUTBOX); 

} 

if (image_status.line) 

{ 

delete_retained_segment(DRAWLINE); 
create_retained_segment( DRAWLINE); 
drawlineO; 

close_retained_segment(DRAWLINE); 

} 

1 


) 
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draw_cut_plane(n_point) int n_point; 

(inti; 

move_abs_3(cut_plane[0].x,cut_plane[0].y,cut_plane[0].z); 

for (i*0;i<n_point;i++M 

xlistfi] * cut_plane[i]jc; 
ylistfi] = cut_plane[i].y; 
ziistfi] = cut_plane[i].z; 


xlist(i] = cut^plane[0]jc^Ust(i] = cut_plane(0].y;zlist[i] = cut_planefO].z; 
set_lineJndex(MENU_TEXT_COLOR); 

poiyline.abs.SCxlistjlist^list^^point+l); 

} /* end of draw_cut_plane */ 


#include <math.h> 

#include H model.h" 

/* Cross-product, dot_product and vector normalization */ 


/* ************ Dot-product ******************/ 
/* Input the points of vector A and B */ 

/* Return value : cos( theta) 


float dot_product(pva,pvb) VECTOR *pva,*pvb; 

{ 

float cos-theta; 

normalize(pva); 

normalizefpvb); 

cos_theta=((pva->u)*(pvb->u) + (pva->v)*(pvb->v) + (pva->w) (pvb->w)); 
retum(cos-theta); 

} /* end of dot_product */ 

,*************** Normalize***"*************/ 

/* Input the point of the vector which will be normalized */ 

normalize(pv) VECTOR *pv; 

{ 

float m ’ / x*/ xx 

m= sqrt((pv->u)*(pv->u) + (pv->v)*(pv->v) + (pv->w)*(pv->w)); 
pv->u / =m; 
pv->v /=m; 
pv->w /=m; 

} /* end of vector_normalize */ 

/************ Cross.product **************/ 
cross_product(pvc,pva,pvb) 

VECTOR *pva,*pvb,*pvc; 

{ 

/♦normal ize(pva); 
normalize(pvb);*/ 
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pvc->u ■ (pva->v) * (pvb->w) - (pva->w) * (pvb->v); 
pvc->v m (pva->w) * (pvb->u) - (pva->u) * (pvb->w); 
pvc->w - (pva->u) * (pvb->v) - (pva->v) * (pvb->u); 

ifl(pvc->u**0) && (pvc->v*»0) && (pvc->w**0)) 

{ 

if (pvb->u*»0) {pvc->u * 1.0; pvc->v=0.0; pvc->w*0.0;} 

else 

l 

if ( pvb->v*»0) { pvc->u a -1.0/pvb->w; pvc->v»0.0; pvc->w « 1.0/pvb->u;l 

else 

{ 

pvc->u = 1.0; 
pvc->w » 1.0; 

pvc->v » -( pvb->u + pvb->w) / pvb->v; 

} 

} 

normalize(pvc); 

} 


*/ 

} /* end of cross_product */ 


matinvfptr) 

float *ptr; 


short index[3][2], ipivot(3]; 
float pivot(3],mtxout(3][3]; 
short row, colum; 
float max; 
short i,j, k, 1; 


for (i=0;i<3;i++) 

for (j=0;j<3;j++) 

mtxoutfilj] = *(ptr+i*3+j); 

for (j = 0; j < 3; j++) 
ipivotfj] = 0; 


for (i = 0; i < 3; i++) { 
max = 0.0; 

/* Search for pivot element */ 
for (j = 0; j < 3; j++) { 

if ( ipivot(j] = 1) 
continue; 

for Gt a 0; k < 3; k++) { 

if(ipivot[k] as 1) 
continue; 
if (ipivotfk] > 1) 

retum( 0); 

if (fabs(max) < fabs(mtxout(j][k])) { 
row = j; 
colum = k; 
max a mtxoutfjllk]; 
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} 


J 

} 

f* Row intercahnge */ 
ipivottcolum] +* 1; 
if (row !» colum) { 

for (1 =» 0; l < 3; 1++) { 

max = mtxoutjrow][l]; 
mtxout[row][ll = mtxouttcolumlOl; 
mtxout[colum][l] = max; 


} 

indexii][0] = row; 

indexfi][l] * colum; 

pivotfi] » mtxout(colum][ colum]; 

mtxout(colum][colum] = 1.0; 
for (1 = 0; 1 < 3; 1++) 

mtxout[colum][l] /= pivotfi]; 


for (j =0; j < 3; j++) 

if (j 1= colum) { 

max = mtxoutfjllcolum]; 
mtxout{j][columj « 0.0; 
for (1=0; 1<3; 1++) 

mtxout(j][l] -= mtxoutfcolumlQ] * max 


} 


} 


for (i = 0; i < 3; i++) { 

1 = 3 • 1 - i; 

if (indexfl][0] != index[l][l]) ( 
row = index[l][0]; 
colum = index[l][l]; 
for (k=0; k<3; k++) { 

max = mtxout(k][row]; 
mtxout[k][row] = mtxout[k][colum]; 
mtxout[k][colum] = max; 

} 


} 


) 

for (i=0;i<3*4++) 

for (j=0;j<3;j++) 

*(ptr+i*3+j) = mtxout[i][j]; 


/* 

* Module: load_file() 

* Functions: 

1. Check data format; 

2. Allocate memory; 

3. Setup scaling parameters; 
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4. Genarate object list; 

5. Compute the parameters of each poly. 

6. Compute the average normal at each vertex; 

7. Free the temporary storage. 


* Input: 

data hie: objects data, polygons data. 
♦Output: 

object list, polygons list, vertices list; 

* Last change 10/8/88 
♦/ 


#include <math.h> 
#include <stdio.h> 
#include "model.h” 


extern int 
extern int 
extern float 
extern short 
extern float 

extern float 


npoly,face_remove,num_shadeJevel; 
npvertQ,*pvertptrQ,pvertG,indxl is tQ number; 
bboxOI2],planeqQ[4],verticesn[3]jiormal[]I3]; 
mycolorl,cindex[],plan_infoG,clineD; 
xmax,xmin,ymax,ymin,zmax,zmin t xcent 

,ycent,zcent, length, emin,emax, scale; 
dxlistO,dyiistO,dzlistG 

,xhstGjlistO,zhstQ,n_vertO[3],hnejfertQ3]; 


/* line variables *! 

extern int nline,nlvert,linejndexjist0^ipline0; 


extern struct image_status_format image_status; 


/* Externals */ 

extern nobj,npoly,nvert; 

extern free_vert,free_poly,free_obj,free_pvert; 

extern old_vert,old_poly,old_pvert; 

extern float bbox[3][2]; 

/* Externals */ 

extern struct obj *objlst; 
extern struct poly *polylst; 
extern struct vert *vertlst; 
extern struct pvert ♦pvertlst; 

load_file(filename) 
char *filename; 

{ 

int i j; 

int vtmp,vl,v2,v3; 
float ftmp,maxd,offset{3]; 
float x,y,z,x0,y0, length; 

FILE *fptr; 
struct count! 

short n 
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}; 

struct count ^normal count; 

/* 1. Check data format; •/ 

if ((fptr * fopen(filename, V)) »* NULL) { 

printfCCant open file: %s\n", filename); 
retum(l); 

l 

fscanflfptr, ’'%d%d%d m , &nobj,&npoiy,&nvert); 
if (nobj > MAXOBJ ) { 

printfl"Too many objects \n *); 
retum(l); 

} 

if (npoly > MAXPOLY ) { 

printfTToo many polygonsXn"); 
retum(l); 

} 

if (nvert > MAXVERT ) { 

printfTToo many verices\n"); 
retum(l); 

} 

if (nobj<=0 I I npoly <=0 I I nvert <=2 )( 
printfTCheck your data! \n M ); 
retum(l); 

} 

/* 2. Allocate memory *f 

if (objlst !=s NULL ) free(objlst); 

objlst = (struct obj *)malloc(nobj * sizeofl struct obj)); 

if (polylst != NULL ) free( polylst); 

polylst = (struct poly *)malloc((npoly+MAXNEWPOLY) * sizeoffstruct poly)); 
if(vertist != NULL ) free(vertlst); 

vertlst = (struct vert *)malloc((nvert+MAXNEWVERT) * sizeoftstruct vert)); 
if (pvertlst != NULL ) free(pvertlst); 

pvertlst ={struct pvert *) malloc((MAXPVERT + MAXNEWPVERT) * sizeoflshort)); 
/* 3. Setup scaling parameters */ 

fscanflfptr, &bboxf0][0], <StbboxfO][l], &bbox(l][0],&bbox[ 1][1], 

&bbox[2][0], &bbox[2][l]); 

maxd = 0.0; 

for (i = 0; i < 3; i++) { 

offsetfi] = (bbox{i][0] + bbox£i][l]) / 2.0; 
bbox{i][0] offset[i]; 
bboxfilf 1] -= offsetli]; 
if (bbox£i][0] > bbox[i][l]) ( 
ftmp = bboxfi][0]; 
bbox[i][0] = bboxfi][l]; 
bbox[i][l] = ftmp; 

} 

if (maxd < bbox£i][l]) 

maxd a bboxfill]; 

) 

scale « 1000.0 / maxd; 
for (i a 0; i < 3; i++) { 
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bboxtilO] *= scale; 
bbox[i][l] *= scale; 
) 


free_vert=free_pvert=free_poly»free_pvert=0; 

/* 4. Genarate object list */ 

if (normalcount != NULL ) free(normalcount); 

normalcount* (struct count *)malloc((nvert )* sizeoftstruct count)); 
if (normalcount =■ NULL ) { 
printfTNnMerory too small!! I n M ); 
exit(l); 

} 

for (free obj»0;free_obj<nobj;free_obj++) /* objects input *1 

T 

fscanflfptr, "9kl%d^objlst[free_obj].npoly, &objlstffree_obj].nvert); 
objlst[free_obj].index= free_poly; 

/* if (normalcount != NULL ) free( normal count); 

normalcount* (struct count *)malloc(objlst[free_obj]. avert 
* sizeoflstruct count));*/ 

f* vertices input */ 

for (i = 0; i <objlst(free_obj].nvert; i++) 

{ 

fscanflfptr, H W %f\ &x, &y,&z); 
vertlst[free_vert].vertexjc = (x-offset{0])*scale; 
vertlst[free_vert].vertex.y = (y-offset[l])*scale; 
vertlst{free_vert].Yertex.z = (z-offset[2])*scale; 
vertlstffree_vert].normal.x = vertlst[free_vert].normal.y 

= vertlst{free_vert].normal.z =0.0; 
normalcount(free_vert++].n = 0; 

}/* end of vertices input */ 

/* polygons input */ 

for (i = 0; i <objlst{free_obj].npoly ; i++) 

( 

fscanflfptr, '%d’\&vtmp ); 
poly lst(free_poly].np vert = vtmp; 

if ((old_pvert + polylst[free_poly].npvert) > MAXPVERT) 

{ 

printfTCheck your data !! RUN OVER PVERTLISTAn"); 
retum(3); 

} 

fscanflfptr, "%d”, &vtmp); 
polylst[free_poly].info=vtmp; 

poly lst(free_poly].p vert .ptr = &pvertlst(free_pvert].index; 
for (j = 0; j <polyist(free_poly].npvert ; j++) 

{ 

fscanflfptr, ” %d'\ &vtmp); 
pvertlst(free_pvert++] .index = vtmp - 1; 

} 

polylst[free_poly].planq.a = polylst[free_poly].planq.b 
= polylst[free_poly].planq.c = polylst(free_poly].planq.d = 0.0; 

vl = pvertlst{free_pvert-l]. index; 
v2 m pvertlst{free_p vert-2]. index; 
v3 = pvertlst{free_p vert-3], index; 


Q-3. 



/* 5. Compute the parameters A,B»C and D of each poly. */ 

for(j«0;j<3;j++) 

{ 

polylst(free_poly].planq.a +■ vertlst(vl].vertex-y* 

(vertlst[v2].vertex^ - vertlst[v3].vertex.z); 
polylsttfree_poly].planq.b +» vertLst[vl].vertexjc * 

(vertlst(v3].vertex.z - vertlst[v2].vertex.z); 
polylstf free_poly].planq.c +* vertlst(yl].vertex.x * 

(vertlst{v2].vertex.y - vertlst[v3].vertex.y); 
polylst(free_poly].planq.d += vertlstfvlJ.vertexjc * 

((vertlst[v3].vertex.y * vertlst(v2].vertex.z) - 
(vertlst[v2].vertex.y * vertlst(v3].vertex.z)); 
vtmp = vl; vl = v2; v2 = v3; v3 = vtmp; 

) 

/* 6. Compute the average normal at each vertex */ 

x * polylst(free_poly].planq.a; 
y = polylst{free_poly].planq.b; 
z = polylst[free_poly].planq.c; 
length = sqrt( x*x + y*y + z*z); 

for(j = 1; j <=* polylst[free_poly].npvert 3 ++) 

{ /* accum normls */ 

vtmp = pvertlst(free_pvert-j]. index; 
vertlst(vtmp].normaljc += x/iength; 
vertlst{vtmp].normal.y += y/length; 
vertlst[vtmp].normal.z += z/length; 
normaicount[vtmp].n++; 

} 


free_poly++; 


} /* end of poly input */ 

j=0; 

for (i = old_vert; i < free_vert; i++) 

{ 

vertlst[i].normaljc/= normalcount{j].n; 
vertlst{i].normaLy /= normal count(j].n; 
vertlst[i].normal.z/= normal countfj++].n; 
} 

/* 7. Free the temporary storage */ 


old_poly=free_poly; 

old_pvert=free_pvert; 

old_vert=free_vert; 

} /* end of object input */ 

f* 7. Free the temporary storage */ 
free(normalcount); 

fclose(fptr); 
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retum(NULL); 


} /* end of load JileO */ 


tin elude <usercore.h> 

#include <sun/fbio.h> 
finclude <stdio.h> 
findude <math.h> 

#include "model.h" 

extern int n.button; 

extern float menu_x,menu_y > menu_h > menu_w,menu_f > window_fact; 
extern struct menu_table_format menu_table[]; 

extern struct view_parameters_format my_view_parameter«; 

extern struct image_status_format image_status; 


int menu_select() 

{ 

float mx,my,fx,fy,fz; 
float x,y,z; 

int ijXbutton^nycase; 

set__primiti ve - attributes(&PRIMATTS) ; 

mycase * 0; 

set_echo( LOCATOR, 1 ,1); 
for (;;) 

{ 

do 

a wait_any_button_get_locator„2(2 000000 Ojl^&buttonj&mXj&my); 

while(button == 0); 

if ((button==l)&&(mx>0.9)&&(my>0.01)&&(mx<0.98)&&(my<0.66)) 

{ 

mx={mx-.9)*1000.0; 
my=my*767.0/0.75; 
for (k^;k<n J>utton;k-H-) 

if ((mx<menu_tabletk].maxx)&&(mx>menu_table[k].minx)&& 

(my<menu_table[k].maxy)&&(my>menu_table[k].miny)) 

mycase=menu_table{k].button_id; 

break; 

} 

1 

set_echo(LOCATOR,l,0); 

retum(mycase); 

} /* end of menu_select */ 

int call_menu(raenuname,flag) 
int menuname,flag; 

{ 
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inti; 

inquire_viewing_parameters(&my_.view_parameters); 

setjmenu_vw<); 

set jwgment_viaibility(menunaine, TRUE); 
i * menu_selectO; 

setjBegment_visibility(menunaine,flag); 

get_viawing_parameter8(&my_view_parameter8); 

retumti); 

} /* end of call menu */ 


set_menu_vw() 

{ 


set_view_reference_point(0.0,0.0,0.0); 

set_view_plane_normal(0.0,0.0, -1.0); 
set_view_plane_distance(0.0); 
set_projection(PARALLEL,0.0,0.0,1.0); 
setjvi e w_up_3(0. 0, 1 .0 , 0.0); 
set_window(0.0, 100.0,0.0,767 .0); 
set_view_depth(0.0,1.0); 
set_window_clipping(FALSE); 
set_viewport_3(0.92,1.0,0.0,.75,0.0,1.0); 
} /* end of set_menu_view */ 


init_menu() 

{ 

intij,k; 

inquire _viewing_parameters(&my_view_parameters); 
set_menu_vw(); 

create_retained_s egm en K MENU) ; 
set_segmentjvasibility( MENU, 1 TRUE); 

build_menu_table(); 

buiid_menul(); 

/* 

build_menu(); 

*/ 

close_retained_segment(MENU); 

set_viewing_parameters(&my_view_pararneters); 

} /* end of init_menu */ 

buil d_m enu_tabl e{ ) 

{ 

FILE *fptr; 
short i,nj; 
int status; 
char nam[16]; 
float px,py; 

fptr = fopenCmenuJHe", V); 

fscanflfptr, , &menu_x,&menu_y > &menu_w,&menu_h,&n_button,&menu_f); 

px= menu_x; 


63 



py- menu_y; 

for (i«0?<n Jjutton^ ++) 

{ 

for (j*03<163++) nam[j]=4; 

fscanflfptr, "%d\ fcstatus); 

fscanflfctr, ”%d m , &menu_table[i].button Jd); 

fscanflfrtr, "%d\ &menu Jable[i].button_frame); 


&>> 

do 

fscanflfptr, m %c" 9 &nam[j]); 
while (nam[j++]!=10); 


menu_table{i].button_status = status; 
for (j»0;j<16;j++) 

if ((namtj] >=65 && namlj] <=90 ) I I (nam[j] >=97 && nam(j] <=122)) 
menu_table[i].name[j]=nainlj]; 

else 

menu_table[il.namelj]=32; 

if (menu_table[i].button_s tatus = TRUE) 

menu_table[i].maxx=px+menu_w-menu_h/menu_f; 

else 

menu_table[i].maxx=px+menu_w/2.0 -menu_h/menu_f; 

menu_table[i].minx=px+menu_h/menu_f; 

menu_table[i].maxy=py; 

menu - table[i].miny=py-menu_h; 

if (menu_table[i].button_s tatus == TRUE) 
py -= menu_h; 

else 

{ 

if ((px+menu_w/2.0)>= (menu_x+menu_w)) 

{ 

px = menu_x; 
py .= menu_h‘» 

} 

else 

px +=menu_w/2.0; 

) 

}/* end of i loop *! 
fclose(fptr); 

} /* end of build_raenu_table *f 

build_menul() 

{ 

intij^k; 

float px,py, factored 12],y[12],z[12]; 


setjinewidth(0.3); 

set_lineJndex(MENU_BOX_COLOR); 



factor*menu_h/menu_f; 

px=menu_x; 

py»menu_y; 

for (i-0?<n_buttonp++) 

( 

x[ 0 ]*menu_tableCi].minx-factor; 

ytO]»menu_table[i].miny+factor, 

xtl]=menu_table(i].ininx; 

yt l]=menu_table{i].ininy; 

xt2]*menu_table(i].maxx; 

y[2]=menu_table{i].ininy; 

x[3]=x[2]+factor; 

yt3]-yC0]; 

xt4]=x£3]; 

yt4]=menu_table[i].maxy-factor; 

x[5]»menu_table(i].maxx; 

y{5]*menu_table[i].maxy; 

x£6]=x£l]; 

y[6]=menu_table[i].maxy; 

xf7]=x£0]; 

yC7]=yf4]; 

x[8]=xt0]; 

yt8]=y[0]; 

move_abs_3(x[7],y(7],0.5); 
if (menu_table{i].button_frame == 0) 

( 

polyline_abs_2(x,y,9); 

move_abs_2(x[0],y[0]+factor*1.8); 

set_text_index(MENU_TEXT_COLOR); 

text(&menu_table[i].name[0]); 

} 

else 

{ 

set_fill_index(menu_tabl e[ i ] .bu tton_frame); 
polygon_abs_2(x,y ,9 ); 
move_abs_2(x[0],y[0]+factor*1.6); 
set_text_index( 129); 
text(&menu_table[il.name[0]); 

} 

if (menu_table[i].button_status== TRUE) 

1 

py -=menu_h; 
px = menu_x; 

1 

else 

( 

if ((px+menu_w/2.0)>= (menu_x+menu_w)) 

l 

px = menu_x; 
py -= menu_h; 

1 

else 

px +=menu_w/2.0; 

) 

1 

set_linewidth(0.0); 

} /* end of build_menu *1 
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#in elude <usercore.h> 

#in elude <eun/fbioii> 

#include <*tdio.h> 
finclude <math.h> 

#indude "model.h" 

extern float redO, grna,bluadot,T3,L3, window Jact; 
extern struct vwsurf *our_surface; 

shut_down_corelO 

{ 

terminate_device(KEYBOARD,l); 
deselect_view_surface(our_surface); 
terminate_view_surface(our_surface); 
termin ate_core( ); 

} /* end of shut_down_corel */ 


start_up_coreO 

l 

inti; 

FILE *fptr, 

initialize_core{BASIC, SYNCHRONOUS, THREED); 
our_surface*>cinapsize = 256; 

our_surface->cmapname[0]= , \0 , ; 

iflinitialize_view_surface(our_surface,TRUE)) exitf 1); 
initialize_device(BUTTON , 1); 
initialize_device(BUTTON ,2); 
initialize_device<BUTTON ,3); 
select_view_8urface(our_surface); 
set Jight_direction( -0.45, 0.45, -0.45); 

set - shading_parameters(.5,.5,.5,0.5,7.,0,0); 

initialize_device( KEYBOARD, 1); 

set_echo_surface(KEYBOARD,l,our_surface); 

set_keyboard(l,80/'\l); 

initialize _device{LOCATOR, 1 ); 

initialize_deviee{ PICK, l)j 

set_pick(l, 0.001); 

set_echo(LOC ATOR, 1 ,0); 

set_echo_surface(LOCATOR,l,our_surface); 

/*set_output_clipping(TRUE);*/ 

fptr = fopenCnew_color.dat", V); 

inquire_color_indices(our_surface,0,255,red,grn,blu); 

fscanflfptr, &dot,&T,&B,&L,&R); 

for (i=0d<256p++) 

fscanflfptr, "%f\ &red[i]); 
for (i=0d<256d++) 

fscanflfptr, M %f\ &gm[i]); 


66 



for(i*0;i<256;i++) 
fscanflffctr, "%T 9 &blu£ij>; 

define_colorJndices(our_8urface, 0,255, red, gm,blu); 
fclo»e(fptr); 

} /* end of start_up_core ♦/ 


setvwpo(vx, vy, vz,bbox) 
float vx, vy, vz,bbox{3][2]; 

{ 

inti; 

float diag, del, objdist, near; 

set_view_reference_point(vx, vy, vz); 
setjview _plane_normal(-vx, -vy, -vz); 
set_projection(PERSPECTTVE, 0., 0., 0.); 
set_view_plane_distance(256.0); 
if ((vx =» 0.0) && (vz =» 0.0)) 

set_view_up_3(0.0, 0.0, vy); 

else 

setjview_up_3(0.Q, 1.0, 0.0); 

set_window(-80.0*window_fact, 80.0*window_fact, -65.0*window_fact, 65.0*window_fact); 

diag a 0.0; 

for(i = 0; i < 3; i++) { 

del a bbox(i][l] - bbox[i][0]; 
diag +* del * del; 

} 

diag = sqrt(diag) / 2.0; 

objdist = sqrt( vx*vx + yy*vy + vz*vz); 

near a (diag >a objdist) ? objdist;2.0 : objdist-diag; 

set_view_depth( near, objdist + diag); 

set_window_ciipping(THUE); 

setJront_planejdipping(TRUE); 

set_back_piane clipping(TRlJE); 

/* 

set_viewportJ3(.17, .92, 0., .75, 0.0, 1.0); 

V 

}f* end of setvwpo */ 

setvwpvl) 

{ 

set_view_reference_point(0.0,0. 0 , 0 . 0 ); 

set_view_plane_normai( 0 . 0 , 0 . 0 ,- 1 . 0 ); 

set_view_plane_distance(0.0); 

set_projectiomPARALLEL, 0 . 0 , 0 . 0 , 1 . 0 ); 

setjview_up_3(0.0,1.0,0.0); 

set_windowt0.0, 1023.0, 0.0, 767.0); 

setjview_depth(0.0,1.0); 

set_window_clipping(FALSE); 

set__viewportJ3(Q.0,1.0,0.0,.75,0.0,1.0); 

} /* end of setvwpv */ 
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